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APRESENTAÇÃO 


A investigação científica está ligada ao uso de análises técnicas com a finalidade 
de encontrar respostas para determinados questionamentos. Nas ciências agrárias não é 
diferente, e se torna cada vez mais importante, tendo em vista que as novas tecnologias 
são obtidas a partir de estudos que visam melhorar técnicas e, até mesmo, acabam por 
inserir formas de execuções inovadoras para alguns processos agrícolas. 

Com o constante crescimento da população mundial, o setor agrícola tende a 
necessitar de meios mais eficazes de produção para suprir as demandas alimentícias 
mundiais futuras, e estas exigências acabam por gerar um conjunto de questionamentos 
que só podem ser desvendados através de investigações precisas. 

O grande desafio da agropecuária mundial hoje é produzir mais e melhor, ocupando 
menos espaço de forma sustentável, e para isso há muitos anos pesquisas vêm sendo 
realizadas com a finalidade de contribuir para melhorias das ações, proporcionando 
ambientes de produção equilibrados e que permanecem em constantes melhorias. 

Portanto, é notória a importância dos questionamentos gerados no processo de 
investigação e mais importantes ainda são as respostas resultadas através dele, que 
acabam por resultar em soluções inovadoras para substituição total ou parcial dos métodos 
confrontados. 

Neste sentido, a presente obra reúne pesquisas inovadoras para a difusão de ideias 
importantes e com impacto direto no setor em questão, visando entregar informações de 
alto valor e relevância para o leitor e atualizando-o das tecnologias e inovações que são 


cada vez mais comuns do setor agrário internacional. 
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RESUMEN: Se seleccionaron de forma 
divergente dos líneas de conejos para aumentar 
y disminuir la variabilidad ambiental del tamahio 
de la camada al nacimiento. La disminución de 


Investigación, tecnología e innovación en ciencias agrícolas 


CONEJOS 


la variabilidad del tamafio de la camada genera 
hembras más resilientes, con menos sensibilidad 
al estrés y las enfermedades, siendo un criterio 
de selección útil para mejorar la sensibilidad 
ambiental. Los ácidos grasos modulan la función 
de las células inmunitarias. Los ácidos grasos 
saturados (SFA) tienen un efecto inhibidor sobre 
la proliferación de linfocitos, los ácidos grasos 
monoinsaturados (MUFA) ejercen un efecto 
protector y antiinflamatorio sobre los macrófagos, 
y los ácidos grasos poliinsaturados n-3 (PUFA) 
afectan la respuesta de los linfocitos por medio 
de las citoquinas IL-1, IL-2, IL-6 y TNF, así como 
de la prostaglandina E2 y leucotrieno B4. El perfil 
de ácidos grasos plasmáticos se evaluó en 10 
hembras de la línea homogénea y 12 hembras 
de la línea heterogénea de la 12º generación de 
selección. La línea homogénea mostró niveles 
más altos de SFA (+3,98 ng / ml P = 0,90 para 
C14:0; +2,30 ng / ml P = 0,98 para C15:0; +54 
ng/ml P = 0,90 para C16:0 y +29 ng / ml P = 
0,90 para C18:0) y de MUFA (+12,0 ng / ml P 
= 0,98 para C16:1 y +53 ng / ml P = 0,90 para 
C18:1n9c) que la línea heterogénea. Además, 
esta línea también tenía una mayor cantidad de 
PUFA n-3 (+2,18 ng / ml P = 0,90 para C18:3n3 
y +1,91 ng / ml P = 0,90 para C20:5n3) y una 
menor cantidad de PUFA n-6 (-3,66 ng /ml P = 
0,96 para C20:3n6 y -0,28 ng / ml P = 0,90 para 
C20:4n6) que la heterogénea. En conclusión, la 
selección por sensibilidad ambiental muestra una 
respuesta correlacionada en el perfil de ácidos 
grasos plasmáticos. 

PALABRAS CLAVE: MUFA, PUFA n-3, PUFA 
n-6, linfocitos, resiliencia, sensibilidad ambiental. 
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ABSTRACT: Two rabbit lines were divergently selected for increasing and decreasing 
environmental variability of litter size at birth. Decreasing litter size variability generates more 
resilient females with less sensitivity to stress and diseases, being a useful selection criterion 
to improve environmental sensitivity. Fatty acids modulate the immune cell function. Saturated 
fatty acids (SFAs) have an inhibitory effect on lIymphocyte proliferation, monounsaturated fatty 
acids (MUFASs) exerts a protective and anti-inflammatory effect on macrophages, and n-3 
polyunsaturated fatty acids (PUFASs) affect response of lymphocytes by mean of IL-1, IL-2, 
IL-6, TNF as well as prostaglandin E2 and leukotriene B4. Plasma fatty acids profile was 
assessed in 10 females from the homogeneous line and 12 females from the heterogeneous 
line ffom the 12th generation of selection. The homogeneous line showed higher levels of 
SFAs (+3.98 ng/ml P=0.90 for C14:0; +2.30 ng/ml P=0.98 for C15:0; +54 ng/ml P=0.90 for 
C16:0 and +29 ng/ml P=0.90 for C18:0) and MUFAs (+12.0 ng/ml P=0.98 for C16:1 and +53 
ng/ml P=0.90 for C18:1n9c) than the heterogeneous line. Besides, this line had also a higher 
amount of n-3 PUFAs (+2.18 ng/ml P=0.90 for C18:3n3 and +1.91 ng/ml P=0.90 for C20:5n3) 
and a lower amount of n-6 PUFASs (-3.66 ng/ml P=0.96 for C20:3n6 and -0.28 ng/ml P=0.90 for 
C20:4n6) than the heterogeneous one. In conclusion, selection for environmental sensitivity 
shows a correlated response in the plasma fatty acids profile. 

KEYWORDS: MUFA, PUFA n-3, PUFA n-6, lymphocytes, resilience, environmental sensitivity. 


INTRODUCCIÓN 


Se ha realizado con éxito un experimento de selección divergente por sensibilidad 
ambiental en conejos en la Universidad Miguel Hernández de Elche (Blasco et al., 2017). 
Este caracter se midió como la variabilidad ambiental del tamafio de la camada al parto. 
La selección para reducir la variabilidad del tamafio de la camada puede ser una forma útil 
de mejorar la resiliencia de la hembra, que se define como la capacidad de las hembras 
de verse mínimamente afectadas por las perturbaciones, o su capacidad para volver 
rápidamente al estado anterior a la exposición a una perturbación (Colditz y Hine, 2016). 
En este sentido, la línea seleccionada para reducir la variabilidad del tamafio de la camada 
mostró menos sensibilidad al estrés y las enfermedades que la heterogénea (Argente et 
al., 2019). La diferente susceptibilidad a la enfermedad entre líneas estaría relacionada 
con una respuesta inmune diferente. Se sabe que los ácidos grasos desempefian diversas 
funciones en las células inmunitarias, modulando su respuesta; por ejemplo, se ha 
encontrado que regulan la fagocitosis, la producción de especies reactivas de oxígeno, la 
producción de citocinas y la migración de leucocitos, y también interfieren con la producción 
de antígenos por parte de los macrófagos (Yagoob y Calder, 2007). Por lo tanto, las líneas 
divergentes para la sensibilidad ambiental podrían mostrar diferencias en los niveles de 
ácidos grasos en plasma. 

El objetivo de este estudio fue evaluar la respuesta correlacionada a la selección 
para la variabilidad del tamafio de la camada en el perfil de ácidos grasos en el plasma, con 
el fin de identificar biomarcadores específicos para la sensibilidad ambiental. 
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MATERIALES Y MÉTODOS 


Animales y disefio experimental 


Los animales de este estudio provienen de la duodécima generación de un 
experimento de selección divergente para la sensibilidad ambiental. La selección se basó 
en la varianza fenotípica del tamafio de la camada de cada hembra, después de corregir 
este carácter por el afio-estación y el orden de parto-lactación (Blasco et al., 2017). Todos 
los animales fueron alojados en la granja de la Universidad Miguel Hernández de Elche 
(Espafia). Los conejos fueron alimentados con una dieta comercial estândar (Nutricun Elite 
Gra€, de Heus Nutrición Animal, La Corufia, Espafia). Se les proporcionó comida y agua 
ad libitum. Las hembras se alojaron en jaulas individuales (37,5 x 33 x 90 cm3) bajo un 
fotoperíodo constante de 16 h de luz continua: 8 h de oscuridad continua, y con ventilación 
controlada durante todo el experimento. Las hembras se aparearon por primera vez a las 
18 semanas de edad. 

Se tomó una muestra de sangre de la vena de la oreja a la segunda monta en 10 
hembras no lactantes de la línea homogénea y en 12 hembras no lactantes de la línea 
heterogénea en tamafio de la camada. El segundo apareamiento es un punto clave y 
representativo en la vida reproductiva de la hembra. Las muestras se recogieron en tubos 
con K3-EDTA, y se centrifugaron inmediatamente a 4000 rpm durante 15 min. El plasma 
se almacenó a -80 ºC hasta el análisis de lípidos. Todos los procedimientos experimentales 
con animales fueron aprobados por el Comité de Ética en Investigación de la Universidad 
Miguel Hernández de Elche (Número de referencia 2017/VSC/PEA/00212), de acuerdo con 
las Directivas del Consejo 98/58/CE y 2010/63/UE. 

Análisis químicos 

Se tomaron 200 ql de muestra de plasma en un tubo de vidrio con tapa scrum. La 
muestra de plasma se procesó siguiendo el método de Shirai et al. (2005). En el último 
paso, después de la evaporación con isooctano bajo una corriente de nitrógeno gaseoso, 
se afiadieron 200 gl de hexano. Los ácidos grasos se midieron mediante un cromatógrafo 
de gases (GC-17A, Shimadzu, Kyoto, Japón) con detector de ionización de llama (FID), 
equipado con un CP-Sil 88 para columna capilar FAME (100 m x 0,25 mm x 0,36 mm; 0,20 
um espesor de película; Agilent Technologies, Madrid, Espafia). El gas portador fue helio 
(flujo de 1,2 ml / min) con una inyección de 1:1. La rampa de temperaturas fue de 45 “ºC 
durante los primeros 4 min; de 13 “C/min hasta 175 “ºC a continuación; 4 ºC /min hasta 
215 “C mantenidos durante 30 min; temperatura del inyector, 250 ºC; y temperatura del 
detector, 260 ºC. Los ácidos grasos se identificaron comparando los tiempos de retención 
con el estândar FAME MIX (CRM47885, Supelco, Espafia). 
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Análisis estadístico 


El modelo utilizado para analizar el perfil de ácidos grasos plasmáticos incluyó los 
efectos de línea (línea homogénea y heterogénea) y mes (marzo, abril y mayo). Todos 
los análisis se realizaron utilizando metodologia bayesiana (Blasco, 2017). Consideramos 
que había suficiente evidencia de una diferencia entre la línea heterogénea (H) y la línea 
homogénea (L) si la probabilidad (P) de la diferencia era mayor de 0,90. 


RESULTADOS Y DISCUSIÓN 


En comparación con la línea heterogénea (Tabla 1), la línea homogénea mostró una 
mayor cantidad de ácidos grasos saturados (SFA) (+3,98 ng / ml P = 0,90 para C14:0; +2,30 
ng/ml P = 0,98 para C15:0 ; +54 ng / ml P = 0,90 para C16:0 y +29 ng /ml P = 0,90 para 
C18:0) y ácidos grasos monoinsaturados (MUFA) (+12,0 ng/ml P = 0,98 para C16:1 y +53 
ng/ml P= 0,90 para C18:1n9c). En relación a los ácidos grasos poliinsaturados (PUFA), la 
línea homogénea también exhibió una mayor cantidad de PUFAs n-3(+2.18 ng /ml P = 0,90 
para C18:3n3; +1.91 ng/ml P = 0,90 para C20:5n3). Con respecto a los PUFA n-6, la línea 
homogénea tuvo una mayor cantidad de C22:2n6 (+0,55 ng /ml P = 0,96) y menor cantidad 
de C20:3n6 (-3,66 ng / ml P = 0,96) y C20:4n6 ( -0,28 ng / ml P = 0,90) que la heterogénea. 











Línea Línea Homogénea (L) 
Heterogénea (H) 
Carácter (ng/ml) Media Media Da HPD9,,, P 
c8:0 7.14 4.09 3.65 -4.38, 11.9 0.81 
C10:0 0.87 0.16 0.71 -0.59, 2.05 0.85 
C11:0 0.07 0.35 -0.28 -0.79, 0.24 0.87 
C12:0 1.70 1.29 0.41 -1.06, 1.93 0.71 
C14:0 12.4 16.4 -3.98 -11.16, 2.44 0.90 
C15:0 2.34 4.64 -2.30 -4.48, -0.20 0.98 
C16:0 244 298 -54 -169, 45 0.90 
C17:0 5.25 6.46 -1.21 -3.65, 1.31 0.84 
C18:0 121 150 -29 -73 17 0.90 
C21:0 0.80 0.05 0.75 -0.51, 2.05 0.84 
C22:0 1.07 0.10 0.97 -0.70, 2.40 0.79 
C24:0 1.59 0.45 1.14 -0.70, 0.15 0.85 
SFAs 400 483 -83 -249, 68 0.90 
C14:1 1.76 1.55 0.21 -1.72, 2.10 0.59 
C15:1 2.35 2.47 -0.12 -1.78, 1.54 0.55 
C16:1 13.7 25.17 -12.0 -283.7, 0.1 0.98 
C17:1 1.79 2.49 -0.70 -2.10, 0.72 0.85 
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C18:1n9c 176 229 -58 -135, 32 0.90 


C20:1 1.16 0.94 0.22 -0.98, 1.41 0.65 
MUFAs 197 262 -65 -159, 31 0.92 
C18:2n6c 271 297 -26 -145, 95 0.67 
C18:3n6 0.67 3.48 -2.81 -5.69, 1.17 0.84 
C18:3n3 5.57 7.75 -2.18 -6.19, 0.98 0.90 
C20:2 2.31 2.28 0.02 -1.46, 1.55 0.51 
C20:3n6 4.44 0.76 3.66 -0.34, 7.80 0.96 
C20:3n3 0.87 0.90 -0.03 -1.48, 1.47 0.51 
C20:4n6 0.52 0.24 0.28 -0.12, 0.82 0.90 
C22:2n6 0.36 0.91 -0.55 -1.16, 0.07 0.96 
C20:5n3 1.59 3.50 -1.91 -7.59, 1.62 0.90 
C22:6n3 1.77 0.36 1.41 -1.81, 3.71 0.81 
n-3 PUFAs 9.87 11.97 -2.40 -9.59, 1.75 0.90 
PUFAs 286 317 -31 -150, 98 0.67 





D,,,: media de la diferencia entre líneas heterogéneas y homogéneas. HPD95%: región de densidad 
posterior más alta al 95%. P: probabilidad de que la diferencia sea> O cuando D,, > 0, y probabilidad 
de que la diferencia sea <O cuando D,, <0. SFA: ácidos grasos saturados. MUFA: ácidos grasos 
monoinsaturados. PUFA: Ácidos grasos poliinsaturados. 


Tabla 1. Perfil de ácidos grasos en plasma a la monta en las líneas heterogénea y homogénea para la 
variabilidad del tamafio de la camada. 


Los ácidos grasos participan en la regulación y activación de la respuesta inmune 
innata y adaptativa mediante la producción y síntesis de citocinas proinflamatorias, así 
como prostaglandinas, leucotrienos, tromboxanos y lipoxinas (Pompéia et al., 2000 en 
humanos). Por ejemplo, se ha encontrado un efecto negativo del ácido palmítico (C16:0) y 
el ácido esteárico (C18:0) sobre la proliferación de linfocitos y un efecto protector del ácido 
palmitoleato (C16:1) y el ácido oleico cis-9 (C18:1n9c) (Gorjão et al., 2007; Carrillo et al., 
2012; Chan et al., 2015 en humanos), mientras que ácido a linolénico (ALA, C18:3n3), 
ácido eicosapentaenoico (EPA, C20:5n3) y docosahexaenoico ácido (DHA, C22:6n3) han 
mostrado un efecto antiinflamatorio por medio de la disminución en la producción de IL- 
1, IL-2, IL-6, TNF así como de la prostaglandina E2 y el leucotrieno B4 (Katayama et al., 
2003 y Kelley et al., 1999 en humanos; Rodríguez et al., 2019 en conejos). La proteína 
C reactiva es una proteína de la fase aguda secretada por los hepatocitos durante la 
inflamación, en respuesta a las citocinas proinflamatorias (Pepys y Hirschfield, 2003), 
siendo un biomarcador útil en procesos inflamatorios. En un estudio anterior, Argente et al. 
(2019) detectaron una menor concentración de linfocitos y proteína C reactiva en la línea 
homogénea; por tanto, esta línea parece ser menos sensible a enfermedades y procesos 
inflamatorios que la heterogénea. Estos hallazgos concuerdan con una mayor cantidad de 
ácido palmítico, ácido esteárico, ácido a- linolénico y ácido eicosapentaenoico en la línea 
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homogénea. 


CONCLUSIÓN 


La selección por sensibilidad ambiental muestra una respuesta correlacionada en 
el perfil de ácidos grasos plasmáticos, corroborando una menor respuesta inflamatoria y 


menor sensibilidad a enfermedades en la línea homogénea. 
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RESUMEN: En Calvilo, Aguascalientes, se 
produce guayaba dese 1824. Con el objetivo de 
favorecer el crecimiento de la planta y mejorar la 
calidad de frutos, durante 2015-2017 se condujo 
el experimento. Para ello, se elaboró el abono 
de frutas, Té de estiércol y Té de composta, 
cada uno se aplicó al suelo, al follaje y al 
suelo y follaje más un testigo, generândose 10 
tratamientos, cada uno se repitió 4 veces, y se 
usó un arreglo completamente al azar, la unidad 
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experimental fue un árbol de dos afios de edad 
de la variedad “China”. Las aplicaciones (1 litro/ 
árbol) se realizaron cada 15 días. Los abonos 
en el suelo causaron diferencias significativas 
(P<0.05) en Da, MO y N. El P se incrementó con 
los tratamientos aplicados al suelo y foliares y El 
K, Ca y Mg, resultaron excesivos. En la solución 
del suelo el TÉCS aumentó en 3.5 mg/kg, el N y 
el Na en 16.45 mg/kg. En hojas el Na incrementó 
en 57.4 mg/kg y en fruto los tratamientos TéEF, 
TéCFS y AFF aumentaron su tamafio. Se 
concluye que los abonos orgánicos líquidos son 
una buena alternativa ecológica para mejorar la 
fertilidad del suelo y la nutrición del guayabo. 
PALABRAS CLAVE: Psidium guajava, Fertilidad, 
Nutrición, Abonos Orgánicos. 


ABSTRACT: In  Calvilo, Aguascalientes, 
guayaba from 1824 is produced. With the 
objective of favoring the growth of the plant and 
improving the quality of fruits, during 2015-2017 
the experiment is carried out. For this purpose, 
the fruit allowance was made, with manure and 
compost, each one was applied to the soil, to the 
leaves and to the soil and to the soil more than one 
testigo, generating 10 treatments, each one was 
repeated 4 times, and it was used a completely 
random arrangement, the experimental unit was 
a tree from two years of age of the “China” variety. 
The applications (1 liter/ftree) will be carried out 
every 15 days. The allowances in the soil caused 
significant differences (P<0.05) in Da, MO and N. 
El P if increased with the treatments applied to 
the soil and foliar and El K, Ca and Mg, resulted 
in excessive. In the suelo solution the TéCS 
increased to 3.5 mg/kg, the N and the Na to 16.45 
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mg/kg. Today, the TéEF, TéCFS and AFF treatments increased in size by 57.4 mg/kg. It is 
concluded that liquid organic allowances are a good ecological alternative to improve soil 
fertility and guayabo nutrition. 

KEYWORDS: Psidium guajava, Fertility, Nutrition, Organic Allowances. 


INTRODUCCIÓN 


El guayabo (Psidium guajava L.) pertenece a la familia Myrtaceae, su área ecológica 
se encuentra entre los paralelos 30º de latitud al norte y sur del ecuador (Alfonso, 2014). 
México, es considerado centro de diversidad de la especie por la amplia variabilidad 
existente (Singh, et al., 2015). Se adapta a diversos suelos y climas, puede comportarse 
como árbol caducifolio y perennifolio a la vez, según la disponibilidad de agua. 

El guayabo fue introducido al municipio de Calvillo, Aguascalientes en el afio 1824 
y en 1948 al Cafión de Juchipila, Zacatecas, actualmente se dedican al cultivo 2500 
productores y se cosechan 95 mil toneladas de fruta al afio, convirtiêndose en la zona 
productora de guayaba más importante de México con 62.6% de la superficie plantada en 
el âmbito nacional (Arif, et al., 2015). 

La degradación del suelo es consecuencia directa de la labranza equivocada, como 
el uso desproporcionado de fertilizantes químicos, eliminación de residuos orgánicos y en 
consecuencia empobrecimiento del suelo en materia orgánica (Gerardo, et al., 2005). La 
degradación química consiste en la disminución de la fertilidad del suelo y se produce 
cuando la entrada de nutrientes es menor a su salida por las cosechas, los incendios y las 
lixiviaciones (IA Ciampitti, et al., 2008). La degradación física del suelo es la destrucción de 
la estructura por la ruptura de los agregados del suelo, los poros se sellan con material fino, 
dificultando la infiltración del agua (Rucks et al., 1994). 

Existen testimonios que dan fe del efecto benéfico de los abonos orgánicos líquidos 
en el mejoramiento de la reserva nutritiva del suelo (Galindo, et al., 2007). Con la aplicación 
de abonos orgánicos líquidos en el cultivo del arroz (Oriza sativa L.) se han logrado 
excelentes resultados en la calidad y rendimiento, también en el cultivo de frijol (Phaseolus 
vulgaris L.), en altura de planta a los 30 días (Ochoa, et al., 2009). 

En condiciones controladas en el cultivo del tomate (Lycopersicum esculentum L.) 
aplicaron humus líquido al suelo y a las hojas, encontrando una respuesta positiva en cuanto 
al contenido de materia seca en las plantas, número de foliolos, volumen radicular y el peso 
seco de las raíces (Castellanos & Pratt, 1981). También (Cabrera, et al., 2009) sefialan que 
las enmiendas húmicas favorecen el enraizamiento, ya que desarrollan y mantienen un 
sistema radicular joven y vigoroso durante todo el ciclo de cultivo. El desarrollo radicular 
de la planta con aporte de enmiendas húmicas es enorme, y esto hace que el crecimiento 
sea mucho más rápido, debido a que absorbe más elementos nutritivos, y aumenta la 
producción. En este mismo sentido varios autores (Rosemary, et al., 2017) sefialan el efecto 
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estimulante de los ácidos húmicos en la formación de raíces al acelerar la diferenciación del 
punto de crecimiento. Por ello, en la investigación se fijó como objetivo mejorar la fertilidad 
del suelo, favorecer el crecimiento de la planta y mejorar la calidad de frutos del guayabo 
(Psidium guajava L.), en respuesta a la aplicación de abonos orgánicos líquidos. 


MATERIALES Y MÉTODOS 


Localización y características del lugar 


El experimento se realizó durante 2015, 2016 y 2017 en una huerta de guayabo 
localizada en San Tadeo del municipio de Calvillo, Aguascalientes, México, con coordenadas 
21º55” Latitud Norte y 102º41" Longitud Oeste, a una altitud de 1801 msnm. El clima es 
semi-cálido, con una precipitación media anual de 660 milímetros, de ocurrencia irregular, 
la temperatura media anual fluctúa entre 18 y 22 “C y una evapotranspiración potencial que 
alcanza valores de 2500 milímetros por afio. Tradicionalmente, en esta zona se ha utilizado 
la extracción de agua subterránea mediante bombeo con fines de riego, pero el recurso 
hídrico cada afio es escaso y de mediana calidad. 


Selección de árboles 


En la investigación se utilizaron 40 árboles de 2 afios y medio de edad, plantados 
a una distancia de 7m x 7m con riego por micro-aspersión con dos aspersores por árbol, 
el gasto por emisor fue de 4 litros por hora y una frecuencia de aplicación de 1 vez por 
semana. En la huerta, se tomaron 10 hileras contiguas de árboles, a cada una se le asignó 
un tratamiento, se etiquetaron cuatro árboles uniformes y cada árbol fue una unidad 
experimental, se realizó la poda y al inicio se midieron las variables indicadas (Tabla 1). 








Hileras Am DT cm DR cm NF NC NFr 
1 1.92 2.5 1.8 1 9 6 
2 1.91 2.8 2.0 1 17 1.5 
3 1.94 1.8 1.8 un 12 0 
4 2.00 2.3 1.7 2 26 4 
5 1.80 2.0 1.8 1.8 17 0 
6 2.00 3.0 3.0 0 8 8 
7 1.90 2.6 2.1 1 39 1 
8 1.80 2.5 2.7 1.5 16 36 
9 1.90 2.0 1.9 0 35 4 
10 2140 2.8 3.9 1.5 21 12 





Los valores indicados son el promedio de cuatro plantas. A=Altura (m), DT= Diámetro de tronco (cm), 
DR=Diámetro de ramas (cm), NF=Número de flores, NC=Número de cruceros y NFr=Número de frutos. 


Tabla 1. Caracterización de plantas al inicio de la investigación. 
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Disehio de tratamientos 


Fueron definidos 10 tratamientos, cada tratamiento se repitió 4 veces y se usó un 
disefio completamente al azar, la unidad experimental fue un árbol (Tabla 2). 








Tratamiento Clave Descripción 

1 Té ES Té de estiércol aplicado al suelo 

2 Té CS Té de composta aplicado al suelo 

3 AFS Abono de frutas aplicado al suelo 

4 Té EF Té de estiércol aplicado al follaje 

5 Té CF Té de composta aplicado al follaje 

6 AFF Abono de frutas aplicado al follaje 

7 Té ESF Té de estiércol aplicado al suelo y follaje 
8 Té CSF Té de composta aplicada al suelo y follaje 
9 AFSF Abono de frutas aplicado al suelo y follaje 
10 T Testigo (sin aplicación) 





Tabla 2. Tratamientos, clave y descripción. 


Muestreo de suelo 


Fueron tomadas cuatro sub-muestras de suelo en el área que cubre la copa de los 
árboles a la mitad de distancia del tronco a la orilla de copa, con orientación N-S y E-O 
a una profundidad de O a 20 cm, se utilizó una barrena tipo Oakfield para muestreo de 
suelo (38 mm diámetro por 216 mm largo), se depositó el suelo extraído sobre una manta, 
se mezcló manualmente y mediante cuarteos se dividió eliminando dos cuartos opuestos 
hasta dejar aproximadamente un kilo de muestra compuesta. En el laboratorio las muestras 
se extendieron sobre papel absorbente para su secado a temperatura ambiente, con buena 
ventilación y a la sombra. Una vez secas, se pasaron por un tamiz marca Retsch de malla 
2 mm, para obtener un tamafio homogéneo de partícula. A cada muestra se le determinó 
las variables: pH, con medidor de mesa modelo SM-3BW rango de pH 0.00-14.00 exactitud 
+0.01, en una solución suelo: agua 1:2 (p/v); CE, con conductímetro, modelo Amprobe WT- 
20 en término de sales solubles dS/m; CIC, por el método de acetato de amonio 1N a pH 
7, en referencia a la capacidad que tiene un suelo para retener y liberar iones positivos. 
La MO por el método de Walkley-Black, N con Destilación Semimicro-kjeldahl, P método 
Olsen, K se cuantificó por espectrofotometría de emisión de flama, Ca y Mg se cuantificaron 
por Acetato de amonio 1.0 N pH 7 y espectrofotometría de absorción atómica; Fe, Cu, Zn y 
Mn por DTPA y Absorción Atómica y B por Fotocolorimetría de Azometina-H. 


Muestreo de solución de suelo 


Previo a la instalación de la sonda en campo, se le provocó vacio hasta 70-80 
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centibares con bomba de succión. Se sumergió la capsula en un recipiente con agua y se 
dejó durante 12 horas. La sonda consta de un cilindro de PVC unido a una cápsula porosa 
de porcelana de forma cilíndrica con terminal semiesférica, a través de la cual penetra la 
solución del suelo al aplicar una diferencia de presión mediante vacío. El cilindro de PVC 
está a su vez sellado por un tapón de goma, al que atraviesa un tubo de diámetro pequefio 
y paredes semirrígidas que se conectan al sistema de vacío. La muestra se recogió en 
botellas de plástico de 240 ml con tapa de rosca. 

La toma de muestras consistió en cargar la sonda, creando una depresión de 70 cbar 
24 horas antes del siguiente riego y recogiendo la muestra un momento antes del mismo. 
El mecanismo como la solución de suelo pasa a la cápsula, es parecido al mecanismo 
de absorción pasiva de las plantas, esto hace que se pueda considerar como método de 


referencia cuando interviene un cultivo. 


Muestreo foliar 


El muestreo foliar se realizó 6 meses después de la poda y abonado orgánico para 
garantizar un estado nutrimental adecuado de las plantas, así como apariencia vigorosa 
y aspecto sano; ya que no se debe utilizar plantas que hayan estado sometidas a alguna 
clase de estrés; para ello, se tomaron de cada planta 20 hojas recientemente maduras 
de la periferia y parte media de la copa de los árboles. Las muestras se colocaron en 
bolsas plásticas para su transporte al laboratorio en una hielera, luego, se lavaron con 
agua destilada y secaron al ambiente, con la finalidad de eliminar cualquier residuo que 
pudiera interferir con el análisis químico. Posteriormente, se pesaron en una balanza 
Mettler Pc 4400 (el peso se expresó en gramos + 0.01) para obtener la masa fresca (MF), y 
se procedió a colocarlas en una estufa a 65ºC por un tiempo de 24 a 48 horas hasta obtener 
un peso constante y luego se molieron. Enseguida, se tomó 1 gramo de la muestra obtenida 
de cada planta, conformando una muestra compuesta para cada grupo, la preparación de 
este material (muestra molida) y los métodos analíticos empleados se realizaron según 
los procedimientos sefialados por la A.O.A.C. Se determinó nitrógeno (N) por el método 
Kjeldahl, fósforo (P) por el método de Meta-vanadato-molibdato de Amonio, potasio (K), 
calcio (Ca), magnesio (Mg) y sodio (Na) por el método de Espectro Fotometría de Llama. 
Los resultados se expresaron como la concentración de minerales por masa seca (MS) 
(gramos de mineral/100 gramos de materia seca). 

Para el análisis estadístico se utilizó el paquete Minitab 17, se realizó análisis 
de varianza y pruebas de comparación de medias para las variables con resultados 
significativos (P<0.05), las variables físico-químicas del suelo mediante la prueba LSD os 
y para las variables químicas de la solución del suelo y foliares, se utilizó la prueba de 
Tukey, 0.05. 
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RESULTADOS Y DISCUSIÓN 


El análisis químico de los abonos líquidos mostró que el té de estiércol resultó con 
mayor concentración de N, P, Ca, Na, K y Cu (mg/l); le siguió el abono de frutas con 
contenidos elevados de N, P, Ca, Mg, K, Fe, Mn y Zn (mg/1) y finalmente el Té de composta 
con excepción del contenido de Na (mg/l) fue el más pobre (Tabla 3). 








Abono N P Ca Mg Na K Fe Cu Mn Zn 
AF 9.7 10.9 201.3 29.5 71.5 13001 60.1 0.1 4.3 1.3 
TéC 2.3 0.2 49.1 13.2 125.8 591.4 1.1 0.1 0.2 0.4 


TéE 56.2 12.6 270.5 28.0 536.7 3285.7 5.2 0.7 0.4 0.7 





AF=Abono de frutas, TéC=Té de composta, TéE=Té de estiércol. N=Nitrógeno, Micro-Kjenhdal (mg/1); 
P=Fósforo, Molibdato de amonio (mg/l); Ca=Calcio, Mg=Magnesio, Na=Sodio, K=Potasio, Fe=Fierro, 
Cu=Cobre, Mn=Manganeso y Zn=Zinc por Absorción atómica (mg/l). 


Tabla 3. Análisis químico de abonos orgánicos líquidos. 


Variables físico-químicas del suelo 


El suelo pertenece a la clase textural arcillosa (Tabla 4); la textura es una propiedad 
que afecta directamente la densidad aparente, ya que a medida que los suelos se compactan 
disminuye la porosidad y aumenta la densidad, como lo sefiala (Cogger, et al., 2016). 








Tratamientos Arena Limo Arcilla Clase Da MO pH CE 
% % % textural (g/em?) (1.6-3.5%) — (6.6-7.3 moles) dS/m 
TéEF 14.4 240 61.6  Arcilosa 1.080b 2.66 bcde 74a 2.249 ab 
TéES 18.4 220 59.6  Arcilosa 1.112b 3.53 abc 74a 2.274a 
TéEFS 12.4 16.0 71.6  Arcilosa 1.129ab 3.33 abed 75a 2.201 ab 
TéCF 20.4 16.0 63.6  Arcilosa 1.112b 407 a 74a 1.491 b 
Técs 14.4 18.0 67.6  Arcilosa 1.119ab 3.66 ab 74a 1.921 ab 
TéCSF 12.4 16.0 71.6 Arcilosa 1.119ab | 2.52cde 75a 2.192 ab 
AFF 6.4 22.0 71.6  Arcilosa 1.119ab 1.98e 76a 2.268 a 
AFS 20.4 20.0 59.6  aArcilosa 1.168a 2.49 cde 75a 1.506 b 
AFFS 18.4 20.00 61.6  Arcilosa 1.169a 2.32 de 75a 1.884 ab 
Testigo 14.4 220 63.6  Arcilosa 1.168a 2.42 de 76a 1.967 ab 
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Significancia id Ei Ns ns 





Tratamientos: Té de estiércol foliar (TéEF); Té de estiércol suelo (TéES); Té de estiércol foliar y suelo 
(TéEFS); Té de composta foliar (TéCF); Té de composta suelo (TéCS); Té de composta suelo y foliar 
(TéCSF); abono de frutas foliar (AFF); abono de frutas suelo (AFS); abono de frutas foliar y suelo 
(AFFS); testigo, sin abono. 


Variables físico-químicas: Arena %, Limo % y Arcilla %, Bouyoucos, Da=Densidad aparente (g/cm?) 
probeta, MO=Materia Orgánica % Walkey y Black, pH relación 1:2 potenciómetro, CE=Conductividad 
Eléctrica dS/m conductímetro. Sinificancia entre tratamientos: ns= No significativo, *= Significativo (F>p 
0.05) y **=Altamente significativo (F>p 0.01). 


Valores con la misma letra dentro de las columnas son estadísticamente iguales, de acuerdo al método 
LSD de Fisher y una confianza de 95%. 


a=0.05 


Tabla 4. Análisis físico-químico del suelo. 


Los tratamientos causaron diferencias significativas (P<0.05) en densidad aparente 


(gr/em?) y la prueba de comparación de medias LSD 


1005 Sefialó tres grupos; dentro de cada 


grupo los tratamientos son estadísticamente iguales y entre grupos son diferentes. El grupo 
(a) se integró por los tratamientos AFS, AFFS y Testigo, el grupo (ab) por los tratamientos 
TéEFS, TéCS, TéCSF, AFF y el grupo (b) por los tratamientos TéEF, TéES, TéCF. Lo 
que evidenció alta densidad aparente con el abono de frutas aplicado al suelo AFS, así 
como aplicado al suelo y al follaje AFFS y testigo, mientras que el resto de los abonos 
la redujeron (Tabla 4; Figura 1). Los valores bajos de densidad aparente son propios de 
suelos porosos, con buena circulación de aire y agua y fácil penetración de raíces. Como lo 
menciona (Lordan et al., 2015) la densidad aparente es una propiedad dinámica que varía 
con las condiciones estructurales del suelo, la agregación en los suelos tiende a aumentar 
el espacio poroso por tanto disminuye la densidad aparente. Los valores altos de densidad 
aparente son propios de suelos compactos y poco porosos, con aireación deficiente e 
infiltración lenta del agua. 
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Da g/cm3 









F=2.80 
p=0.016 

S=0.0353862 
R2= 45.63% 








AFF AFFS AFS TéCF TéCs  TéCSF  TéEF  TéEFS  TéES Testigo 


Tratamientos 


Figura 1. Variación de la densidad aparente causada por los tratamientos. 


Los abonos causaron diferencias significativas (P<0.01) en el porcentaje de MO; 


se realizó la comparación de medias con el método LSD 


1005: fesultó la formación de 


ocho grupos (Tabla 4, Figura 2), el primero (a) por el tratamiento TéCF; el segundo (ab) 
por el tratamiento TéCS, el tercero (abc) por el tratamiento TéES, el cuarto (abcd) por 
el tratamiento TéEFS, el quinto (abcde) por el tratamiento TéEF, el sexto (cde) por los 
tratamientos TéCSF y AFS, el séptimo (de) por los tratamientos AFFS y Testigo, finalmente 
el octavo (e) por el tratamiento AFF. Este resultado evidenció que con excepción de los 
tratamientos abono de frutas aplicado al follaje AFF y abono de frutas aplicado al follaje 
y al suelo AFFS, en comparación al testigo los tratamientos incrementaron el porcentaje 
de MO en el suelo, destacando el té de composta aplicado al follaje TéCF que registró un 
incremento en 1.65% (Tabla 4; Figura 2). Así mismo lo confirma (Horrocks, et al., 2015) y 
concuerda con la (Prado & Natale, 2004), que establece un contenido medio de MO para 
valores que se encuentren entre 1.6- 3.5%. 
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F=3.39 
p=0.006 
S=0.751521 
R2=50.42% 


MO % 





AFF AFFS AFS Técr Técs Técsr TéEF TéEFS TéES  Testigo 


Tratamientos. 


Figura 2. Variación de la materia orgánica causada por los tratamientos. 


En pH (moles/l) no se encontraron diferencias (P>0.05) causadas por los tratamientos 


(Tabla 4); sin embargo la prueba de comparación de medias LSD evidenció que en los 


q=0.05 
abonos té de estiércol aplicado al follaje TéEF, té de estiércol aplicado al suelo TéES, té 
de composta aplicado al follaje TéCF y té de composta aplicado al suelo TéCS el pH, se 
encuentra en el rango de neutralidad 6.5 — 7.4 (moles/l) y en el resto de los tratamientos 
ligeramente alcalino 7.5 — 7.9 (moles/l) (Horrocks, et al., 2015). 

En CE dS/m no se encontraron diferencias (P>0.05) de tratamientos, la prueba de 
comparación de medias LSD «os identificó tres grupos. El grupo (a) se integró por los 
tratamientos TéES y AFF, el (ab) por los tratamientos TéEF, TéEFS, TéCS, TéCFS, AFFS 
y Testigo y el (b) por los tratamientos TéCF y AFS y se demostró que el té de estiércol 
aplicado al follaje TéEF, incrementó la CE en 0.282 dS /m , té de estiércol aplicado al 
suelo TéES en 0.307 dS/m, Té de estiércol aplicado al follaje y al suelo TÉEFS en 0.234 
dS/m, té de composta aplicada al suelo y al follaje TéCSF en 0.235 dS/m y abono de frutas 
aplicado al follaje AFF en 0.301 dS/m (Tabla 4). Estos resultados indican que los abonos 
no causaron incrementos en las sales que puedan originar algún problema en la planta. El 
análisis de la CE en suelos se hace para establecer si las sales solubles se encuentran en 
cantidad suficiente como para afectar la germinación normal de las semillas o la absorción 
de agua por parte de las plantas. La acumulación de sales solubles en el suelo se atribuye 
principalmente a problemas de drenaje y a la aplicación de riegos continuados, seguidos de 
evaporación y sequia como lo sefiala (Horrocks, et al., 2015). 
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Contenido de nutrientes en el suelo 


Los abonos causaron diferencias significativas (P<0.05) en el contenido de N, la 


comparación de medias LSD, , os formó cinco grupos (Tabla 5, Figura 3), el (a) se integró por 


=0.05 
el tratamiento abono de frutas aplicado al follaje y al suelo AFFS, el (ab) por el tratamiento 
abono de frutas aplicado al suelo AFS, el (abc) por los tratamientos abono de frutas foliar 
AFF y té de estiércol foliar TéEF;, el (bc) por el tratamiento te de estiércol al suelo TéES y el 
(c) por los tratamientos té de estiércol foliar y suelo TéEFS, té de composta foliar TéCF, té 
de composta suelo TéCS, té de composta suelo y foliar CSF y Testigo; en comparación al 
testigo se evidenció que el abono de frutas aplicado al follaje y al suelo AFFS, lo incrementó 
en 18.27 mg/kg, abono de frutas aplicado al suelo AFS, en 16.32 mg/kg, abono de frutas 
aplicado al follaje AFF en 5.0 mg/kg, el Té de estiércol foliar TÉEF en 8.46 mg/kg, el Té de 
estiércol aplicado al suelo TéES en 5.2 mg/kg y el Té de composta foliar (TéCF) en 1.03 
mg/kg. Estos resultados sugieren que es promisoria la producción orgánica (Lordan et al., 











2015). 
TRAT N (10-20 P(15-25  K(175-280 Ca (1000- Mg (60-180 Na 
mg/kg) mg/kg) mg/kg) 2000 mg/kg) mg/kg) mg/kg 
TéEF 17.09abc  12.54ab 1426 ab 3095 ab 562.8 ab 121.6a 
TéES 12.21 bc 5.78b 1305 ab 2877 ab 504.4 b 107.7 a 
TéEFS 6.24c 14.52 ab 1392 ab 3636 a 543.8 ab 152.2a 
TéCF 9.66 c 13.86 ab 1089 b 2710b 517.1 ab 107.3a 
Técs 841 c 7.42 ab 1300 ab 2835 ab 551.1 ab 158.9a 
TéCSF 8.68 c 15.34 ab 1267 b 2895 ab 536.1 ab 141.5a 
AFF 13.83 abc 11.22 ab 1861 a 3230 ab 611.1a 177.9a 
AFS 24.95 ab 11.22 ab 1332 ab 2812 ab 463.7 b 1184 a 
AFFS 26.9 a 16.66a 1367 ab 2867 ab 539.6 ab 109.0 a 
Testigo 8.63c 11.39 ab 1488 ab 3470 ab 536.5 ab 107.9 a 
Significancia x ns ns ns Ns ns 
.... Cont. Cuadro 
TRAT CiC(Cmol/ Fe(2.5- Mn(>2mg/ Zn (>1.5 mg/ Cu (>2.0 B (0.5-2.0 
kg) 5.0 mg/ kg) kg) mg/kg) mg/kg) 
kg) 
TéEF 24.27 ab 9.26 ab 8.987 abc 0.3005 b 0.1925 ab 0.2800 b 
TéES 22.33 ab 11.09 ab 10.50 ab 0.3985 ab 0.1895 ab | 0.4000 ab 
TéEFS 26.87 a 12.31 ab 9.319 abc 0.604 a 0.2140 a 0.300 b 
TéCF 21.05 b 14.67 a 12.23a 0.4080 ab 0.1710 a 0.640 a 
Técs 22.71 ab 7.28b 6.29 be 0.2535 b 0.1405 bc 0.2400 b 
TéCSF 22.15 ab 6.676b 4.924 c 0.2340 b 0.13000 c 0.2400 b 
AFF 26.51 a 6.60 ab 4.88 c 0.2200 b 0.1450 bc 0.320 b 
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AFS 22.38 ab 747b 6.05 bc 0.2500 b 0.1775abc 0.1200b 


AFFS 22.54 ab 10.51 ab 8.51 abc 0.3225 b 0.1920 ab 0.1600 b 
Testigo 25.82 ab 8.73 ab 8.48 abc 0.2980 b 0.1735 abc 0.220 b 
Significancia ns ns * ns ns ns 


Tratamientos: Té de estiércol foliar (TéEF); Té de estiércol suelo (TéES); Té de estiércol foliar y suelo 
(TéEFS); Té de composta foliar (TéCF); Té de composta suelo (TéCS); Té de composta suelo y foliar 
(TéCSF); Abono de frutas foliar (AFF); Abono de frutas suelo (AFS); Abono de frutas foliar y suelo 
(AFFS); Testigo, sin abono. 


Significancia entre tratamientos: ns=No significativo, *= Significativo (F< p 0.05). Valores con la misma 
letra dentro de las columnas son estadísticamente iguales, de acuerdo al método LSD de Fisher y 
una confianza de 95 % 


a=0.05 


N= Nitrógeno total Micro-Kjenhdal (mg/kg), P=Fosforo disponible Olsen (mg/kg), K, Ca, Mg y Na 
intercambiable (mg/kg). CIC=Capacidad de intercambio catiónico (Cmol/kg), 


Tabla 5. Valoración de nutrientes en el suelo con la influencia de los abonos líquidos. 





F=2.44 


p=0.032 
S=9.20188 
R2=42.24% 


N. Inorgánico mg/kg 





AFF AFFS AFS TécF TéCs  TéCSF  TéEF  TéEFS  TéES Testigo 


Tratamientos. 


Figura 3. Variación de nitrógeno inorgánico en el suelo causado por los tratamientos. 


El análisis de varianza no detectó diferencias (P>0.05) en contenido de P, K, Ca, Mg 
y Na en el suelo, sin embargo para estos elementos se realizó la comparación de medias 
LSD, oos Y para P se formaron los grupos (a) con el tratamiento AFFS, el (ab) con los 
tratamientos té de composta suelo y foliar TéCSF, té de estiércol suelo y foliar TéEFS, té de 
composta foliar TéCF;, té de estiércol foliar TÉEF, abono de frutas foliar AFF, abono de frutas 
suelo AFS, té de composta suelo TéCS y Testigo; y el (b) por el tratamiento té de estiércol 
suelo TéES. Las diferencias encontradas en comparación con el testigo fueron: bono de 
frutas foliar y suelo AFFS en 5.27 mg/kg, el té de composta foliar y suelo TéCSF en 3.95 
mg/kg, el té de estiércol foliar y suelo TÉEFS en 3.13 mg/kg, el té de composta foliar TÉCF 
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en 2.47 mg/kg y el té de estiércol foliar TéEF en 1.15 mg/kg (Tabla 5). Estos resultados se 
respaldan con los reportados por (Prado & Natale, 2004). 

Para K se formaron los grupos (a) con el tratamiento abono de frutas foliar AFF, el 
(ab) con los tratamientos Testigo, té de estiércol foliar TéEF, té de estiércol foliar y suelo 
TéEFS, abono de frutas foliar y suelo AFFS, abono de frutas foliar AFF, té de estiércol suelo 
TéES, abono de frutas suelo AFS y té de composta suelo TéCS; y el (b) por el tratamiento 
té de composta foliar TéCF; y se evidenció con respecto al testigo, que el abono de frutas 
foliar AFF, lo incrementó en magnitud de 373 mg /kg (Tabla 5). 

Para Ca, se formaron los grupos (a) por el tratamiento TéEFS, el (ab) por los 
tratamientos AFF, TéEF, Testigo, TéCSF, TéCS, TéES, AFS, y el (b) por el tratamiento 
TéCF; y se evidenció que el tratamiento TéEFS lo incrementó en 163 mg/kg (Tabla 5). 

Para Mg, se formaron los grupos (a) por el tratamiento AFF, el (ab) por los tratamientos 
TéEF, TéCS, TéEFS, AFFS, TéCSF, Testigo, TéCF y el (b) por los tratamientos AFS y TéES; 
se evidenció que el AFF lo incrementó en 74.6 mg/kg (Tabla 5). 

Para Na formó el grupo (a) al que pertenecen todos los tratamientos y se evidenció 
que el tratamiento AFF lo incrementó en 70 mg/kg, el TéCS en 51 mg/kg y el TéCSF en 
33.6 mg /kg (Tabla 5). 

Los tratamientos aumentaron el contenido de Mn (P<0.05) en el suelo, la prueba de 


comparación de medias LSD evidenció que el tratamiento TéCF lo incrementó en 3.75 


a=0.05 
mg/kg, también aumentó el Zn en 0.3 mg/kg y en base a un valor de suficiencia >1.5 mg/ 
kg, se aprecia insuficiencia en el suelo, aunque el Cu aumentó ligeramente, su contenido 
es pobre ya que no fue >2, mientras que el B con rango de suficiencia de 0.5-2, el resultado 


sefiala que se encuentra en el límite inferior de suficiencia. 


Nutrientes en la solución del suelo 


La concentración de nutrientes en la solución del suelo es generalmente baja (Tabla 
6), como lo refiere (lbekwe et al., 2017) quienes evaluaron una concentración uM promedio 
de nutrientes en la solución en la capa arable (0-20 cm) de un suelo agrícola en: pH 7.7, K 
510, Ca 1650, Mg 490, N-NH, 48, N-NO, 3100, S-SO, 590, P-PO, 1.59, Zn 0.48 y Mn 0.002. 








TRAT N (mg/) Na (mg/l) Cu(mg1) Zn (mg/l) pH (molesf) 
TéEF 14.875 be 32.95 b 0.010 a 0.010 a 6.875 b 
TéES 14.000 be 38.70 ab  0.002a 0.000 a 7.350 ab 
TéEFS 12.070 c 48.85 a 0.005 a 0.005 a 7.675 a 
TéCF 16.270 abc 36.25 ab  0.005a 0.002 a 7.525 a 
Técs 20.300 a 43.42ab  0.010a 0.000 a 7.400 ab 
TéCSF 16.275 abc 42.20 ab  0.005a 0.000 a 7.675 a 
AFF 15.400 abc 35.42 b 0.005 a 0.000 a 7.500 a 
AFS 15.925 abc 33.98 b 0.002 a 0.000 a 7.850 a 
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AFFS 18.030 ab 35.00 b 0.000 a 0.000 a 7.675 a 


Testigo 16.800 abc 32.40 b 0.000 a 0.007 a 7.650 a 
Significancia e RA e E a 





Tratamientos: Té de estiércol foliar (TéEF); Té de estiércol suelo (TéES); Té de estiércol foliar y suelo 
(TéEFS); Té de composta foliar (TéCF); Té de composta suelo (TéCS); Té de composta suelo y foliar 
(TéCSF); Abono de frutas foliar (AFF); Abono de frutas suelo (AFS); Abono de frutas foliar y suelo 
(AFFS); Testigo, sin abono. 


N= Nitrógeno (mg/l), Na=Sodio (mg/l), Cu=Cobre (mg/]), Zn= Zinc (mg/l), y pH (moles/]). Significancia,*= 
Significativo; **=Altamente significativo (F< p 0.05). Valores con la misma letra dentro de las columnas 


son estadísticamente iguales, de acuerdo al método LSD os de Fisher y una confianza de 95 %. 


Tabla 6. Contenido de nutrientes en la solución del suelo. 


Los tratamientos contribuyeron en el contenido de N (P<0.01), la prueba de 
comparación de medias LSD os evidenció que el tratamiento TÉCS, con respecto al testigo 
incrementó la concentración en 3.5 mg/|, y Na también aumentó (p<0.01), sefialando que el 
tratamiento TéEFS lo incrementó en 16.45 mg/l (Tabla 6, Figuras 4 y 5). 


p=0.001 
S=2.12513 
R2=56.74 


ab 


Nitrógeno mg/1 





AFF AFFS AFS TéCF  TéCFS  TéCs TéEF TéEFS TéES  Testigo 


Tratamientos. 


Figura 4. Contenido de Nitrógeno en la solución del suelo. 
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55 
p=0.001 


S=5.23177 
so)  R2=55.64% 


ab 


307 








AFF AFFS AFS TéCF  TéCFS  Técs TéEF TéEFS  TéES  Testigo 


Tratamientos. 


Figura 5. Contenido de Sodio en la solución del suelo. 


Los nutrientes son transportados por el agua hacia la cercanía de las raíces por 
flujo de masa, la fuerza motriz es la transpiración de la planta. La cantidad de nutriente 
transportada es el resultado del volumen de solución (V,) multiplicado por la concentración 
del nutriente en la solución (C)) FM = V,* C,. Para que este mecanismo sea eficiente, la 
cantidad de nutriente requerida debe ser pequefia o la concentración del nutriente en la 
solución debe ser muy alta. El Mg, el N-NO y el S-SO,2 se absorben mayoritariamente por 
este mecanismo. Una vez que Ilegan a la superficie puede ocurrir que la planta absorba el 
agua pero no los nutrientes, los cuales se acumulan y pueden precipitar. 


Efecto de los tratamientos en el contenido nutrimental en hojas de guayabo 


Las plantas responden a la aplicación foliar de abonos orgánicos líquidos, al ser 
aplicados al follaje de los cultivos, permite aumentar la cantidad de raíces e incrementar 
la capacidad de fotosíntesis de las plantas, mejorando sustancialmente la producción y 
calidad de la cosecha (Prado & Natale, 2004). 

(Chetri, et al., (1999) encontraron que el contenido foliar de N y K fue bajo en el 
transcurso de la estación de crecimiento y se incrementaron en Agosto llegando al máximo 
de Septiembre a Diciembre, sin efectos significativos en contenido foliar de P, con algunas 
excepciones, el contenido de Mg exhibió un comportamiento similar al Ca. Las hojas de 
Allavad Safeda mostraron contenidos marcadamente altos de K, Ca y Mg en cv Luknow-49. 

En guayabo (Psidium guajava L.), como en otras especies, la composición 
nutrimental de las hojas varía con la edad y la posición de las mismas en los brotes, el 
tipo de brote (vegetativo o reproductor) y el estado fenológico del árbol. Por ejemplo, 
los valores de N, P K, Zn y Cu en hojas de guayabo disminuyen a medida que la hoja 
aumenta de edad, y Ca, Mg y Mn se incrementan, pero se estabilizan a los cinco meses de 
edad (Hundal, et al., 2007). Con la finalidad de realizar un diagnóstico preliminar y derivar 
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niveles críticos y rangos de suficiencia de nutrientes en la hojas de guayabo, (de Oliveira, 
et al., 213) encontraron que los nutrientes limitantes por deficiencia en orden decreciente 
fueron: N>Cu>P=K>Mn>Fe=Zn>S>B=Mg>Ca, y los elementos limitantes por exceso en 
orden decreciente fueron: B>Ca>Fe>Mn>S>Mg>Cu>P>Zn>N=K. Los rangos de índices 
apropiados DRIS fueron: N(24-48), P(2.4-3.1), K(21-29), Ca(6-8), Mg(1.9-2,9) y S(1.9- 
2,3) g/kg respectivamente (Antunes, et al.,2016) . Otros investigadores reportaron que 
los rangos de suficiencia de los nutrientes derivados de las normas DRIS fueron: N (1.41- 
1.65%), P (0.10-0.17%), K (0.51-0.97%), Ca (1.16-2.12%), Mg (0.31-0.51%) y S (0.18- 
0.28%) (Hundal et al., 2007). 








Tratamientos P Ca Mg K Fe Cu Mn 

% % % % (mg/kg) (mg/kg) (mg/kg) 

TéEF 0.17ab 1.60a 021a 1.09a 83.30a 4.55a 48.78 a 
TéES 0.17ab 1.10ab 0.16a 1.05a 80.88a 4.27a 36.15 a 
TéEFS 0.15b 0.96ab Ol7a 1.18a 78.35a 5.07a 36.80 a 
TéCF 0.14b 1.10ab 0.19a 0.9 ab 106.5a 5.25a 59.4 a 
Técs 0.13b 0.89ab 017a 1.02ab 82.1a 3.95a 42.13a 
TéCSF 0.19ab 0.90ab O17a 1.05a 72.10a 4.95a 48.45 a 
AFF 0.15b 0.76b 0.18a 0.98ab 73.13 4.25a 39.83 a 
AFS 0.15b 074b 018a 0.69b 93.3a 4.62a 40.50 a 
AFFS 0.19ab 0.86ab O17a 0.68b 129.4a 4.25a 54.85 a 
Testigo 0.23a 1.16ab 0.20a  0.68b 78.03a 4.82a 59.63 a 





Tratamientos: TéEF=Té de Estiércol Foliar; TéES=Té de Estiércol Suelo; TéEFS=Té de estiércol Foliar 
y Suelo; TéCF=Té de Composta Foliar; TéCS=Té de composta Suelo; TéÉCSF=Té de Composta Suelo 
y Foliar; AFF=Abono de Frutas Foliar; AFS=Abono de Frutas Suelo; AFFS=Abono de frutas Foliar y 
Suelo; Testigo=Sin abono. P (Vanado-molibdato); Ca, Mg, K, Fe, Cu y Mn (Absorción Atómica). 


Rangos de suficiencia: P (0.168-0.236 %); Ca (0.60 - 1.27 %); Mg (0.35 - 0.50 %); K (1.20 to 1.67 %); 
Fe (114 - 178 mg/kg); Cu (6-12 mg/kg). Mn (34 - 77 mg/kg). 


Tabla 7. Efecto de los abonos líquidos en el contenido nutrimental en hojas de guayabo. 


Los abonos no incrementaron el contenido de P en las hojas, los resultados sefialan 
porcentajes inferiores al testigo, sin embargo se detectaron diferencias significativas 
(F, 
tres grupos: (a), (ab) y (b). De acuerdo al diagnóstico del balance de nutrientes DRIS, en 


=3.05>p=0.010) entre los tratamientos y la comparación de medias (Tukey, formó 


0.05 005) 


los grupos (a) y (ab) el contenido foliar de P se encontró en el rango de suficiencia (0.168 
- 0.236%) con los tratamientos: Testigo, AFFS, TéCSF, TéES y TéEF e insuficiente en el 
grupo (b) con los tratamientos TéEFS, TéCF, TéCS, AFF y AFS (Tabla 7). 

En Ca foliar su contenido se encontró dentro del rango de suficiencia (0.60-1.27%) 
reportado por (Anjaneyulu & Raghupathi, 2009), destacando el tratamiento TéEF bajo el 
cual se registró un contenido excesivo, al superar en 44% al testigo. 
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En Mg, no se evidenció efecto de tratamientos, en todos los casos fue insuficiente 
para satisfacer los requerimientos dela planta. 
Los abonos no incrementaron el contenido de K en las hojas, sin embargo se 


detectaron diferencias significativas (F,,.=6.80>p=0.000) entre los tratamientos y la 


0.05 


comparación de medias (Tukey, formó tres grupos: (a), (ab) y (b). De acuerdo al 


o, 
diagnóstico del balance de nutrientes DRIS, sefialado por Anjaneyulu & Raghupathi, 2009, 
el contenido foliar fue inferior al óptimo (1.20 to 1.67 %). 

En la tabla 7, se aprecia que en el contenido foliar de P Mg, Fe, Cu y Mn no 
contribuyeron los abonos. En P, la prueba de comparación de medias de Tukey, 0.05, 
sefiala diferencias entre tratamientos, pero el valor más alto corresponde al testigo (sin 
abono), este resultado informa que en las hojas los valores de concentración se encuentran 
próximos al límite inferior del rango de suficiencia (0.14-0.18 %). En Ca el tratamiento té de 
estiércol aplicado al follaje (TéEF) aumentó en 0.44 %. En K el tratamiento té de estiércol 
aplicado al follaje (TéEF), incrementó en 0.41 %, el tratamiento té de estiércol aplicado al 
suelo (TéES) en 0.39 %, el té de estiércol aplicado al suelo y al follaje (TéEFS ) en 0.50 
%, el té de composta aplicado al follaje (TéCF) en 0.26 %, el té de composta aplicado al 
suelo (TéCS) en 0.34 %, el té de composta aplicado al follaje y al suelo (TéCSF) en 0.37 
% y el abono de frutas aplicado al follaje (AFF) en 0.30 %. En Na el tratamiento abono de 
frutas aplicado al follaje (AFF), incrementó en 57.4 mg/kg. Estos resultados en apariencia 
contradictorios se deben a la época de muestreo, la que se realizó en la etapa de floración 
ya que en la época de fructificación del guayabo, los valores de K son bajos, lo que se 
atribuye a su movilización hacia los frutos, también se han observado tendencias similares 
de N, Py Ken plantas de esta especie, porque sus valores son más altos en los brotes 
fructíferos que en los vegetativos. El contenido foliar de K, Ca y Mg también varía entre 
fechas, de manera que al disminuir K, se incrementa Ca y, en menor grado, Mg (Da Silva, 
et al., 2016). Para guayabo, se recomienda muestrear entre el tercer y quinto nudo a partir 
del ápice, mientras que (Bharat & Das, 2017)2017 indicaron que deben muestrearse hojas 
de cuatro a ocho meses de edad de la parte media de los brotes vegetativos. Por otro lado, 
los resultados obtenidos se fundamentan en los valores óptimos de nutrimentos, pero sin 
considerar las condiciones de clima, suelo y fuente nutrimental con la que se generaron los 
intervalos, por lo que estas recomendaciones pueden variar entre sitios y autores. 
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Referencia N% P% K% Ca% Mg% 





1 1.41-1.65 0.10-0.17 0.51-0.97 1.16-2.12 0.31-0.51 
2 1.42-2.31 0.092-0.134 0.32-0.66 1.70-3.40 0.19-0.32 
3 1.69-2.19 0.168- 0.236 1.20- 1.67 0.60-1.27 0.35-0.50 
4 1.8-2.0 0.12-0.16 1.46-2.08 1.13-1.69 0.25-0.31 





1= (Hundal et al., 2007); 2=(Beyhan, Bozkurt, & Boysal, 2011); 3= Anjaneyulu & Raghupathi, 2009; 
4=(Kotur, Ramkumar, & Singh, 1997). 


ngaadar Continuación cuadro anterior. 





Referencia S% Fe mg/kg Mn mg/kg Zn mg/kg Cu mg/kg 
1 0.18-0.28 105-153 58-110 15-29 6-16 
2 - 70-148 18-63 6.70-11.10 1.32-2.88 
3 0.29-0.47 114- 178 34-77 29-41 6-12 


1= (Hundal et al., 2007); 2=(Beyhan, Bozkurt, & Boysal, 2011); 3= Anjaneyulu & Raghupathi, 2009. 


Tabla 8. Rangos de suficiencia nutrimental derivados de las normas DRIS reportados por varios 
autores. 





F=2.39 
p=0.035 
S=0.328309 
R2=41.71 


AB 


AB 


Calcio% 





AFF AFS AFSF T TéCF  TéCS  TéCSF  TéEF TéES  TéESF 
Tratamientos. 


Figura 6. Contenido de Calcio en las hojas. 





Tratamientos PY% Ca%  Mg% K% Femg Cumg/ Mnmg/ 
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AFFS 
Testigo 


(6) 
(6) 


(o) 
(6) 


B 
B 


B 
B 


(6) 
B 


B 
B 


(6) 
(6) 





B=bajo (valor más pequefio al límite inferior del rango de suficiencia), D=deficiente (valor en el límite 
inferior del rango de suficiencia), O=óptimo (valor dentro del rango de suficiencia), A=alto (valor en el 


límite superior del rango de suficiencia) y E=excesivo (valor superior al rango de suficiencia). 


Tabla 9. Efecto de los abonos líquidos en el balance nutrimental en hojas de guayabo (P. guajava). 








Trat. FH MH FA FB LP PH LH AH L/A L*A 
TéEF 1.65ab 3.63a 3.9ab 1.6bcde 0.6b 1.553a 4a 5.4a 22 a 62.3 a 
TéES 1.75ab 411a 41a 24abcd 0.7b 1.032 d 10.0 cd 4.5c 22 a 46.7 d 
TéEFS 1.56ab 3.90a 41a 1.2e 0.6b 1.201 cd 9.6d 49abc 21a 47.9 cd 
TéCF 1.59ab 3.90a 39ab 2.6abc 0.6b 1.307 bc 10.3 bed 4.9 bc 21a 52.4 bed 
Técs 2.82ab 3.88a 3.8ab 20abcd 0.6b 1.170 cd 10.4 bed 4.8 bc 23a 50.2 cd 
TéCSF 2.10a 3.65a 34b 1.7bcde 1.7a 1.481 ab tab 5.2 ab 2la 59.1 ab 

AFF 2.25b 356a 39ab 2.6 ab 0.7 b 1.245 cd 10.8abc 5.1ab 22a  55.8abc 

AFS 1.71ab 3.80a 3.5ab 2.8a 0.6b 1.218 cd 10.7 abc  4.9bc 22a 52.9bcd 
AFFS 1.65ab 3.61a 3.6ab 1.5cde 0.8ab 1.236cd 10.4bcd 4.7 bc 22 a 50.6 cd 
Testigo 2.14a 3.68a 34b 1.4 de 0.7b 1.363abc 10.9ab 5.0 ab 22a 55.7 abc 





FH=Forma de hoja, (1. Elíptica 2. Oblonga, 3. Lanceolada, 4. Oval, 5. Obovada, 6. Trapezoide); 
MH=Margen de hoja (1=Ondulada, 2= Sinuada, 3=Equilibrada, 4=Recta o plana, 5=Ligeramente 


curva, 6=Doblada); FA=Forma del ápice (Obtusa, 2. Apiculada, 3. Acuminada, 4. Aguda, 5. Redonda); 


FB=Forma de la base 1= Redonda, 2= Oblicua, 3=Aguda, 4=Atenuada, 5=Cordada), LP=Largo del 


pecíolo (cm). PH=peso de hoja (9), LH=Longitud de hoja (cm), AH=Ancho de hoja (cm), L/A=Relación 
largo/ancho de hoja, L*A= Area foliar (cm?). 


Variables morfológicas en hojas de guayabo. 


Tabla 10. Efecto de los abonos líquidos en la forma de las hojas de guayabo. 


La forma de la hoja (FH) fue elíptica en los tratamientos: Te de estiércol foliar (TéEF), 


Té de composta Foliar (TéCF) y abono de frutas foliar (AFF) y Oblonga en los tratamientos: 


Té de estiércol al suelo (TéES), Te de estiércol Foliar y suelo (TéEFS), Te de composta 


al suelo (TéCS), Te de composta suelo y Foliar (TéCSF), Abono de frutas al suelo (AFS), 


Abono de frutas foliar y Suelo (AFFS) y Testigo (sin abono). El margen de hoja (MH) fue 


plano en todos los tratamientos. La forma del ápice (FA) fue aguda. La forma de la base 


(FB) redonda y el largo del pecíolo (LP) resultó inferior a 1 cm. El peso de hojas (PH) varió 
en el rango de 1.032 g (TéES) a 1.553 g (TéEF). La longitud de la hoja (LH) fue grande 
(> a 10 cm), existiendo correspondencia con el ancho de la hoja (AH), el cociente de la 


relación largo-ancho (L/A) de hoja se tomó como indicador de la forma de la hoja, quedando 


en la categoría de mediana (mayor o igual a 1.5 cm y menor o igual a 3.0 cm) resultando 


mediana en todos los casos. El área foliar se estimó multiplicando el largo por ancho (L*A), 


resultando que quedaron comprendidas en el rango de 46.7 (TéES) a 59.1 cm? (TéCSF). 
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Tamafio del fruto 


Para la clasificación por tamafio del fruto se usó la siguiente categorización: Extra (O 
>5.3 cm), Primera (O >4.3y O <5.3 cm.), Segunda (0 >3 cm., y O <4.2 cm.) y Tercera 
(O <3 cm.). En el cuadro 7, se aprecia que los tratamientos Té de estiércol foliar (TéEF), 
Té de composta suelo y foliar (TéCSF) y abono de frutas foliar (AFF), produjeron frutos de 
tamafio extra. Los tratamientos: abono de frutas filiar y suelo (AFFS), Testigo (sin abono), 
abono de frutas al suelo (AFS) y Té de estiércol foliar y suelo (TéEFS) produjeron frutos de 


tamafio primera. 








Tratamiento P/F (9) DP (cm) DE (cm) Categoría 
TéEF 100.1 a 6.250 a 5.525 a Extra 
TéES ND ND ND ND 

TéEFS 58.6 a 5.225 a 4.694 ab Primera 
TéCF 38.0 a 3.800 a 3.700 b Segunda 
Técs ND ND ND ND 
TéCSF 116.9a 6.771 a 5.786 a Extra 
AFF 108.5 a 6.612a 5.738 ab Extra 
AFS 62.7 a 6.300 a 4.733 ab Primera 
AFFS 740a 5.633 a 5.067 ab Primera 
Testigo 73.3 a 5.793 a 4.962 ab Primera 





Té de Estiércol Foliar (TéEF) Té de Estiércol Suelo (TéES); Té de estiércol Foliar y Suelo (TéEFS); Té 
de Composta Foliar (TéCF; Té de composta Suelo (TéCS; Té de Composta Suelo y Foliar (TéCSF); 
Abono de Frutas Foliar (AFF); Abono de Frutas Suelo (AFS); Abono de frutas Foliar y Suelo (AFFS y 

Sin abono (Testigo). Peso por fruto P/F (9), diámetro polar DP (cm), diâmetro ecuatorial DE (cm). ND= 

No determinado. 


Tabla 11. Efecto de los abonos líquidos en el tamafio del fruto de guayaba. 
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RESUMEN: La planta de tomillo es una hierba 
aromática que puede llegar a crecer a una altura 
máxima de 50 cm. En la actualidad tiene distintos 
usos, en fresco o en seco, siendo uno de los 
principales como condimento gastronómico y 
la extracción de su aceite esencial. El objetivo 
de la investigación fue realizar un análisis de 
crecimiento y acumulación de biomasa de 
plantas de tomillo, mismas que se establecieron 
en condiciones de ambiente natural en la parcela 
orgánica del campo agrícola de la Universidad 
Autónoma de Baja California Sur. Para el estudio 
se realizaron 6 muestreos con intervalos de 
tiempos de 20 días. El experimento se desarrolló 
en un suelo Yermósol Háplico (FAO-UNESCO, 
1999). De acuerdo a la clasificación climática 
de Kôppen, modificada por García (1981) para 
condiciones de México, la ciudad de La Paz, 
presenta un clima BW(h” h w (e), es decir, seco 
desértico, cálido. Los parámetros del crecimiento 
evaluados fueron el área foliar (AF), índice de 
área foliar (IAF), tasa relativa de crecimiento 
(TRC), tasa de asimilación neta (TAN) y de masa 
seca y área de proyección. Se realizaron análisis 
estadísticos mediante el Statistica 6.0, con base 
en un disefio de bloques completos al azar. 
Se determinaron las variables de crecimiento 
teóricas, donde el IAF alcanzó un crecimiento 
determinado a los 70 días, por su parte la TAN 
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máxima se alcanzó a los 90 días con 8.8 g cm? dia!, a su vez la TRC, logró su máxima 
expresión a los 90 días, por lo que la máxima ganancia de biomasa en el área de superficie 
ocupada por la planta fue de 7 g dia! a los 90 días, mientras que la máxima expresión del 
cociente del área foliar se logró a los 50 días, a asimismo el cociente de peso foliar del 
peso de se incremento hasta los 90 días, por otro lado la máxima relación vástago/raíz se 
logró a los 90 días, mientras que la DAF logró su máxima expresión a los 90 días después 
del trasplante, por lo que podemos concluir que el tomillo por tener su sistema de vástago 
permaneció en condición de ambiente natural, en un lugar donde los factores y elementos del 
clima que intervienen en su crecimiento oscilan en rangos diferentes. 

PALABRAS CLAVE: Zonas áridas, biomasa, crecimiento, producción. 


EFICÁCIA AGROBIOLÓGICA DO FEIJÃO DOLICHOS (Lablab purpureus L.) NA 
CULTURA DE TOMILHO (Thymus vulgaris L.) 


RESUMO: O tomilho é uma erva aromática que pode atingir uma altura máxima de 50 cm. 
Atualmente tem diversos usos, frescos ou secos, sendo um dos principais como condimento 
gastronômico e na extração de seu óleo essencial. O objetivo da pesquisa foi realizar uma 
análise do crescimento e acúmulo de biomassa de plantas de tomilho, que foram estabelecidas 
em condições ambientais naturais na parcela orgânica do campo agrícola da Universidade 
Autônoma de Baja California Sur. Para o estudo, foram retiradas 6 amostras com intervalos 
de tempo de 20 dias. O experimento foi desenvolvido em solo de Yermósol Háplic (FAO- 
UNESCO, 1999). De acordo com a classificação climática de Kôppen, modificada por García 
(1981) para as condições do México, a cidade de La Paz tem um clima BW(h” h w(e), ou 
seja, desértico quente e seco. Os parâmetros de crescimento avaliados foram área foliar 
(AF), índice de área foliar (IAF), taxa de crescimento relativo (TRC), taxa de assimilação 
líquida (TAN) e massa seca e área de projeção. As análises estatísticas foram realizadas 
usando Statistica 6.0, com base em um delineamento de blocos completos ao acaso. Foram 
determinadas variáveis teóricas de crescimento, onde o IAF atingiu um certo crescimento aos 
70 dias, enquanto o TAN máximo foi atingido aos 90 dias com 8,8 g cm-2 dia-1, por sua vez o 
CRT, atingiu sua expressão máxima aos 90 dias, então o ganho máximo de biomassa na área 
de superfície ocupada pela planta foi de 7 g dia-1 aos 90 dias, enquanto a expressão máxima 
da razão de área foliar foi alcançada aos 50 dias, da mesma forma a razão de peso de folha 
do peso de aumentou até 90 dias, por outro lado, a relação caule/raiz máxima foi alcançada 
aos 90 dias, enquanto o DAF atingiu sua expressão máxima aos 90 dias após o transplante. 
Pelo que podemos concluir que o tomilho, devido ao seu sistema de caule, manteve-se em 
condição de ambiente natural, num local onde os fatores e elementos do clima que intervêm 
no seu crescimento oscilam em diferentes amplitudes. 

PALAVRAS-CHAVE: Zonas áridas, biomassa, crescimento, produção. 


AGROBIOLOGICAL EFFECTIVENESS OF DOLICHOS BEANS (Lablab purpureus 
L.) IN THE CROP OF THYME (Thymus vulgaris L.) 


ABSTRACT: Plant Thyme is an aromatic herb that can grow to a maximum height of 50 cm. 
Currently, it has various uses, in fresh or dry, is one of the main and gourmet seasoning and 


Investigación, tecnología e innovación en ciencias agrícolas Capítulo 3 FEJES 


extraction of essential oil. The research aimed to analyze the growth and accumulation of 
biomass plant thyme, the same that were established in conditions of the natural environment 
in the organic plot of an agricultural field. O The experiment was conducted in a soil yermosol 
Háplico (FAO-UNESCO, 1999). According to the Koppen climate classification, modified by 
García (1981) for conditions of Mexico, the city of La Paz, it presents a BW (h 9 h w (e), 
ie desert dry, warm weather. Growth parameters evaluated were leaf area (AF), leaf area 
index (LAI), relative growth rate (TRC), net assimilation rate (NAR), and dry mass and 
projection area. Statistical analyzes were performed using Statistica 6.0, based on a design 
of a randomized complete block. the variables of theoretical growth, where the IAF reached 
a certain growth at 70 days, by your TAN maximum part was reached at 90 days with 8.8 g 
cm-2 day-1, turn the TRC, it determined achieved its maximum expression at 90 days, so the 
maximum gain of biomass in the surface area occupied by the plant was 7 g day-1 to 90 days, 
while the maximum expression of the ratio of leaf area was achieved at the 50 days, also 
the ratio of leaf weight was increased to 90 days, on the other hand, the maximum stem/root 
ratio was achieved at 90 days, while the DAF achieved its maximum expression at 90 days 
after transplant, so we can conclude that thyme to have their stem system condition remained 
natural environment, in a place where the climate factors and elements involved in its growth 
vary in different ranges. 

KEYWORDS: Arid zones, biomass, growth, production. 


INTRODUCCIÓN 


El tomillo (Thymus vulgaris L.) es una especie perteneciente a la familia Lamiaceae, 
arbusto enano perenne (Omidbaigi y Arjmandi, 2002) cuya parte útil son las hojas y tallos 
(Naghdi Badi et al., 2004). 

El uso de las plantas aromáticas ofrece agradables aromas y sabores en los 
alimentos, alivian el dolor y curan enfermedades (Craker, 2007). Este interés se debe a 
las características aromáticas, terapéuticas y de conservación (Sangwan et al., 2001). Las 
plantas aromáticas son de gran importancia económica y comercial para el mundo. La 
planta de tomillo tiene gran demanda para uso culinario tanto en fresco como en seco, 
así como la utilización de sus aceites esenciales en la industria (Acevedo et al., 2004). La 
importancia económica de las plantas aromáticas se debe a la generación de capital por su 
producción y mercado a nivel mundial y nacional (CCI, 2007). 

El conocimiento de las plantas medicinales y aromáticas se extiende a cualquier 
parte del mundo donde el hombre tradicionalmente ha necesitado de ellas para curar sus 
dolencias o enfermedades. Las diferentes culturas a través de la historia han creado todo 
un conocimiento de remedios vegetales que ha constituido la base de la medicina moderna 
(Font Quer, 1979). Considerando que el cultivo de plantas aromáticas y medicinales supone 
una alternativa a los cultivos tradicionales, es necesario poner a punto las técnicas de cultivo 
y transformación de aquellas especies que tienen una mayor demanda en el mercado. En 
las especies del género Thymus su importancia reside en el contenido elevado que poseen 
en ácidos fenólicos y flavonoides (Martínez, 1980). 


Investigación, tecnología e innovación en ciencias agrícolas Capítulo 3 EEN 


En México, el tomillo se cultiva en los Estados de Puebla, San Luis Potosí y Baja 
California Sur, con una superficie cultivada de 11.25 ha (2.0, 4.5 y 5.25 ha, respectivamente) 
con rendimientos que no superan las 2.68 t ha! de producto seco (SIAP, 2009). Datos 
recientes sefialan que Baja California Sur tiene una superficie sembrada de tomillo de 3 ha 
con un rendimiento de 5.67 t ha”! de producto fresco (SIAP, 2012). 

La utilización de leguminosas como 

Los factores que afectan la producción de biomasa y la acumulación de aceites 
volátiles en las especies aromáticas (McConkey et al., 2000) son los ambientales y la 
variación geográfica (Figueriredo et al., 2008), así como la época de cosecha, la edad de la 
planta y la densidad del cultivo (Kaloustian et al., 2005). 

Por otro lado, los eventos fotosintéticos se desencadenan cuando la clorofila capta 
la energía de la luz incidente en el aparato fotosintético, la energia utilizada para lograr 
la excitación puede llevar a la oxidación de los pigmentos fotosensibles ópticos. Las 
reacciones fotoquímicas que ocurren allí han sido bien estudiadas, tanto in vivo como in 
vitro, y ahora tenemos una idea bastante precisa de los fenómenos que ocurren. Salisbury 
(2000) establece que las tasas fotosintéticas relativas de algunas gramíneas herbáceas y 
especies de hoja ancha se pueden definir en función de la irradiancia que cae por unidad 
de área de la hoja. Teóricamente, sin embargo, la eficiencia fotosintética de las plantas 
nunca debería superar el 18% de la radiación fotosintética activa (RAF) absorbida, y hasta 
la fecha no existe ningún mecanismo que pueda aumentarlas. 

Hunt (2003) muestra que el crecimiento, en el contexto de plantas individuales, 
puede definirse como un aumento irreversible en el tiempo y que estos cambios pueden 
ser de tamafio, forma y cantidad; Asimismo, afirma que, en general, los análisis de 
crecimiento realizados con materia seca y que las plantas anuales y perennes cultivadas 
en condiciones normales muestran, en sus etapas iniciales, la misma curva típica de 
crecimiento. Sin embargo, el cambio en la materia seca total durante estos períodos puede 
variar considerablemente debido a varios factores y se pueden observar pequefios cambios 
durante estos períodos. Finalmente, el autor concluye que una transformación logarítmica 
de los datos puede ser beneficiosa para el investigador, desde el punto de vista visual, 
estadístico y matemático. 

Raffo e Iglesias (2004 ) demostraron que la radiación fotosintética activa absorbida 
por las plantas determina la cantidad de materia seca producida y la calidad de la 
producción, afectando el crecimiento y el rendimiento. Sin embargo, y aunque los efectos 
de los diferentes niveles de radiación en el tomillo no han sido reportados de manera 
comparable a este trabajo, se han realizado investigaciones en otras especies como la 
habichuela (Jarma et al. 1999), que concluyeron que la tasa de crecimiento, tasa de 
absorción neta (TAN), área foliar específica (AFE), relación de área foliar (RAF) y área 
foliar (AF)) mostraron valores mayores en los tratamientos de sombreado, que involucraron 


un tratamiento con radiación incidente del 100%. Por su parte, el índice de cosecha (IC) 
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mostró valores similares para los tratamientos sombreados, aunque de estos, el índice con 
mayor sombreado fue el más efectivo para contribuir al sombreado lleno de semillas. 

Con base en la inexistencia de trabajos que permitan conocer el comportamiento del 
tomillo ante las condiciones edafoclimáticas de Baja California Sur se planteo el objetivo 
de determinar la curva de crecimiento y acumulación de biomasa en plantas de tomillo 
(Thymus vulgaris L.) con fertilización de abonos verdes. 


La investigación se realizó en el ciclo primavera-verano, en el campo agrícola 
experimental de la Universidad Autónoma de Baja California Sur (UABCS) localizado 
en el km 5.5 de la carretera al sur, en la Ciudad de La Paz, ubicada en el Valle de La 
Paz, municipio del mismo nombre, Estado de Baja California Sur, México, situada en las 
coordenadas 24º 10” latitud norte y 110º 19º longitud oeste, a 18.5 metros sobre el nivel del 


mar (Figura 1). 





Figura 2. Plantas de tomillo en la parcela orgánica de la UABCS y toma de datos de peso seco y peso 
fresco en el estudio de análisis de crecimiento y acumulacióin de biomassa. 
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Condiciones de suelo y clima del área experimental 


El experimento se desarrolló en un suelo Yermósol Háplico (FAO-UNESCO, 1999). 
De acuerdo con la clasificación climática de Kôppen, modificada por García (1981) para 
condiciones de México, la ciudad de La Paz presenta un clima BW(h” h w (e), es decir, 
seco desértico, cálido. 

El abono verde fue el frijol Dolichos (Lablab purpúreos), el cual se evaluó desde el 
punto de vista de su utilidad práctica como sustituto de fertilizantes químicos nitrogenados 
y aportes de materia orgánica al suelo, para ello se tuvo en cuenta el criterio seguido 
por diversos autores que definen las especies aptas para estos fines. Para el calculo de 
materia verde, se tomó al azar un metro cuadrado por cada uno de los seis tratamientos y 
posteriormente se determinó el peso mediante una balanza de reloj, finalmente se calculó 
el peso obtenido por hectárea. La incorporación se realizó cortando y depositando sobre 
el terreno los residuos cuando las plantas presentaban alrededor de 30 % de floración, 
esto se estableció tomando como base lo recomendado por Beltrán-Morales et a! (2004) y 
Beltrán-Morales (2006). La cantidad de abono verde aplicada al suelo fue en promedio de 
50 t ha”, mientras que la producción de materia seca del frijol dolichos fue de 13 t ha”. De 
Por otra parte también se determinó la aportación de nutrimentos del frijol dolichos (Cuadro 
1) utilizando la metodologia propuesta por Alcántar y Sandoval (1999). 








Nutrimento Aporte nutricional 
kg ha” 
Ca 211.39 
Mg 75.14 
K 233.80 
Na 123.58 
Fe 16.80 
Mn 2.34 
Zn 0.39 
Cu 0.39 
P 36.22 
cl 208.46 
NO3-N 2.12 
N 364.03 
B 2.62 





Cuadro 1. Aportación de nutrimentos en kg ha” del frijol dolichos (Dolichos lablab). 


Los parâmetros del crecimiento evaluados fueron el área foliar (Af), índice de 
área foliar (IAF), tasa relativa de crecimiento (TRC), tasa de asimilación neta (TAN) y de 
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masa seca y área de proyección. Se realizaron evaluaciones fisiológicas del estado de 
desarrollo del cultivo y se tomaron muestras a diferentes intervalos con un total de 6 fechas 
de muestreo para la biomasa fresca y seca. Para cada análisis se tomaron de la parcela 
experimental diez plantas al azar. 

Las evaluaciones que se realizaron fueron: 


a. Biomasa total: se determinó la biomasa de hojas, tallos y raíz fraccionados 
para cada órgano. Para determinar la biomasa seca de las plantas muestrea- 
das, se secaron colocándose en estufa a 70 ºC durante 48 horas. Se sumaron 
las biomasas (masas secas) de tallos y hojas para determinar la biomasa de la 
parte aérea expresada en g planta”. La Figura 2, muestra algunas plantas de 
tomillo a las cuales se les realizó el análisis de crecimiento y acumulación de 
biomasa. 


b. Área foliar expresada en cm?. Para esta determinación se tomaron diez plan- 
tas de cada parcela y a cada una de las hojas se les calculó el área según la 
metodologia establecida por Ruiz et al. (2007). 


c. Índice de área foliar: se determinó como la relación entre la suma del área 
foliar de las diez plantas, lo cual se realizó cada 15 días. El área de proyección 
se midió sobre el área que las plantas ocupan. La masa fresca del tallo o peso 
fresco se le realizó a cada planta separando la parte de tallo y hojas de la raíz, 
expresada en g planta”, las demás variables se realizaron teóricamente con la 
metodologia propuesta por Coombs et al. (1998). En la Figura 3 se observa la 
forma en que se tomaron las lecturas del peso fresco y seco de las plantas de 
tomillo. 


Determinación del análisis fisiológico del cultivo 
Los parámetros se calcularon de la siguiente manera: 
Índice de área foliar (IAF) 
IAF=AF/AP | donde: AF=área foliar (cm?) por planta 
AP= área de proyección (cm?) 
Tasa de asimilación neta (TAN) 


La tasa de asimilación neta (TAN, g cm? día!) se determinó con la ecuación: 


W2-W1, In A2-Ln Ailn A2-Ln Ai 


TAN = AZ-AL t2-t1 t2-t1 


=g cm? día” 
Dónde: W= peso seco total (gr) 

t= tiempo (días, semanas) 

Ln= logaritmo natural 

A= área foliar (cm?) 
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Tasa relativa de crecimiento (TRC) 
La tasa relativa de crecimiento (TRC, g día”) se determinó con la ecuación: 


qro =!"W2=InWi jo, 
= = =gr.ídias 


Dónde: W= peso seco total 
t= tiempo (días, semanas) 
Ln= logaritmo natural 
Tasa de crecimiento del cultivo (TCC) 


La tasa de crecimiento del cultivo (TCC, g día!) se determinó con la ecuación: 


Dónde: W=peso seco total (9) 
t= tiempo (días, semanas) 


Cociente de área foliar (CAF) 


El cociente de área foliar (CAF, cm? día!) se determinó con la ecuación: 
AF : 
CAF=—=cm' /dia 
W 
Dónde: W=peso seco total en g. 
AF=área foliar por planta cm? 
Cociente de peso foliar (CPF) 
El cociente de peso foliar (CPF, g g!) se determinó con la ecuación: 
CPF=PF/W=gr/ gr 
Dónde: W=peso seco total (9) 
PF=peso seco foliar (9) 
Área foliar específico (AFE) 
El área foliar específico (AFE, cm? g'!) se determinó de la siguiente forma: 
AFE=AF/PF=cm'/ gr 
Dónde: AF=área foliar (cm?) 
PF=peso seco foliar (g) 
Relación vástago-raíz (W/r) 


La relación vástago-raíz (W/R, g g') se determinó de la siguiente forma: 
v/r= peso seco vástago/peso seco raíz = g g' 
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Duración del área foliar (DAF) 
La duración del área foliar (DAF, cm?*día) se determinó con la ecuación: 


- (IAFA+ IAF2)+ (t2—t1) 
= a 


DAF =cm'* dia 


Dónde: IAF=índice de área foliar (cm?) 
t= tiempo (días, semana) 
Tasa foliar unitaria (TFU) 
La tasa foliar unitaria (TFU, gr cm2/cm?2 día!) se determinó con la ecuación: 
o (W2 — W1)+ (Ln AF2 — Ln AF1) 


TFU 
(AF2— AF1)+ (t2—t1) 


=gr+*cm'/cm'+ dias 


Rendimiento (R) 


Estimación aproximada ya que AF no es muy exacta 


R = DAF*TFU= gr 


Número, fechas e intervalos de muestreos 


En total se realizaron seis muestreos, mismos que se efectuaron en diferentes 
fechas e intervalos, para un total de 90 días de duración del cultivo establecido en campo 








(Cuadro 2). 
Numero de Fecha de Intervalo Total de días 
muestreo muestreo de tiempo 
1 F1 0 0 
2 F2 20 20 
3 F3 20 40 
4 F4 20 60 
5 F5 20 80 
6 F6 20 100 





Cuadro 2. Número, fechas de muestreos e intervalos de tiempo de muestreo que se realizaron para el 
estudio de la biomasa en la planta de tomillo. 


Análisis estadístico 


Se realizó un análisis descriptivo de los parámetros del crecimiento evaluados, 
empleando el programa STATISTICA 6.0. 
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RESULTADOS 


Peso fresco y seco de tallo y de raíz 


En el Cuadro 3 se observan los valores del peso (9) promedio de biomasa fresca y 


seca del sistema de vástago y sistema radicular de la planta de tomillo de cada muestreo 


de 10 plantas tomadas al azar. Con estos resultados de biomasa se calcularon los índices 


de crecimiento. 








Número de muestreo PFT (9) PST (9) PFR (9) PSR (9) 
1 0.042 0.002 0.001 0.0 
2 0.463 0.084 0.098 0.019 
3 1.961 0.463 0.139 0.060 
4 4.449 1.061 0.312 0.121 
5 4.405 0.945 0.294 0.103 
6 46.581 15.076 1.193 0.69 





PFT= peso fresco de tallo. PST= peso seco de tallo, PFR= peso fresco de raíz. PSR= peso seco de 


Cuadro 3. Promedio de biomasa fresca y seca de tallo y raíz de plantas de tomillo. 


Longitud de tallo y de raíz 


raíz. 


El Cuadro 4 muestra el crecimiento promedio de tallo y raíz de 10 plantas de tomillo 


(Thymus vulgaris L.). 





Número de muestreo Longitud de tallo (cm) 


Longitud de raíz 





(em) 
1 6.69 5.82 
2 12.23 7.53 
3 18.07 8.28 
4 21.55 15.36 
5 21.59 15.36 
6 40.68 15.7 





Cuadro 4. Crecimiento promedio de longitud de tallo y de raíz de plantas de tomillo. 
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Valores promedio de los parámetros de crecimiento 


En el Cuadro 5 se presentan las fechas de muestreo y los valores promedio 


calculados para cada uno de los parâmetros del crecimiento evaluados. 








Fecha IAF TAN TRC TCC CAF CPF AFE VIR DAF TFU R 
F1 0 1.3 0.09 0.0025 0.62 1.5 0.53 2 0.0639 0.079 0.05 
F2 0.00047 3.4 0.02  0.0072 019 1.22 0.33 4.42 0.0591 0.832 0.041 
F3 0.0005 3.36 0.03  0.0476 0.42 1.12 0.48 4.71 0.03824 2.071 0.067 
F4 0.0003 0.099 0.09 0.237 0.26 1.11 0.22 8.83 0.0347 8.809 0.030 
F5 0.00031 2.77 0.39  7.1271 0.5 4.77 0.27 9.17 0.0153 1.791 0.027 
F6 0.00073 1.5 0.05  0.6565 0.5 3.13 0.16 21.84  0.0919 1.567 0.144 





IAF: índice de área foliar (cm'?), TAN: tasa de asimilación neta (g cm? día!), TRC: tasa relativa de 
crecimiento (g día!), TCC: tasa de crecimiento del cultivo (g día), CAF: cociente de área foliar (cm? 
día!), CPF: cociente de peso foliar (9), AFS: área foliar específico (cm? g!), W/R: relación vástago-raíz 
(9), DAF: duración del área foliar, TFU: tasa foliar unitaria (g cm? día”!), R: rendimiento (9). 


Cuadro 5. Fechas de muestreo y valores promedio para cada parámetro del crecimiento evaluado. 


DISCUSIÓN 


Índice de área foliar (IAF) 


El índice de área foliar es un importante atributo del crecimiento aéreo de las plantas, 
debido a que las hojas son dominantes en el proceso fotosintético y por consiguiente para la 
producción de asimilados (Nobel, 1999). EI IAF es la sumatoria de todas las áreas de hojas 
por una unidad de suelo y la estimación directa consiste en la cosecha de material vegetal 
en un área conocida (Wilhelm et al., 2000), pero también se determina indirectamente por 
medio de la intercepción de la radiación solar (Cortés, 2003). 

EI IAF, el cual alcanzó un crecimiento determinado a los 50 días, pero en la última 
fecha se logró un número máximo de IAF, esto debido a que la plantación se le realizó 
una poda y se dejaron pasar 20 días a la fecha de lectura anterior para luego volverse a 
muestrear. 

Tasa de asimilación neta (TAN) 


La TAN también denominada tasa foliar unitaria (TFU) o índice de eficiencia 
fotosintética, representa la ganancia neta en peso seco por unidad de área foliar y es una 
medida indirecta de la fotosíntesis (Hunt, 1982). Es definida como el incremento de material 
vegetal por unidad de tiempo y se expresa en g'm?dia! (Beadle, 1988). La TAN es una 
medida de la eficacia del follaje, el cual constituye la principal fuente de fotoasimilados en 
la producción de materia seca e indica también la velocidad de fotosíntesis neta en un lapso 
relativamente largo, entre dos muestreos (Escalante y Kohashi, 1993). 

Aa los 90 días se logró la máxima TAN (8.8 g cm? dia), posteriormente esta tiende 
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a disminuir por el envejecimiento de la planta, el cual es un resultado que se compara con 
lo expresado por Gardner et al. (1990) quien menciona que este parámetro no es constante 
con el tiempo y muestra una tendencia a disminuir con la edad de la planta. La disminución 
se acelera por un ambiente desfavorable y la ganancia de materia seca por unidad de área 
foliar disminuye en la medida que brotan nuevas hojas, debido al sombreo reciproco. Por 
su parte, Hunt (1982) menciona que la TAN, indica la eficiencia fotosintética promedio, 
individual o en una comunidad de plantas. La capacidad de la planta para incrementar 
su masa seca en función del área asimiladora en periodos cortos a lo largo del ciclo de 
crecimiento depende del área foliar, de la disposición y edad de las hojas y de los procesos 
de regulación interna relacionados con la demanda de los asimilados. 


Tasa relativa de crecimiento (TRC) 


La TRC, la cual logró su máxima expresión a los 90 días. Según Pedroza et al. 
(1997) esta es un índice de eficiencia que expresa el crecimiento en términos de una tasa 
de incremento en tamafio por unidad de tamafio y tiempo; lo que representa la eficiencia 
de la planta como productor de nuevo material y depende de la fotosíntesis total y de la 
respiración (Sivakumar y Shaw, 1978); además se propone como una medida que integra 
el comportamiento fisiológico de las plantas (Radford, 1967). Es una medida del balance 
entre la capacidad potencial de fotosíntesis y el costo respiratorio (Archila et a!., 1998), pero 
además expresa el incremento en masa seca de la planta en un intervalo de tiempo dado, 


tomando como referencia el valor inicial de la masa seca producida y acumulada. 
Tasa de crecimiento del cultivo (TCC) 


Según Hunt (1982) la TCC mide la ganancia de biomasa vegetal en el área de 
superficie ocupada por la planta. Es aplicable a plantas que crecen juntas en cultivos 
cerrados. La máxima TCC ocurre cuando las plantas son suficientemente grandes o densas 
para explotar todos los factores ambientales en mayor grado. En ambientes favorables, la 
máxima TCC ocurre cuando la cobertura de las hojas es completa y puede representar el 
máximo potencial de producción de masa seca y de tasas de conversión en un momento 
dado (Brown, 1984). De acuerdo a lo expuesto anteriormente, los resultados encontrados 
en tomillo, se puede decir que la máxima ganancia de biomasa en el área de superficie 
ocupada por la planta fue de 7 g dia” a los 90 días. 


Cociente de área foliar (CAF) 


El área foliar y el índice de área foliar son parámetros ampliamente utilizados en 
estudios de ecofisiologia de cultivos (Coombs y Hall, 1982). La máxima expresión del 
cociente del área foliar se logró a los 50 días, para luego lograr la disminución del cociente 
de área foliar en tomillo. Estos resultados concuerdan con lo expuesto por Kvet y Marshall 
(1971), Astegiano y Favaro (1987) quienes mencionan que las funciones que relacionan 
estas dimensiones varían claramente entre especies, aunque también afecta el estudio de 
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desarrollo, genotipo, estación de crecimiento y edad de las hojas. 
Cociente de peso foliar (CPF) 


EL cociente de peso foliar, donde al nivel del peso de las plantas se tuvo un 
incremento hasta los 70 días y después disminuyó. Este se considera la relación entre el 
área foliar total de la planta y la materia seca total de las hojas de la planta. En condiciones 
normales y para la mayoría de los cultivos, el espesor de las hojas aumenta con la edad, 
en el caso del frijol aunque estos valores fluctúan, permanecen dentro de un rango en el 
cual el aumento en espesor es mínimo; sin embargo, esta parece ser la situación normal 
para plantas de frijol desarrolladas en condiciones de soluciones nutritivas. Jones (1968) 
cita valores de AFP para frijol con condiciones de campo y de invernadero; de acuerdo con 
este autor, el AFP de las plantas desarrolladas en soluciones nutritivas es hasta tres veces 
mayor que el de plantas desarrolladas en el campo lo cual se debe a la exuberancia del 
crecimiento en estas condiciones (Jones, 1968). 

Las variaciones en el CPF tiene gran influencia en el índice de crecimiento relativo 
(ICR) ya que debido a la relación que existe entre el ICR y la relación de área foliar (RAF) 
y entre esta y el AFE, cualquier cambio en el área foliar especifica hará variar el índice de 
crecimiento relativo a través de la razón de área foliar (Ascencio, 1972). 


Área foliar especifica (AFE) 


La respuesta del área foliar específica (AFE) en función del tiempo. Durante todo el 
periodo evaluado, la máxima AFE se presentó al inicio del trabajo; lo anterior sugiere que al 
considerar la expansión del área foliar como una respuesta plástica altamente influenciada 
por la radiación, es probable que las hojas más pequefias y delgadas observadas en los 
niveles donde alcanza menor radiación sean la respuesta de un autoajuste de la planta 
para captar el máximo posible de energía solar disponible para la asimilación clorofílica; 
sin embargo, esta situación no representó una mayor eficiencia en la producción de 
fotoasimilados por unidad de área foliar, tal como se demostró con la TAN. 

Según Páez et al. (2000) indican que las plantas pueden ajustarse a un ambiente 
de menor irradiación aumentando el área foliar específica (AFE); así, en el caso del tomate 
estos autores observaron este efecto solamente hasta los 45 días de crecimiento vegetativo. 
Una tendencia similar se observó en investigaciones realizadas por De Grazia et al. (2001) 
en lechuga (Lactuca sativa L.) en las que el AFE fue afectada por la atenuación del nivel de 
radiación, provocando un incremento en ella. Por su parte, Flórez et al. (2006) afirma que 
el AFE, es una medida de la superficie foliar de la planta en términos de densidad o grosor 
relativo de la hoja. Se define como la relación entre el área total de la hoja y la masa del 
área foliar de la planta. 


Relación vástago-raíz (v/r) 


La máxima relación vástago/raíz se logró a los 90 días, esta máxima relación se 
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da debido a que la última lectura se realizó 20 días después de la penúltima, esto quiere 
decir que se dejó pasar el tiempo para determinar si realmente se incrementa la relación 
vástago raíz. Según Lincoln et al. (2006) la relación de la biomasa vástago-raíz y la 
fotosíntesis realizada por el vástago, menciona que un brote crecerá tanto como para que 
el factor limitante sea la cantidad de agua incorporada por las raíces; por el contrario, las 
raíces crecerán hasta que su demanda de fotoasimilados al vástago iguale el aporte. Este 
equilibrio funcional cambia si se reduce el aporte de agua. 

Según Helms y Baker (1982) la relación de la biomasa vástago-raíz parece estar 
gobernada por el vástago. De forma sencilla, un brote crecerá tanto como para que el factor 
limitante sea la cantidad de agua incorporada por las raíces; por el contrario, las raíces 
crecerán hasta que su demanda de fotoasimilados al vástago iguale el aporte. 

Por su parte, Hawley y Smith (1972) mencionan que la expansión foliar se ve 
afectada rápidamente al reducir el aporte de agua, mientras que la actividad fotosintética se 
ve mucho menos afectada. La inhibición de la expansión reduce el consumo de carbono y 
energia y una gran parte de los productos asimilados por las plantas pueden ser distribuidos 
por el sistema radicular, donde pueden sustentar un posterior crecimiento, al mismo tiempo 
los ápices radiculares pierden turgencia en suelos secos. 


Duración del área foliar (DAF) 


La DAF es un área, las unidades de cualquier área son el producto de sus 
dimensiones. Puesto que el IAF no posee unidades, la DAF determinada mediante 
integración se expresa en unidades de tiempo, usualmente días o semanas (Robert et al., 
2004). La máxima expresión alcanzada a los 90 días después del trasplante, esto quiere 
decir que se logró la máxima duración en cm? por día. Esta duración del área foliar, aumentó 
en las primeras fechas y para los 90 días disminuyó, quizás por el aumento de las altas 
temperaturas. 

Según Rodríguez y Leihner (2006), la DAF ha demostrado ser un criterio valioso 
para identificar genotipos productivos. La justificación para ello es que las plantas que 
logran mantener sus hojas durante un mayor número de días requieren una menor 
inversión de energía en el mantenimiento del follaje activo. Por el contrario, los genotipos 
con hojas que caen en poco tiempo, deben restituir el área foliar continuamente, lo cual 
implica un costo biológico que repercute negativamente sobre el rendimiento. En el caso 
de la yuca, por ejemplo, los genotipos más productivos son aquellos que mantienen sus 
hojas por más de 100 días. En este sentido, Ginson y Fische (2004) mencionan que las 
variedades tradicionales de porte alto y de larga duración tienen muchas veces, una 
eficiencia fotosintética (o una tasa foliar unitaria) más alta y una duración de la materia 
seca total mayor que las variedades modernas semi-enanas; éstas son de corta duración 


pero su rendimiento es alto. 
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Tasa foliar unitaria (TFU) 


El peso es una de las formas más adecuadas para expresar el rendimiento parcial o 
total de una planta, ya que es una cuantificación más precisa de la producción por ser una 
medida directa. Además, se pueden establecer comparaciones de rendimiento entre sitios 
diferentes, principalmente cuando se expresa en un peso seco al horno (Gutiérrez, 1984). 

La tasa foliar unitaria, que tuvo un constante aumento foliar, logrando 1600 g cm? dia”; 
esta máxima duración se logró a los 90 días, aunque a partir de los 70 días se mantuvo similar 
hasta los 90 días del último muestreo. 


Rendimiento (R) 


Según Challa et al. (1995) el rendimiento de un cultivo viene dado por la capacidad de 
acumular biomasa (materia fresca y seca) en los órganos que se destinan a la cosecha y un 
incremento proporcional de la biomasa destinada a estos órganos garantiza un incremento 
del rendimiento. Así, la distribución de materia seca entre los diferentes órganos de la 
planta tiene un papel fundamental en la producción de un cultivo. Los resultados obtenidos 
demuestran que lo establecido por Challa et al. (1995) se corroboran en tomillo. Como en 


el transcurso de la evaluación se logró a los 90 días el máximo rendimiento. 


CONCLUSIONES 


* La planta de tomillo cuenta con una facilidad de poder adaptarse en una zona 
donde existe una baja humedad en el suelo. 


* | Puesto que en tomillo también su sistema de vástago permaneció en condi- 
ción de ambiente natural, en un lugar donde los factores y elementos del clima 
que intervienen en su crecimiento oscilan en rangos diferentes. Realmente esta 
planta no quiere de gran demanda de agua. 


* | Para obtener una buena producción de la planta de tomillo de forma orgánico, 
se requiere de contar con un suelo con un alto contenido de materia orgánica 
y/o incorporar de alguna otra fuente donde se obtenga una fertilización de ma- 
nera orgânica para el crecimiento y la producción de esta planta. 


* | De las variables de crecimiento teóricas, se encontró que el IAF alcanzó un 
crecimiento determinado a los 70 días. 


*— La TAN máxima se alcanzó a los 70 días con 8.8 g cm? dia", mientras que la 
TRC, logró su máxima expresión a los 90 días, con una ganancia por la planta 
de 7 g dia! a los 90 días. 


* El máximo cociente del área foliar se logró a los 50 días, a asimismo el cociente 
de peso foliar del peso de las plantas logró un incremento hasta los 70 días, por 
otro lado la máxima relación vástago/raiz se logró a los 90 días. 


* La DAF logró su máxima expresión a los 90 días después del trasplante, 
por lo que podemos concluir que el tomillo por tener su sistema de vástago 
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permaneció en condición de ambiente natural, en un lugar donde los factores 
y elementos del clima que intervienen en su crecimiento oscilan en rangos 
diferentes. Realmente esta planta no quiere de gran demanda de agua por lo 
que la máxima área foliar se logra a los 75 días. 
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RESUMEN: Una nueva línea maternal de 
conejos resiliente seleccionada por variabilidad 
del tamafo de camada ha sido creada en la 
Universidad Miguel Hernández de Elche (línea 
A). En cunicultura, el gazapo de engorde se 
obtiene del cruzamiento entre líneas maternales 
con líneas paternas seleccionadas por caracteres 
de crecimiento. El objetivo de este trabajo es 
comparar el tamafio de camada y los caracteres 
ponderales de los gazapos al predestete en 
animales de la línea A y en animales cruzados. 
Un total de 156 conejas de la línea A fueron 
utilizadas. Las hembras fueron apareadas 
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con machos de la misma línea A y machos 
de una línea seleccionada por velocidad de 
crecimiento para obtener los gazapos cruzados. 
Los caracteres estudiados fueron el número de 
gazapos nacidos totales y vivos al nacimiento y 
a los 21 días. El peso de la camada total y la de 
los gazapos vivos al nacimiento y a los 21 días 
de edad, también se calculó el peso individual. 
Los análisis estadísticos se realizaron con 
metodologia bayesiana. El tamafio de camada y 
el peso individual de los gazapos fue similar en 
ambos grupos de animales. Sin embargo, el peso 
de la camada fue superior para los animales 
cruzados tanto para el total de la camada con una 
diferencia de 29 g (P=0.99) como para la camada 
nacida viva (46 9; P=0.99). El peso de la camada 
a los 21 días tiende a ser superior en los animales 
cruzados (1129; P=0.86). En conclusión, el 
cruzamiento de la línea maternal resiliente con 
machos de una línea de crecimiento mejora los 
caracteres ponderales de la camada desde el 
nacimiento a los 21 días de edad. 

PALABRAS CLAVE: Conejo, cruzamiento, peso 
de la camada, prolificidad, resiliencia. 


ABSTRACT: A new maternal line of resilient 
rabbits selected for variability in litter size has 
been created at the Miguel Hernández University 
of Elche (line A). In rabbit farming, the fattening 
kitten is obtained from the crossing between 
maternal lines with paternal lines selected for 
growth traits. The objective of this work is to 
compare the litter size and the weight characters 
of the kits at pre-weaning in line A animals and in 
crossbred animals. A total of 156 rabbits from line 
Awere used. The females were mated with males 
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of the same line A and males of a line selected by growth speed to obtain the crossed kits. The 
characters studied were the number of kits born total and alive at birth and at 21 days. The 
weight of the total litter and that of the live kits at birth and at 21 days of age, the individual 
weight was also calculated. Statistical analyzes were carried out with Bayesian methodology. 
The litter size and individual weight of the kits was similar in both groups of animals. However, 
litter weight was higher for crossbred animals both for the total litter with a difference of 29 g (P 
= 0.99) and for the live-born litter (46 9; P = 0.99). Litter weight at 21 days tends to be higher 
in crossbred animals (1129; P = 0.86). In conclusion, the crossing of the resilient maternal line 
with males of a growth line improves the weight characteristics of the litter from birth to 21 
days of age. 

KEYWORDS: Rabbit, crossbreeding, litter weight, prolificacy, resilience. 


INTRODUCCIÓN 


El avance de la producción animal está ligado a la adaptación de los animales a las 
nuevas condiciones de producción tanto ambientales, producidas por el cambio climático, 
como de manejo debido por ejemplo a la reducción en el uso de antibióticos o la utilización 
de subproductos en la alimentación. Se hace por tanto necesario, disponer de un material 
genético animal que sea resiliente, es decir, que sea capaz de amoldarse a desafios 
ambientales sin mermar su capacidad productiva. 

Para dar solución a este problema, se está Ilevando a cabo un programa de mejora 
genética en conejo por variabilidad del tamafio de camada en la Universidad Miguel 
Hernández de Elche. El criterio de selección es la varianza fenotípica del tamafio de 
camada dentro de hembra, una vez corregido el tamafio de camada por los efectos del 
afio-estación y estado fisiológico de la hembra (Blasco et al., 2017). La línea seleccionada 
para incrementar la variabilidad del tamafio de camada (línea H) presenta más variabilidad 
(4,4 gazapos?) que la línea seleccionada para disminuir la variabilidad (línea A; 2,7 
gazapos”; Blasco et al., 2017). La variabilidad del tamafio de camada está relacionada 
con la capacidad de adaptación de la hembra a los cambios ambientales adversos; es 
decir con su sensibilidad al estrés y a las enfermedades, o lo que es lo mismo en su 
resiliencia. En este sentido, la línea A muestra una menor concentración basal de cortisol 
y una menor respuesta al estrés que la línea H (Argente et al., 2019; Beloumi et al., 2020). 
Se sabe que el estrés crónico conduce a la desregulación del sistema inmune y aumenta 
la predisposición a enfermar (Glaser y Kiecolt-Glaser, 2005). Esto está de acuerdo con 
el menor porcentaje de eliminación involuntaria encontrado en la línea A (Argente et al., 
2019). Por tanto, disminuir la variabilidad del tamafio de camada tiene implicaciones en la 
mejora del bienestar de las hembras. 

Además, Garcia et al. (2018) observaron que la línea A presenta una mayor 
capacidad de movilizar reservas energéticas al parto que la línea H. El parto es un momento 
de gran demanda energética y aquellas hembras que no sean capaces de recuperar sus 
reservas corporales tras el parto tenderán a una menor fertilidad (Fortun-Lamothe, 2006) 
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y a una disminución de su defensa inmune (Castellini et al., 2010), lo que está de acuerdo 
con el menor porcentaje de eliminación y mortalidad encontrado en la línea A (Argente et 
al., 2019). Por lo tanto, la selección para disminuir la variabilidad del tamafio de camada 
produce hembras que manejan adecuadamente sus reservas corporales y gestionan la 
movilización energética correctamente, y en consecuencia presentan menor riesgo de 
morir o ser eliminadas (García et al., 2018; Argente et al., 2019). 

En cunicultura, la producción del gazapo de engorde está basado en el cruzamiento 
de hembras maternales con machos seleccionados por caracteres de crecimiento (García 
y Argente, 2020). Por tanto, para poder ser utilizada las hembras de la línea A como 
hembras maternales en las explotaciones, es necesario conocer su comportamiento en el 
cruzamiento. El objetivo de este trabajo es comparar el tamafio de camada y los caracteres 
ponderales de los gazapos desde el nacimiento hasta los 21 días de edad en animales de 
la línea A y en animales cruzados. 


MATERIAL Y MÉTODOS 


Todos los procedimientos experimentales que implican animales han sido aprobados 
por el comité ético de investigación de la Universidad Miguel Hernández de Elche (Número 
de Referencia 2021/VSC/PEA/0211). 


Material animal y disefio de experimento 


Se utilizaron un total de 156 conejas pertenecientes a la línea de baja variabilidad 
por tamafio de camada (línea A) desde Junio de 2021 hasta Noviembre de 2021. Los 
machos pertenecían a la línea A para producir el gazapo AxA y a la línea N para producir 
el gazapo AxN. La línea N se selecciona por velocidad de crecimiento en el periodo de 
engorde. Las asignaciones de los machos de la línea A con hembras de la misma línea se 
realizaron de forma que no tuvieran en común ningún abuelo en común. 

El manejo reproductivo se realiza en bandas semanales y con montas a 12 días 
post-parto. El diagnóstico de gestación se realiza mediante palpación abdominal a los 12 
días de gestación. Los partos tuvieron lugar a los 31 días de gestación. No se realizaron 
adopciones en el periodo de lactación. Los gazapos fueron destetados a los 28 días de 
edad. 

En el momento del parto, se contabilizaron los gazapos nacidos vivos y nacidos 
muertos, la fecha del parto y el orden del parto. Se pesaron los gazapos vivos y muertos y 
se comprobó si más del 80% de los gazapos habían mamado, mediante la identificación de 
una marca de leche en el estómago de los gazapos que es fácilmente visible (Argente et al., 
1999). A los 21 días de edad, se contabilizaron y se pesaron los gazapos. 

Los caracteres estudiados fueron: 


*- La prolificidad: medida como el número de gazapos nacidos totales, vivos y 
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muertos y el número de gazapos a 21 días. 


* | Peso total de la camada (9): al nacimiento de los gazapos totales y de los ga- 
zapos vivos y a los 21 días 


* | Peso individual del gazapo (9): al nacimiento de los gazapos nacidos totales y 
de los gazapos vivos, y a los 21 días 


Análisis estadísticos 
Para el número de nacidos totales, vivos y muertos se utilizó el siguiente modelo: 
Ya SU + C+ E + OP, + EP + Csa 
Donde Yu 
cruce AxN; E, es la estación (estación templada; estación calurosa); OP, es el orden de 


es el carácter, | es la media, C, es el tipo de animal (A: cruce AxA; N: 


parto (primer parto; segundo parto; tercer parto o más); EP es el efecto permanente de una 
hembra sobre sus partos y Cu ES el error del modelo. 

Para el número de gazapos a 21 días y los pesos de la camadas y pesos individuales 
de los gazapos hasta los 21 días se utilizó el modelo anterior, pero se incluyó los nacidos 
totales como covariable y el efecto del comportamiento maternal con dos niveles (correcto 
e incorrecto). 

Todos los análisis se realizaron utilizando metodologia bayesiana, con aprioris 
planos acotados para todos los parámetros desconocidos. Las distribuciones marginales 


posteriores de las diferencias entre líneas se estimaron usando muestreo de Gibbs. 


RESULTADOS Y DISCUSIÓN 


En la tabla 1 se muestran los parámetros de las distribuciones marginales posteriores 
de las diferencias entre los dos tipos de cruzamientos estudiadas, animales AxA y AxN.. 

El tamafio de camada fue similar entre ambos cruzamientos desde el parto hasta los 
21 días de edad. El número de gazapos nacidos totales y vivos fue similar al encontrado 
en otros estudios de esta misma línea (Agea et al., 2019) y en otras líneas maternales en 
Argelia (Belabbas et al., 2021) y en Arabia Saudi (Al-Saef et al., 2008). Estos resultados 
son esperables, pues el tamafio de camada es un carácter que se adscribe a la hembra, y 
en este caso las hembras de ambos cruzamientos pertenecían a la línea seleccionada por 
homogeneidad del tamafio de camada. 

El peso de la camada fue un 7% superior en la camada total y un 13% superior 
en los nacidos vivos para los animales cruzados (AxN) que los animales puros (AxA) 
(P=0.99). Este mayor peso de la camada en los gazapos AxN tiende a mantenerse a los 
21 días de edad (D = -112 9; P = 0.86). Estos resultados se han obtenido corrigiendo por 
la covariable tamaio de camada, es decir las diferencias encontradas en los dos tipos de 
animales no son debidas a un mayor tamafio de camada, pues se han estimado a igualdad 
de tamafio de camada. Los padres de los gazapos N son machos pertenecientes a una 
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línea seleccionada por velocidad de crecimiento (Agea et al., 2018). Los pesos obtenidos 
al nacimiento son similares a las líneas maternales V y Saudi Gabali (Al-Saef et al., 2008) 
pero son superiores a los 21 días para estas mismas líneas. 

El peso individual de los gazapos fue similar en ambos tipos de animales. El peso 
de los gazapos vivos es de alrededor de 53.5 g y es un peso similar a gazapos de otras 
líneas (Argente et al., 1999). Se estima que los gazapos con pesos superiores a 50 g al 
nacimiento tienen más de un 90% de probabilidad de sobrevivir al nacimiento y más de 
un 85% de probabilidad de sobrevivir a los 4 días después del parto, si han ingerido leche 
en las primeras horas tras el parto. Mientras que estos porcentajes se reducen un 10% 
aproximadamente si los gazapos no han ingerido leche (Agea et al., 2009; Argente et al. 














1999). 
Carácter A N DN HPDoso, P 

= Nacidos totales 8.21 8.72 -0.51 -1.57; 0.53 0.84 
8 Nacidos vivos 7.45 7.89 -0.44 -1.62; 0.73 0.77 
E Nacidos muertos 0.74 0.81 -0.06 -0.86; 0.71 0.57 
a Gazapos a 21 días 6.46 6.99 -0.53 -1.55; 0.55 0.84 
E s Total nacimiento 410 439 -29 -50; -4 0.99 
8 E o Vivos nacimiento 366 412 -46 -883; -6 0.99 

Õ A21 días 2055 2167 “112 -317; 90 0.86 
E Total nacimiento 58.5 54.4 -0.9 -3.6; 2.1 0.71 
e 
E Ss Vivo nacimiento 53.5 53.8 -0.3 -3.8; 2.6 0.57 
ê A 21 días 345 352 -7 -47; 30 0.64 





Õ 


an: Mediana de la diferencia entre animales A y N. HPD95%: Región de alta densidad posterior al 
95%. P: Prob (D>0) cuando D>0 y Prob (D<0) cuando D<o. 


Tabla 1. Parâmetros de las distribuciones marginales posteriores de las diferencias entre animales A y 
N para las características desde el nacimiento al destete. 


CONCLUSIONES 


El cruzamiento de la línea maternal resiliente con machos de una línea de crecimiento 
mejora los caracteres ponderales de las camadas desde el nacimiento a los 21 días de 
edad. 
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RESUMEN: En los últimos veinte afios, el 
Paraguay ha realizado un avance importante en 
la producción de información estadística macros. 
Contrariamente, informaciones sobre aspectos 
microeconómicos y especificamente agrícola no 
han seguido el mismo patrón, siendo una de las 
razones la dificultad para acceder a información 
primaria en terreno. La experiencia muestra que 
es frecuente enfrentar situaciones en la cual 
el agricultor se rehúsa (y con todo derecho) 
a proporcionar informaciones económicas y 
financieras de su finca. Como alternativa se tiene 
la metodologia desarrollada por The International 
Farm Comparison Network, (COLUNGA, ET AL., 
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2009), y la Companhia Nacional de Abastecimento 
(CONAB) (CONAB, 2010), que costea lo que 
denomina una “granja típica”, también conocida 
como “modal, patrón o estândar”. Por otro lado, 
Ghida, (2009); organizó una metodologia para la 
evaluación económica, financiera y patrimonial de 
una finca agrícola. En lo referente a la evaluación 
económica, parten del cálculo del ingreso 
bruto, al cual le van descontando los costos y 
remuneraciones de factores hasta llegar a la 
rentabilidad. La investigación aquí presentada 
buscó contribuir a la generación de información 
microeconómica agrícola, de ahí que su objetivo 
haya sido la evaluación económica de cuatros 
fincas modales. Se concluyó que concluye 
que las fincas cubrieron sus costos directos e 
indirectos y remuneraron mano de obra familiar, 
pero solamente tres hicieron lo propio con capital 
y tierra. Dos fincas tuvieron una rentabilidad 
superior a la inflación de EEUU. 

PALABRAS CLAVE: Finca modal, evaluación 
económica, Paraguay. 


THE PROFITABILITY OF AN 
AGRICULTURAL FARM FROM 
DIFFERENT PERSPECTIVES: THE CASE 
OF THE 2016/2017MARKETIMG YEAR 
ABSTRACT: Over the last twenty vears, 
Paraguay has made an important progress in 
the production of macro statistical information. 
Conversely, information on microeconomic and 
specifically agricultural aspects have not followed 
the same pattern. One of the reasons behind is 
the difficulty to access primary information from 
the field. Experience shows that it is frequent to 
face situations in which the farmer refuses (it is 
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his right) to provide economic and financial information about his farm. As an alternative, there 
is the methodology developed by The International Farm Comparison Network, (COLUNGA, 
ET AL., 2009), and the Companhia Nacional de Abastecimento (CONAB) (CONAB, 2010), 
which analizes a so called “typical farm”, also known as “modal, pattern or standard”. Ghida, 
(2009) organized a methodology for the economic, financial and patrimonial evaluation of an 
agricultural farm. Regarding the economic evaluation, it starts from the calculation of gross 
income, which is discounted by factor costs and remunerations until reaching profitability. The 
research presented here seeks to contribute to the generation of agricultural microeconomic 
information, hence its objective has been the economic evaluation of four modal farms. The 
paper concluded that the farms covered their costs, and the labour factor of production. Three 
farms did the same with capital and land, and two farms outperformed US inflation. 
KEYWORDS: Paraguay, profitability, agriculture. 


INTRODUCCIÓN 


Una finca agropecuaria se concibe como una unidad productiva que mediante la 
combinación de factores de producción clásicos (capital, tierra y trabajo) y, dirigidos por 
el gerenciamiento del propietario (excepto en los casos de grandes empresas), produce 
bienes y servicios para su comercialización y consumo propio (BARBERIS, BONGIOVANNI, 
GILETTA, 2017). Como todo emprendimiento productivo precisa que sus ingresos sean 
superiores a sus egresos. Esta afirmación lleva a la necesidad de evaluar el emprendimiento 
o finca de modo tal a medir el grado en que estos ingresos cubren los costos y generan 
excedentes. La evaluación de una finca o empresa agropecuaria se realiza mediante la 
utilizaciónde indicadores que se agrupan en físicos (o productivos), económicos, financieros 
y patrimoniales (MARTÍNEZ FERRARIO, 1995); COLUNGA, ETAL., 2009; GHIDA, 2009). La 
evaluación económica, especificamente considera además del egreso en efectivo, aquellos 
costos e ingresos no monetizados tales como la depreciación, el costo de oportunidad, 
consumo, cesiones y/0 transferencias internas, diferencia de inventario, producción no 
vendida en elejercicio analizado (stock). GHIDA, (2009) organizó una metodologiapara la 
evaluación económica, financiera y patrimonial de una finca agrícola. En lo referente a la 
evaluación económica, parte del cálculo del ingreso bruto, al cual le van descontando los 
costos y remuneraciones de factores hasta llegar a la rentabilidad. Esta fue la utilizada en 
este trabajo según se describe en el acápite metodologia. La investigación aquí presentada 
tuvo como objetivo la evaluación económica de cuatros fincas modales del Paraguay. 


METODOLOGÍA 

El presente trabajo se considera un estudio de caso de tipo cuantitativo en el que 
se utilizó un disefo no experimenta, intencional no probabilístico. A través del análisis 
descriptivo de un proceso, buscó sentar las bases para investigaciones posteriores El área 
de estudio ubicada en Paraguay fue dividida en cuatro zonas productivas, a saber (i) norte 
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del Departamento de Alto Paraná, denominada región norte (ii) sur del Departamento de 
Alto Paraná, denominada región sur (iii) nordeste del Departamento de Itapuá, denominada 
región centro y, (iv) este del Departamento de Caaguazú, denominada región Campo 9. En 
cada una de ellas y en colaboración con técnicos locales y agricultores, se construyó una 
finca modal. Este tipo de finca es la representación de la unidad productiva de la región 
y del paquete tecnológico, (coeficiente técnico de insumos, maquinarias e implementos, 
servicios,) más común en la citada unidad geográfica (CONAB, 2010). La finca modal en 
cada una de las cuatro regiones tuvo con una superficie determinaday, un nivel tecnológico 
o manejo representativo de la región. Se acordó en las reuniones que las actividades 
mecanizadas (labores agrícolas) serían totalmente tercerizadas, debido a la gran diversidad 
delparque de maquinarias e implementos agrícolas que poseían los agricultores, que hizo 
difícil acordar un costeo representativo en cada zona. Se acordó igualmente que las fincas 
modales sembrarían tres cultivos, soja, trigo, maíz. Si bien existían experienciasy prácticas 
de siembra de abonos verdes y otros rubros, (girasol y canola), la elevada heterogeneidad 
en su manejo intra e interregiones dificultó estructurar una tecnologíamodal, por lo cual se 
optó por no incluirlas. 

Luego de estas exclusiones o ajustes se acordó que cada finca modal siembre en 
toda suextensión soja en primavera, seguido de maiz y trigo en partes iguales para la región 
centro y Campo 9, y con variaciones especificas en la región norte y sur. Esta diferenciase 
debió al hecho que en la región norte luego de la soja se cultiva en mayor cantidad de maíz 
que trigo atendiendo a que la temperatura en el periodo invernal no es muyfavorable para el 
cultivo del trigo. En el sur ocurre lo contrario, donde la temperatura no es muy favorable a la 
producción de maiíz de segunda. Con estos ajustes se buscaba representar lo más cercano 
posible la realidad de la agricultura en esos lugares. 

Los datos de costos e ingresos correspondieron a la zafra 2016/2017 para soja y, 
2017 para maiíz y trigo. El precio de los productos agrícolas fue obtenido de consultas con 
técnicos de cada zona, y delos silos de acopio. Se utilizó el precio más frecuentemente 
citado. Los rendimientos fueron los más frecuentes de acuerdo con el nivel tecnológico 
aplicado según los agricultores consultados. Los costos directos se dividieron en insumos 
técnicos e insumos físicos expresados en dólares americanos. Los ingresos igualmente se 
midieron en la misma moneda. 

El análisis económico de la finca se midió a través de los siguientes indicadores 
(Ghida, 2009). 

(i) Margen bruto: fue calculado como la diferencia entre ingreso bruto y costo directo. 
La ecuación usada fue la que se presenta más abajo (1 y2), donde MB es margen bruto, P 
es precio del bien, Q es cantidad del bien, IB es igual al ingreso bruto y CD es costo directo. 
Una vez obtenido el margen bruto, se calculó la relación margen bruto y costo directo 
comúnmente usada como un indicador de rentabilidad entre los agricultores, aunque hay 
que tener presente su limitación ante el hecho que el ratio resultante no incluye los costos 
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indirectos o de estructura. También se relacionó el margen bruto con la superficie total 


cultivada. 
MB=(P*Q)-CD......... (1) 
MP=IB-CD ......<.. (2) 


(ii) Resultado operativo: Fue la diferencia entre margen bruto y costo indirecto o de 
estructura según la ecuación (3) de más abajo, donde RO es resultado operativo, MB es 
margen bruto y Cl es costo indirecto o de estructura. EL costo indirecto incluyó, (i) impuesto 
inmobiliario, (ii) un contador, y (iii) los intereses pagados por el crédito correspondiente 
al capital operativo. El primero se calculó como el 1% del valor fiscal de tierra agrícola, 
(MINISTERIO DE HACIENDA, 2016), del principal Distrito de cada región. Este monto fue 
convertido a dólares americanos. La tasa de cambio utilizada fue de 6i/USD 5.619 igual 
al promedio de la zafra 2016/2017 (BCP, 2020). Los intereses pagados son generalmente 
costos directos, sin embargo, durante el proceso de estructuración de loscostos de la finca 
modal, se tropezó con la realidad de que varios insumos fueron utilizados en más de un 
cultivo, creando la consiguiente dificultad de asignar el costo financiero (el interés pagado) 
de dicho insumo a cada cultivo en forma específica. El capital operativo en su totalidad se 
obtuvo vía crédito bancario y fue igual al costo directo de cada cultivo Las condiciones de 
crédito del capital operativo fueron de seis meses de plazo para cada cultivo, y una tasa de 
8,19% anual pagadero en una sola cuota. Esta tasa fue el valor medio de las tasas activas 
para préstamos en dólares durante la zafra setiembre 2016/agosto 2017 según el Banco 
Central del Paraguay (Banco Central del Paraguay, [BCP], 2020). En el caso de la soja, el 
crédito fue desembolsado en agosto y cancelado en abril sumando 273 días, en tanto que el 
maíz fue de enero a julio con un total de 212 días. En el caso específico del trigo se estimó 
desde a julio (92 días). El servicio de contabilidad fue a razón de (i/mes 500.000 (costo del 
servicio másfrecuente, para fincas hasta 500 ha), igual a Gi/afio 6.500.000, incluyendo un 
mes de aguinaldo, que expresado en dólares totalizó USD/afio 1.157, monto que fue similar 
en todas las regiones. 

Luego se calculó la relación resultado operativo y costo total, y resultado operativo por 
hectárea sembrada, a efectos de compararlos con los anteriores ratios y poder dimensionar 
el peso de los costos indirectos. Finalmente, la participación de los costos indirectos en el 
costo total fue calculado. En la literatura no se ha encontrado una recomendación sobre la 
relación ideal entre costos directos e indirectos. Los autores concuerdan en que se debe, 
dentro de las posibilidades, asignar la mayor parte de los costos a la categoría de costos 
directos (Martínez Ferrario, 1995), (Vargas Otte, 2000), (Ghida, 2009), (Casari y Gorziglia, 
2014). 

RO=MB-CI............ (3) 

(iii) Ingreso neto: se calculó restando al resultado operativo las amortizaciones 


indirectasde los activos no utilizados directamente en los cultivos, tales como construcciones, 
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maquinarias, equipos, automotores y otros capitales de trabajo. Se los destinó a remunerar 
los factores de producción (tierra, capital y trabajo no remunerado). En el disefio de las 
fincas modales se acordó, no incluir ninguna mejora o construcción, tampoco vehículos 
automotores, entre los activos, ante la dificultad para acordar una dotación representativa. 
En resumidas cuentas, el ingreso neto fue similar al resultado operativo. En la fórmula (4), 
la utilizada, IN es ingreso neto, RO es resultado operativo y A es amortizaciones. 

IN=RO-A......... (4) 

(iv) Ingreso al capital: es lo que quedó de los ingresos netos luego de remunerar la 
manode obra del productor y su familia que no fue incluida como costo directo, tampoco 
como indirecto. En este caso se utilizó como proxy del valor de la mano de obra familiar 
no remunerada el de un salario mínimo anual incluyendo aguinaldo que totalizó USD/afio 
5.073 (BCP, 2020). En la fórmula (5), IC es ingreso al capital, IN es ingreso neto, MOF es 
mano de obra familiar. El ingreso al capital tuvo dos fines, remunerar el interés del capital 
de la finca, calcular la rentabilidad del activo promedio de la finca. 

IC=IN-MOF............. (5) 

(v) Utilidad líquida: En la formula (6), UL es utilidad líquida, IC es ingreso al capital e 
IK es interés del capital. La utilidad líquida fue el remanente luego de remunerar el capital 
de la finca, formado por el capital de trabajo (efectivo en caja) y la tierra. El costo del 
capital de trabajo fue igual a los intereses pagados por el crédito, mientras que para el 
costo de la tierra se tomó como proxy el precio de alquiler. Se identificaron dos formas de 
pago por el alquiler de la tierra. La primera fue de 1 t/ha y, la segunda fue el equivalente 
el 25% de la producción cosechada entregada en el silo. En este trabajo se optó por el 
segundo por ser el más usual en las regiones. La utilidad líquida está destinada a retribuir el 
trabajo gerencial de la finca. A parte de estas retribuciones, más bien de corto plazo, Ghida 
(2018), sefiala que una utilidad líquida positiva es un indicador de que la finca con el nivel 
tecnológico aplicado y el alquiler pagado por la tierra es viable en el largo plazo. 

UL=IC-IK............. (6) 

(vi) Rentabilidad: fue el cociente entre el ingreso al capital y el activo promedio 
expresado en porcentaje. Representó la tasa de interés que rindió el activo utilizado enla 
finca, que al início de las actividades estuvo formado por la tierra, a un precio de USD/ha 
10.000,y el capital operativo. Al final del ciclo el valor de la tierra se mantuvo sin cambios, al 
que se le sumó el saldo en efectivo (resultado operativo). La fórmula (7), la utilizada, R es 
rentabilidad; IC es ingreso al capital; ACT es activo promedio. 

R=(IC/ACTP)*100 ..................... (7) 


Investigación, tecnología e innovación en ciencias agrícolas Capítulo 5 FER 


RESULTADO Y DISCUSIÓN 


Región Norte 

La finca representativa de la región norte tuvo 500 hectáreas de superficie, en 
las cuales se sembró soja en su totalidad, más 400 hectáreas de maiz en enero y 100 
hectáreas de trigo en mayo, totalizando 1.000 hectáreas sembradas. El Distrito de Minga 
Pora (24º52'34”S 54º56'00"0) fue la ubicación de referencia. 

El ingreso bruto total de la finca fue USD 822.000 del cual la soja participó con el 
65%, el maíz con el 29% y el trigo con el 6%. Por su lado, el costo directo total llegó a USD 
487,867, monto al cual la soja aportó el 49%, el maíz el 42% y el trigo el 9%. La diferencia 
entre ambos dio como resultado un margen bruto de USD 334.133. Del margen bruto total 
correspondió a la soja el 88%, al maíz el 10% y al trigo el 2%. La relación margen bruto y 
costo directo fue de 68%. Lo que indicó que, por cada unidad monetaria de costo directo, 
el retorno fue 68 unidades, con un rango que fue desde 122% en el caso de la soja y 14% 
del trigo, con el maíz en 17%. Por su lado la relación entre el margen bruto y la superficie 
sembrada dio un resultado de USD/ha 334. 

Cubierto los costos directos, se procedió a calcular los costos indirectos o de 
estructura, que restado del margen bruto dio lo que se conoce como resultado operativo. El 
valor fiscal de tierra agrícola del Distrito de Minga Pora, para el 2017 fue de Gi/ha 1.160.320 
(MINISTERIO DE HACIENDA, 2016), totalizando 6 580.160.000 para toda la finca. Sobre 
esta base el impuesto inmobiliario en dólares fue igual a USD 1.032. El contador tuvo un costo 
de USD/afio 1.157. El costo financiero de los créditos dio un total de USD 25.691, según el 
siguiente detalle: soja (USD 14.828), maíz (USD 9.980) y trigo (USD 883), totalizando USD 
25.691. En total los costos indirectos sumaron USD 27.880, que descontados del margen 
bruto dio un resultado operativo de USD 306.253. En otros termino, el ingreso de la finca 
modal fue suficiente para cubrir todos los costos tanto directos como indirectos dejando 
un saldo destinado a remunerar factores de producción. La relación entre el resultado 
operativo y el costo total fue de 59%, frente a los 68% cuando solamente se tuvo en cuesta 
los costos directos. Igualmente se calculó el resultado operativo por hectárea cultivada, 
que dio USD/ha 306, que muestra los efectos de la inclusión de los costos indirectos en el 
cálculo. Además, estimar los costos indirectos permitió calcular el costo total, que fue igual 
a USD 515.747, de los cuales el 5% correspondió a costos indirectos y el saldo a costos 
directos. 

El ingreso neto fue similar al resultado operativo y se utilizó para remunerar 
el factor trabajo o la mano de obra familiar que no formó parte de los costos. El saldo 
restante denominado ingreso al capital, utilizado para remunerar el capital de la finca, y la 
rentabilidad del activo promedio. 

El capital de la finca estuvo formado por el capital de trabajo (efectivo en caja) y la 
tierra. EL costo del capital o la remuneración total del capital fue igual a por un lado al costo 
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de oportunidad del factor tierra, que para esta finca alcanzó USD 133.000. Este monto fue 
calculado de la siguiente manera 500 ha x 3,80 t/ha x 280 USDA x 0,25. Por el otro lado, 
se tuvo el costo del capital de trabajo, que fue igual a los intereses pagados por el crédito, 
que sumaron USD 25.691 tal como ya calculó anteriormente. En total el costo del capital a 
remunerar totalizó USD 158.691, que descontado del ingreso al capital dio como resultado 
una utilidad líquida de USD 142.489, destinado a retribuir el trabajo gerencial de la finca. 

En el cálculo de la rentabilidad se procedió primeramente a calcular el activo 
promedio de la finca compuesto por la tierra y el efectivo en caja. Al inicio en caja se tenía 
el crédito para soja (USD 238.747) y la tierra con un valor de USD/ha 10.000 totalizando 
un capital promedio de USD 5.000.000. Al final del ciclo el valor de la tierra se mantuvo sin 
cambios, al que se le sumó el saldo en efectivo (resultado operativo) igual a USD 306.258. 
De este modo el valor del activo promedio fue de USD 5.272.500. La rentabilidad fue de 
5,71% que estuvo por encima de la la tasa de interés efectivo pasiva, en moneda extranjera 
para un certificado de depósito de ahorro fue 4,14% (BCP, 2021). Atendiendo a que la inversión 
estuvo en su mayor parte dolarizada, se comparó la rentabilidad alcanzada con la tasa de 
inflación de Estados Unidos, que en el 2016 fue de 1,26% (WORLD BANK, 2021), de modo 
que el retorno alcanzado mantuvo el valor de la moneda y rindió en términos reales 4,45%. 
Región Campo 9 

La finca representativa de la región Campo 9 tuvo una superficie de 400 hectáreas, 
en las cuales se sembró soja en su totalidad, más 200 hectáreas de maíz en enero y 200 
hectáreas de trigo en mayo, totalizando 800 hectáreas sembradas. El Distrito de J. Eulogio 
Estigarribia sirvió de referencia (25º22'24"S 55º42/15"0). 

El ingreso bruto total de la finca totalizó USD 558.000 en el cual la soja participó con 
el 58%, el maíz con el 23% y el trigo con el 19%. Por su lado, el costo directo total Ilegó 
a USD 387.240, monto al cual la soja aportó el 50%, el maíz el 28% y el trigo el 22%. La 
diferencia entre ambos dio como resultado un margen bruto de USD 170.760. Del margen 
bruto total correspondió a la soja el 76%, al maíz el 10% y al trigo el 14%. La relación 
margen bruto y costo directo fue de 44%. Lo que indica que, por cada unidad monetaria de 
costo directo, el retorno fue 44 unidades, con un rango que fue desde. 66% en el caso de 
la soja y 16% del maiíz, con el trigo en 29%. De igual modo se calculó la relación entre el 
margen bruto y la superficie sembrada, que dio un resultado de USD/ha 218. 

Cubierto los costos directos, se procedió a calcular los costos indirectos o de 
estructura, que restado del margen bruto dio lo que se conoce como resultado operativo. 
El valor fiscal de tierra agrícola del Distrito de J. Eulogio Estigarribia, para el 2017 fue 
de Gi/ha 1.139.600 (MINISTERIO DE HACIENDA, 2016), totalizando Gi 455.840.000 para 
toda la finca. Sobre esta base el impuesto inmobiliario en dólares fue igual a USD 811. El 
contador tuvo un costo de USD/afio 1.157. El costo financiero de los créditos dio un total de 
USD 19.096, según el siguiente detalle: soja (USD 12.049), maiz (USD 5.253) y trigo (USD 
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1.750). En total los costos indirectos sumaron USD 21.064, que descontados del margen 
bruto dio un resultado operativo de USD 149.695. En otros termino, el ingreso de la finca 
modal fue suficiente para cubrir todos los costos tanto directos como indirectos dejando un 
saldo destinado a remunerar factores de producción. 

La relación entre el resultado operativo y el costo total fue de 37%, frente a los 44% 
cuando solamente se tuvo en cuesta los costos directos. Igualmente se calculó el resultado 
operativo por hectárea cultivada, que dio USD/ha 187. 

Además, estimar los costos indirectos permitió calcular el costo total, que fue igual 
a USD 408.305, de los cuales el 5% correspondió a costos indirectos y el saldo a costos 
directos. 

Tal como se explicó en la metodologia el ingreso neto fue igual al resultado operativo, 
y se utilizó para remunerar el factor trabajo o la mano de obra familiar que no formó parte 
de los costos. El saldo restante denominado ingreso al capital, utilizado para remunerar el 
capital de la finca, y la rentabilidad del activo promedio. 

El capital de la finca estuvo formado por el capital de trabajo (efectivo en caja) y la 
tierra. EL costo del capital o la remuneración total del capital fue igual a por un lado al costo 
de oportunidad del factor tierra, que para esta finca alcanzó USD 81.000. Este monto fue 
calculado de la siguiente manera 400 ha x 3,00 t/ha x 270 USDA x 0,25. Por el otro lado, 
se tuvo el costo del capital de trabajo, que fue igual a los intereses pagados por el crédito, 
que sumaron USD 19.096 tal como ya calculó anteriormente. En total el costo del capital a 
remunerar totalizó USD 100.096, que descontado del ingreso al capital dio como resultado 
una utilidad líquida de USD 44.526, destinado a retribuir el trabajo gerencial de la finca. 

En el cálculo de la rentabilidad se procedió primeramente a calcular el activo 
promedio de la finca compuesto por la tierra y el efectivo en caja. Al inicio en caja se tenía 
el crédito para soja (USD 194.725) y la tierra con un valor de USD/ha 10.000 totalizando 
un capital promedio de USD 4.000.000. Al final del ciclo el valor de la tierra se mantuvo sin 
cambios, al que se le sumó el saldo en efectivo (resultado operativo) igual a USD 149.695. 
De este modo el valor del activo promedio fue de USD 4.107.895. La rentabilidad fue de 
3,52% que estuvo por debajo de la tasa de interés efectivo pasiva, en moneda extranjera 
para un certificado de depósito de ahorro fue 4,14% (BCP, 2021). Atendiendo a que la 
inversión estuvo en su mayor parte dolarizada, se comparó la rentabilidad alcanzada con la 
tasa de inflación de Estados Unidos, que en el 2016 fue de 1,26% (WORLD BANK, 2021), 
de modo que el retorno alcanzado mantuvo el valor de la moneda y rindió en términos 
reales 2,26%. 


Región Centro 


La finca representativa de la región centro tuvo una superficie de 300 hectáreas, 
en las cuales se sembró soja en su totalidad durante la primavera, más 150 hectáreas de 
maiz en enero y 150 hectáreas de trigo en mayo, totalizando 600 hectáreas sembradas. El 
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Distrito de Iguazú fue la referencia geográfica (25º28'00”S 55º00'00"0). 

El ingreso bruto total de la finca totalizó USD 490.125 en el cual la soja participó con 
el 67%, el maíz con el 18% y el trigo con el 15%. Por su lado, el costo directo total Ilegó 
a USD 274.227, monto al cual la soja aportó el 53%, el maiz el 27% y el trigo el 20%. La 
diferencia entre ambos dio como resultado un margen bruto de USD 215.898. Del margen 
bruto total correspondió a la soja el 85%, al maíz el 7% y al trigo el 8%. La relación margen 
bruto y costo directo fue de 79%. Lo que indica que, por cada unidad monetaria de costo 
directo, el retorno fue 79 unidades, con un rango que fue desde. 126% en el caso de la soja 
y 20% del maiíz, con el trigo en 30%. De igual modo se calculó la relación entre el margen 
bruto y la superficie sembrada, que dio un resultado de USD/ha 360. 

Cubierto los costos directos, se procedió a calcular los costos indirectos o de 
estructura, que restado del margen bruto dio lo que se conoce como resultado operativo. 
El impuesto inmobiliario fue USD 619, calculado sobre un valor fiscal de tierra agrícola 
del Distrito de Minga Guazú, que para el 2017 fue de Gi/ha 1.160.320 (MINISTERIO DE 
HACIENDA, 2016), totalizando Gi 348.096.000. El contador tuvo un costo de USD/afio 1.157. 
El costo financiero de los créditos dio un total de USD 13.753, según el siguiente detalle: 
soja (USD 9.070), maíz (USD 3.533) y trigo (USD 1.150). En total los costos indirectos 
sumaron USD 15.530, que descontados del margen bruto dio un resultado operativo de 
USD 200.369. En otros termino, el ingreso de la finca modal fue suficiente para cubrir 
todos los costos tanto directos como indirectos dejando un saldo destinado a remunerar 
factores de producción. La relación entre el resultado operativo y el costo total fue de 69%, 
frente a los 79% cuando solamente se tuvo en cuesta los costos directos. Igualmente se 
calculó el resultado operativo por hectárea cultivada, que dio USD/ha 334, que muestra los 
efectos de la inclusión de los costos indirectos en el cálculo. Además, estimar los costos 
indirectos permitió calcular el costo total, que fue igual a USD 289.756, de los cuales el 5% 
correspondió a costos indirectos y el saldo a costos directos. 

Tal como se explicó en la metodologia el ingreso neto fue igual al resultado operativo, 
y se utilizó para remunerar el factor trabajo o la mano de obra familiar que no formó parte 
de los costos. El saldo restante denominado ingreso al capital, utilizado para remunerar el 
capital de la finca, y la rentabilidad del activo promedio. 

El capital de la finca estuvo formado por el capital de trabajo (efectivo en caja) y la 
tierra. EL costo del capital o la remuneración total del capital fue igual a por un lado al costo 
de oportunidad del factor tierra, que para esta finca alcanzó USD 82.650. Este monto fue 
calculado de la siguiente manera 300 ha x 3,80 t/ha x 290 USDA x 0,25. Por el otro lado, 
se tuvo el costo del capital de trabajo, que fue igual a los intereses pagados por el crédito, 
que sumaron USD 13.753 tal como ya calculó anteriormente. En total el costo del capital a 
remunerar totalizó USD 96.403, que descontado del ingreso al capital dio como resultado 
una utilidad líquida de USD 98.892, destinado a retribuir el trabajo gerencial de la finca. 

En el cálculo de la rentabilidad se procedió primeramente a calcular el activo 
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promedio de la finca compuesto por la tierra y el efectivo en caja. Al inicio en caja se tenía 
el crédito para soja (USD 146.044) y la tierra con un valor de USD/ha 10.000 totalizando 
un capital promedio de USD 3.000.000. Al final del ciclo el valor de la tierra se mantuvo sin 
cambios, al que se le sumó el saldo en efectivo (resultado operativo) igual a USD 200.369. 
De este modo el valor del activo promedio fue de USD 3.173.206. La rentabilidad fue de 
6,15% que estuvo por encima de la tasa de interés efectivo pasiva, en moneda extranjera 
para un certificado de depósito de ahorro fue 4,14% (BCP, 2021). Atendiendo a que la 
inversión estuvo en su mayor parte dolarizada, se comparó la rentabilidad alcanzada con la 
tasa de inflación de Estados Unidos, que en el 2016 fue de 1,26% (WORLD BANK, 2021), 
de modo que el retorno alcanzado mantuvo el valor de la moneda y rindió en términos 
reales 4,89%. 


Región Sur 

La finca representativa de la región sur tuvo 50 hectáreas de superficie, en las cuales 
se sembró soja en su totalidad, más 10 hectáreas de maiíz en enero y 40 hectáreas de 
trigo en mayo, totalizando 100 hectáreas sembradas. El Distrito de San Rafael del Paraná 
(26º39'S 54º56'0) actuó de referencia geográfica. 

El ingreso bruto total de la finca totalizó USD 62.850 en el cual la soja participó con 
el 64%, el maíz con el 10% y el trigo con el 27%. Por su lado, el costo directo total Ilegó 
a USD 43.211, monto al cual la soja aportó el 54%, el maíz el 12% y el trigo el 34%. La 
diferencia entre ambos dio como resultado un margen bruto de USD 19.639. Del margen 
bruto total correspondió a la soja el 85%, al maíz el 4% y al trigo el 11%. La relación margen 
bruto y costo directo fue de 45%. Lo que indica que, por cada unidad monetaria de costo 
directo, el retorno fue 45 unidades, con un rango que fue desde 71% en el caso de la soja 
y 15% del maiíz, con el trigo en 16%. De igual modo se calculó la relación entre el margen 
bruto y la superficie sembrada, que dio un resultado de USD/ha 196. 

Luego se calcularon los costos indirectos o de estructura, que restado del margen 
bruto dio lo que se conoce como resultado operativo. El impuesto inmobiliario fue USD 
140, que se calculó como el 1% del valor fiscal de tierra agrícola del Distrito de San Rafael 
del Paraná, que para el 2017 fue de Gi/ha 1.574.720 (MINISTERIO DE HACIENDA, 2016), 
totalizando À 78.735.00. El contador tuvo un costo de USD/ahio 1.157. El costo financiero 
de los créditos dio un total de USD 7.013, según el siguiente detalle: soja (USD 1.453), maíz 
(USD 2.495) y trigo (USD 3.065). En total los costos indirectos sumaron USD 8.310, que 
descontados del margen bruto dio un resultado operativo de USD 11.329. En otros termino, 
el ingreso de la finca modal fue suficiente para cubrir todos los costos tanto directos como 
indirectos dejando un saldo destinado a remunerar factores de producción. La relación 
entre el resultado operativo y el costo total fue de 22%, frente a los 45% cuando solamente 
se tuvo en cuesta los costos directos. Igualmente se calculó el resultado operativo por 
hectárea cultivada, que dio USD/ha 113, que muestra los efectos de la inclusión de los 
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costos indirectos en el cálculo. Además, estimar los costos indirectos permitió calcular 
el costo total, que fue igual a USD 51.521, de los cuales el 16% correspondió a costos 
indirectos y el saldo a costos directos. 

Tal como se explicó en la metodologia el ingreso neto fue igual al resultado operativo, 
y se utilizó para remunerar el factor trabajo o la mano de obra familiar que no formó parte 
de los costos. El saldo restante denominado ingreso al capital, utilizado para remunerar el 
capital de la finca, y la rentabilidad del activo promedio. 

El capital de la finca estuvo formado por el capital de trabajo (efectivo en caja) y 
la tierra. EL costo del capital o la remuneración total del capital fue igual a por un lado al 
costo de oportunidad del factor tierra, que para esta finca alcanzó USD 10.013. Este monto 
fue calculado de la siguiente manera 50 ha x 3,00 t/ha x 267 USDA x 0,25. Por el otro 
lado, se tuvo el costo del capital de trabajo, que fue igual a los intereses pagados por el 
crédito, sumaron USD 7.013 tal como ya calculó anteriormente. En total el costo del capital 
a remunerar totalizó USD 17.026, que descontado del ingreso al capital dio como resultado 
una utilidad líquida negativa de USD 10.770. Una utilidad líquida negativa indica que la 
empresa agrícola, en esa región, con el actual nivel tecnológico y uso de suelo no es viable 
en el largo plazo. 

En el cálculo de la rentabilidad se procedió primeramente a calcular el activo 
promedio de la finca compuesto por la tierra y el efectivo en caja. Al inicio en caja se tenía 
el crédito para soja (USD 23.392) y la tierra con un valor de USD/ha 10.000 totalizando 
un capital promedio de USD 500.000. Al final del ciclo el valor de la tierra se mantuvo sin 
cambios, al que se le sumó el saldo en efectivo (resultado operativo) igual a USD 11.329. 
De este modo el valor del activo promedio fue de USD 517.360. La rentabilidad fue de 
1,21% que fue menor a la tasa de interés efectivo pasiva, en moneda extranjera para un 
certificado de depósito de ahorro que fue 4,14% (BCP, 2021). Atendiendo a que la inversión 
estuvo en su mayor parte dolarizada, se comparó la rentabilidad alcanzada con la tasa de 
inflación de Estados Unidos, que en el 2016 fue de 1,26% (WORLD BANK, 2021), el retorno 


alcanzado fue levemente inferior. 


CONCLUSIÓN Y RECOMENDACIONES 


Con respecto a objetivo, este trabajo realizó la evaluación y concluye que las fincas 
cubren sus costos directos e indirectos. En la remuneración de los factores de producción, 
todas remuneraron mano de obra familiar, pero solamente tres hicieron lo propio con capital 
y tierra. En lo referente a la rentabilidad, dos fincas, las de mayor extensión, alcanzaron 
valores superiores al costo de oportunidad. Respecto a la comparación con la inflación de 
EEUU, atendiendo a que los costos e ingresos de los cultivos se realizan en la moneda de 
dicho país, todas las fincas, excepto la de la región sur, tuvieron rentabilidad superior a la 
inflación de los EEUU. La recomendación apunta a realizar el trabajo anualmente, e incluir 
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un análisis financiero. 





Indicadores Región norte Región centro Región Campo Región Sur 
Ingreso bruto total 822,000 490,125 558,000 62,850 
Ingreso bruto Soja 532,000 330,600 324,000 40,050 
Ingreso bruto Maiz 242,000 88,275 126,000 6,000 
Ingreso bruto Trigo 48,000 71,250 108,000 16,800 
Costo Directo 487,867 274,227 387,240 43,21 
CD Soja 238,747 146,044 194,725 23,392 
CD Maiíz 206,920 73,260 108,908 5,173 
CD Trigo 42,200 54,923 83,608 14,646 
Margen Bruto 334,133 215,898 170,760 19,639 
MB Soja 293,253 184,556 129,275 16,658 
MB Maíz 35,080 15,015 17,092 827 
MB Trigo 5,800 16,327 24,392 2,154 
Gastos de estructura y administrativos 27,880 15,530 21,064 8,310 
Otros Egresos Indirecto 

Resultado Operativo 306,253 200,369 149,695 11,329 
Amortización indirecta de mejoras 0 0 0 0 
Amortización indirecta de maquinarias y equipos 0 0 0 0 
Amortización de automotores y otros capitales 0 0 0 0 
Ingreso neto 306,253 200,369 149,695 11,329 
Remuneración mano de obra del productora 5,073 5,073 5,073 5,073 
Ingreso al capital 301,180 195,296 144,622 6,255 
Rentabilidad 5.71% 6.15% 3.52% 1.21% 
Utilidad líquida 142,489 98,892 44,526 -10,770 


Fuente: Projecto 12-inv-122 "La rentabilidad de los rubros agrícolas conectados a mercados 
internaciononales en los dos estratos de la agricultura paraguaya: empresarial y familiar. 


Tabla 1. Indicadores del análisis económico (USD). 
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tamento del Cauca en la cadena avícola nacional se 
caracteriza por pequefios sistemas de producción que fortalecen la economía familiar 
especialmente en los municipios que hacen parte de la meseta de Popayán (Carvajal et al, 
2017); estos sistemas abastecen la demanda local que contribuyen a mejorar la calidad de 
vida y la seguridad alimentaria de las comunidades rurales de la región. 

En la alimentación de pollos de engorde, se han identificado especies vegetales útiles 
que están subvaloradas y no son aprovechados adecuadamente (Carvajal et al., 2017), 
tal como, Tithonia diversifolia, conocida en la región como botón de oro, especie versátil 
que se adapta desde O hasta 2400 msnm, soporta diversas temperaturas ambientales, 
y condiciones adversas como la baja fertilidad de los suelos y pobre precipitación anual 
(Argúello et al, 2019; Gutiérrez et al., 2019), no obstante, el mayor potencial productivo de 
la especie se expresa en clima medio. T. diversifolia contiene niveles de proteína cruda del 
25,07%, energia bruta de 4223+4,04 kcal GE/kg (Buragohain, 2016) y alta digestibilidad 
de la materia (Rodríguez et al, 2018). Adicionalmente, en las especies forrajeras presenta 
resultados significativos en producción de biomasa aérea, características que pueden 
potencializar su uso con el objeto de disminuir los costos de producción en la alimentación 
de pollos de engorde, y en general, para especies de interés zootécnico (Zabala, 2021). 

El objetivo de la investigación fue determinar el efecto de tres niveles de inclusión 
de follaje de Tithonia diversifolia (0%, 17%, 21% y 25% en base seca) en parámetros 
productivos de pollos de engorde para la etapa de finalización. 


MATERIALES Y MÉTODOS 


La investigación se realizó en la finca Bella Vista, vereda La Estación del municipio de 
Morales (Cauca, Colombia), localizada a 1.635 msnm., 2º45'37” de latitud norte 76º38'02” 
de longitud oeste, con temperaturas entre 23 - 28ºC y humedad relativa promedio del 80%. 


Animales de estudio 


Se utilizaron 72 pollos de engorde de 27 días de edad con plan de vacunación 
completo, de la línea COOB 500, alimentados con concentrado comercial y con un peso 
promedio de 1950+20 g. Los animales se sometieron a un periodo de acostumbramiento 
de 3 días antes de suministrar la dieta experimental y se mantuvieron 15 días en 
experimentación para una edad de sacrificio de 45 días. 

El manejo de los animales se caracterizó por emplear una cama de aserrín, previas 
medidas de salubridad del galpón, desinfección y flameado de pisos y paredes, régimen de 
luz de 24 horas, temperatura bajo condiciones naturales de la zona, regulando el confort de 
los animales mediante el uso de cortinas, densidad de 11 pollos por m?, alimentación en dos 
raciones (6:00 am y 1:00 pm) y suministro de agua a voluntad en comederos y bebederos 


plásticos convencionales. 
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Obtención de la harina de hojas de T. diversifolia 


Las hojas se cosecharon de plantas en etapa de floración, sin manejo agronómico, 
mediante corte manual en época de bajas precipitaciones; posteriormente, se deshidrataron 
durante tres días en un secador parabólico y se molieron hasta obtener una harina con 
el 25% de materia seca, respecto al contenido de humedad inicial. La evaluación de la 
composición química de la harina de hoja, fue realizada por el laboratorio de nutrición de la 
Universidad de Narifio (tabla 1), en cuanto a los valores del contenido de Lisina y Metionina 
de T. diversifolia de 1,08% y 0,23% utilizados respectivamente, se obtienen de los estudios 
realizados por Grashorn (2017). 


Parâmetro Harina de hoja T. diversifolia Método de análisis en laboratorio 
PC (%) 29,59 Volumétrico Kjeldahl 

EM Mcal/kg 2,54 Calorimetría 
Ca (%) 2,47 Espectrofotometria UV-VIS 
P(%) 0,36 Espectrofotometríia UV-VIS 
FC (%) 6,49 Gravimétrico 


Tabla 1. Bromatológico de harina de hoja de T. diversifolia en base seca 


Fuente: análisis bromatológico de la harina de hojas de T. diversifolia, realizado en el laboratorio de la 
Universidad de Nariho (Pasto-Narifio). PC: proteína cruda; EM: energia metabolizable; Ca: calcio; P: 
fósforo; FC: fibra cruda. 


Elaboración del concentrado 


Para balancear las dietas iso-proteicas e iso-energéticas, se tuvieron en cuenta los 
requerimientos nutricionales de las tablas brasilefias para cerdos y aves de Rostagno et al. 
(2017), empleando maiíz blanco, sorgo, harina de pescado, torta de soya, aceite vegetal, 
pre mezcla de vitaminas y minerales, sal marina, melaza, harina de hojas de T. diversifolia, 
carbonato de calcio, fosfato bicálcico, metionina y monoclorhidrato de lisina; el concentrado 
en forma de pellets se elaboró utilizando un método artesanal con molino. 

Variables evaluadas 


Consumo de alimento (C): se calculó mediante diferencia entre el peso del alimento 
ofrecido y el peso del alimento rechazado, empleando una balanza electrónica digital 
gramera (10kg/19). La cantidad de alimento rechazado fue registrado antes de suministrar 
la ración, según los residuos encontrados en el comedero realizando el monitoreo dos 
veces al día. 

Ganancia de peso (GP): resultado de la diferencia entre el peso final y el inicial 
(después de los 3 días de acostumbramiento), de cada uno de los animales por tratamiento. 
Se empleó una balanza mecánica tipo reloj de 100kg. 
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Conversión alimenticia (CA). Se calculó dividiendo el consumo de alimento de los 
15 días de evaluación entre la ganancia de peso. 


CA=C/GP 


Disefio experimental y análisis estadístico 


Se empleó un disefio experimental completamente al azar, donde el factor tratamiento 
está determinado por cuatro dietas para la alimentación en etapa de finalización, el control 
y tres niveles de inclusión de harina de hojas de T. diversifolia: 17%, 21 % y 25%, definidos 
como tratamientos TDO, TD17, TD21 y TD25 respectivamente, con tres repeticiones 
distribuidas al azar y seis animales por repetición. La información se tabuló en una matriz 
de Excel para cada variable evaluada. El análisis estadístico de la información se realizó 
en el software SPSS v 25, mediante un análisis de varianza (ANOVA) con un margen de 
significancia del 5% (p<0,05) y prueba de promedios de Tukey para identificar el mejor 


tratamiento. 


RESULTADOS Y DISCUSIÓN 


El análisis de varianza arrojó diferencias estadísticas (p<0,05) para las variables 
ganancia de peso y conversión alimenticia. Los resultados muestran que en las dietas 
elaboradas con los diferentes niveles de inclusión (17, 21 y 25%), no se afectó el consumo 
de alimento; resultados similares obtuvieron Buragohain y Kalita (2015), quienes evaluaron 
niveles de inclusión del 2, 4, 6 y 8% de botón de oro; a su vez, Molina (2020) sefiala que 
probablemente, el consumo no se afecte porque las aves poseen pocas papilas gustativas; 
por su parte, Valdivié et al., (2008), indica que no hay diferencias en palatabilidad y 
aceptación por parte de las aves, cuando las dietas se ofertan, en forma de harina. 


Ganancia de peso 


En la Figura 1 se puede observar que existieron diferencias significativas entre 
los tratamientos (p<0,05). En este sentido, la prueba de promedios de Tukey sefiala que 
el tratamiento con mejores resultados fue T2 (21%) con una ganancia de peso de 1222 
gr, seguido del T3 (25%) con pesos de 1111 gr; mientras que, el testigo (TO) y T1 (17%) 
presentaron un comportamiento similar con la menor ganancia, y con pesos de 914 gr y 978 
gr respectivamente. Los resultados con buen desempefio productivo (T2 y T3) se pueden 
atribuir a que la dieta con botón de oro, presentó un contenido adecuado de energia y 
proteína (Gutiérrez y Hurtado, 2019). González et al, (2014) sefialan que el uso de plantas 
forrajeras en la dieta de diferentes especies animales es una de las opciones para la 
producción eficiente y rentable en sistemas agropecuarios; ante estos condicionantes, se 
buscan fuentes proteicas económicas, evaluando follajes de especies arbóreas que, por 
su potencialidad, pueden garantizar dietas nutricionalmente adecuadas para animales en 
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producción. 
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Figura 1. Ganancia de peso de pollos de engorde empleando diferentes niveles de inclusión de T. 
diversifolia en la dieta 


Cabe anotar que la fibra encontrada en el follaje de botón de oro (6,49% tabla 1), 
no fue un limitante para el balance y digestibilidad de las dietas experimentales, y sumado 
al contenido de proteína en la harina de botón de oro (Rodríguez et al., 2018), resultaría 
factible encontrar efectos positivos con las mayores inclusiones evaluadas. 

Es importante mencionar, que el nivel de proteína de harina de botón de oro con 
que se trabajó fue del 29,59% en promedio (tabla 1), resultado que podría aumentar la 
disponibilidad para el aprovechamiento de los pollos (Vega et al., 2019). En este sentido, 
Makinde & Egbekun (2016) y Mahmoud et al. (2017) con 22% de PB observaron mejor 
ganancia de peso, conversión alimenticia, digestibilidad de proteína bruta y rendimiento 
de la canal; por su parte, Guilcapi (2013) en la etapa de crecimiento con 20% de PB, logró 
mejores parámetros productivos en cuanto a ganancia de peso y conversión alimenticia y 
en la etapa de engorde con 18% de PB, consiguieron los mejores promedios productivos 
en peso final, ganancia de peso y conversión alimenticia. 


Conversión alimenticia 


Entre los tratamientos, se presentaron diferencias estadísticas para la conversión 
alimenticia (p<0,05). La prueba de promedios de Tukey mostró que el tratamiento con 
mejor conversión fue T2 (21%) con valor de 1,6, seguido de T3 (25%) con factor de 
conversión de 1,8 (Figura 2), mientras que, el control presentó el mayor factor con 2,2. 
Resultados similares obtuvo Connolly (2017) quien reporta conversión alimenticia de 2.12 
con concentrado comercial y al incluir dietas no convencionales presenta valores de 1,64. 
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Figura 2. Conversión alimenticia de pollos de engorde alimentados con diferentes niveles de inclusión 
de T. diversifolia. 


Fenavi (2020) reporta valores de conversión alimenticia de 1,73 con edad al sacrifico 
de 39 días, comparado a los resultados de la presente investigación, a excepción de TO, 
los valores son más eficientes utilizando diferentes niveles de inclusión de harina de botón 
de oro, que al presentar contenidos de proteína entre 20,37 y 23,37% (0 29% como en esta 
investigación), y carbohidratos solubles totales entre 9,65 y 12,92% (Londofio et al., 2019), 
se contempla como una excelente alternativa para la alimentación de pollos de engorde. Al 
respecto, Holguín et al, (2019) sefiala que, si bien la nutrición animal representa en promedio 
un 70% de los costos de producción, los buenos resultados en las producciones avícolas 
se equilibran con una buena genética del animal, manejo adecuado de enfermedades y el 
desarrollo de sistemas para brindar confort ambiental. 

Savón et al.,(2017) por su parte, mencionan que la digestibilidad in vitro al utilizar 
botón de oro con altos niveles de proteína como en la presente evaluación (29%), oscila 
entre 75,28 y 78,59%, lo que lo convierte en una alternativa de forraje para manejo de 
materia prima no convencional en la alimentación de aves y que al utilizar un nivel inclusión 
del 21% de harina de hojas de botón de oro en la alimentación en etapa de finalización, la 
digestibilidad es alta y hay un buen aprovechamiento de nutrientes. 


Análisis económico 


En la tabla 2, se observa que los costos de producción para el manejo de los niveles 
de inclusión de harina de botón de oro son mayores con respecto al TO; sin embargo, es 
importante resaltar que los insumos utilizados para la elaboración del concentrado en los 
tratamientos T2 y T3, es más económico, probablemente, el uso de botón de oro, minimiza 
los costos. 
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Tratamiento 


Descripción 
To T1 T2 T3 
M. obra 70.000 75.000 78.000 81.000 
Insumos 273.464 274.994 272.197 270.718 
Otros 3.375 3.375 3.375 3.375 
Total 346.839 353.368.5 353.571 355.093 


Tabla 2. Costos de producción de los diferentes tratamientos 


El valor del gramo de pollo en el municipio de Morales tuvo un costo de $7, en este 
sentido, T2 fue el tratamiento que mayor ganancia de peso obtuvo (tabla 3), y generó 
mejores ingresos al momento de vender los pollos en canal, presentando una diferencia de 
aproximadamente $200 frente a los tratamientos TO y T1. 


Tratamiento Peso q. Valor gramo $ Valor Total $ 
TO 2.865 7 20.055 
Ti 2.911 7 20.377 
T2 3.233 7 22.631 
T3 3.050 7 21.350 


Tabla 3. Costo en venta promedio de la carne de pollo. 


En la tabla 4 se observa la diferencia entre costos totales, beneficio bruto y utilidad 
entre T2 y los tratamientos TO y T1; se revelan mayores ganancias al utilizar el 21% de 
nivel de inclusión de harina de hoja de botón de oro en la dieta, y el promedio de peso de 
los pollos a los 45 días de edad fue mayor frente a los demás tratamientos y se encuentran 
diferencias económicas que conllevan a tener un mejor margen de rentabilidad del 15,2%, 
por lo cual se minimizan costos y mejoran las utilidades en una producción de manejo 
de pollos de engorde en etapa de finalización. La relación Beneficio-Costo (B/C) está 
presentada por la relación Ingresos-Egresos, en donde ambos deben ser calculados de 
acuerdo con el flujo de caja (Flórez et al., 2018). 


Costos 


Tratamiento totales $ Beneficio bruto $ Utilidad $ 
TO 346.839 360.990 14.151 
T1 353.369 366.799 13.430 
T2 358.571 407.399 58.828 
T3 355.092 384.300 29.208 


Tabla 4. Utilidad obtenida por tratamiento 
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CONCLUSIONES 


El remplazo parcial de harina de hoja de botón de oro que se realizó destacó que 
a un 21% de nivel de inclusión, como materia prima no convencional, es una alternativa 
de alimentación. Por lo tanto, en sistemas de producción campesina, es factible elaborar 
dietas no convencionales para alimentación de pollos en fase de finalización. 

La elaboración del concentrado con un nivel de inclusión de 21% de harina de 
hojas de botón de oro, disminuyen los costos en los insumos y genera mayores ingresos, 
comparada con los concentrados con 25%, 17% y 0% de inclusión de harina de hojas de 
botón de oro. 
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ABSTRACT: Developing studies for water 
treatment that seek to adapt and balance the 
preservation of the environment is fundamental 
for the present time, which is why studies with 
organic coagulants have become so important. 
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TREATMENT 


Aluminum sulfate is widely used, but high 
concentrations can cause problems to human 
health. The tests for water treatment were carried 
out simulating the processes that occur in a water 
treatment plant, such as coagulation, flocculation, 
sedimentation, and filtration. The parameters 
pH, electrical conductivity, and turbidity were 
analyzed. The coagulants Tannin, Aluminum 
Sulfate, Moringa Oleifera, Ferric Chloride, Bufloc 
5122 and, Bufloc 5122 accompanied by the 
auxiliary Bufloc 5158 were used. The results 
showed that the coagulant extracted from the 
seed of Moringa oleifera stood out after direct 
filtration, showing results similar to ferric chloride, 
both showing results higher than 99% of turbidity 
removal. For pH, the results showed no variation, 
as for conductivity, Moringa oleifera presented 
values above the average due to its saline 
composition. It is concluded, therefore, that there 
is a possibility that an organic coagulant presents 
efficiency equivalent to an inorganic coagulant 
in water treatment, bringing benefits to the 
environment as it is a less polluting component. 
KEYWORDS: Coagulants, filtration, water 
treatment. 


AVALIAÇÃO DA EFICIÊNCIA DE 
DIFERENTES AGENTES COAGULANTES 
ASSOCIADO A SISTEMA DE FILTRAÇÃO 

DIRETA NO TRATAMENTO DE ÁGUA 


RESUMO: Desenvolver estudos para tratamento 
de água que procurem adequar e equilibrar a 
preservação do meio ambiente é fundamental 
para os tempos atuais, por isso estudos 
com coagulantes orgânicos tornam-se tão 
importantes. O sulfato de alumínio é bastante 
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empregado, porém elevadas concentrações podem causar problemas à saúde humana. Os 
ensaios para tratamento de água foram realizados simulando os processos que ocorrem numa 
estação de tratamento de água, como a coagulação, floculação, sedimentação e filtração e 
foram analisados os parâmetros pH, condutividade elétrica e turbidez. Foram empregados 
os coagulantes Tanino, o Sulfato de Alumínio, a Moringa Oleifera, o Cloreto Férrico, o Bufloc 
5122 e o Bufloc 5122 acompanhado do auxiliar Bufloc 5158. Os resultados mostraram que 
o coagulante extraído da semente de Moringa oleifera se destacou após a filtração direta, 
apresentando resultados semelhantes ao cloreto férrico, ambos apresentando resultados 
superiores a 99% de remoção de turbidez. Para o pH, os resultados não demonstraram 
variação e para a condutividade a Moringa oleifera apresentou valores acima da média 
devido sua composição salina. Conclui-se, portanto, que há possibilidade de um coagulante 
orgânico apresentar eficiência equivalente a um coagulante inorgânico no tratamento de 
água, trazendo benefícios ao meio ambiente por ser um componente menos poluente. 
PALAVRAS-CHAVE: Coagulantes, filtração, tratamento de água. 


11 INTRODUCTION 


Brazil has 12% of the world's availability of freshwater. Even so, 80% of Brazilian 
water resources are found in the Amazon region, where the smallest part of the population 
lives and the water demand is lower (AGÊNCIA NACIONAL DE ÁGUAS, 2018), such 
outcome is a consequence of the great social and regional inequality present in the region. 
country (MACHADO, 2003). Therefore, developing new technologies and research in the 
area of water treatment that seeks to facilitate access to treated water is essential for the 
development of a disease-free society with quality of life. 

According to Seckler (2017), the water treatment process can be considered as a 
set of water manipulations in its most different presentations, so that it can be classified as 
suitable for public supply. The processes are divided into coagulation (the stage in which 
suspended particles are agglomerated with the aid of a coagulant); flocculation (agglutination 
of smaller flakes of dirt into larger flakes); decantation (deposition of flocs by gravitational 
action) and filtration (the stage in which impurities are retained in a filter). 

For Libânio (2010) filtration consists of a process capable of correcting failures that 
occurred in coagulation, flocculation, and sedimentation, ergo, filtration is responsible for 
retaining the colloidal particles of the water, thus improving the treatment and guaranteeing 
the quality of the water treated. 

Currently, the most used inorganic coagulant for water treatment is aluminum sulfate 
(AI(SO,),), but high concentrations of aluminum in the environment can be a cause of 
problems to human health (ROSALINO, 2011). Therefore, the use of organic coagulants 
based on seeds and tree bark has been a viable solution, as they have low toxicity and 
greater sustainability (SANTOS, 2011). 

With this in mind, this study aims to approach water treatment in the sustainable 
concept of using organic coagulants as a way of replacing or aiding the use of inorganic 
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coagulants in the process involved. 


21 MATERIAL AND METHODS 


The tests were carried out at the UTFPR Water Resources Laboratory — Londrina 
campus. 12 L of synthetic water based on a kaolinite solution (Al,Si,O,(OH),) were prepared 
to simulate turbidity levels close to 100 uT using the methodology of Mendes (1989). The 
research included six treatments using different coagulants, namely Tannin (T1), Aluminum 
Sulfate (T2), Moringa Oleifera (T3), Ferric Chloride (T4), Bufloc 5122 (T5), and Bufloc 5122 
accompanied by auxiliary Bufloc 5158 (T6). A preliminary assay was carried out in which 
several concentrations of coagulants were randomly tested and it was defined that the most 
efficient concentration for T1, T2, and T4 was 6 mL L', for T3 and T5 a concentration of 2 
mL L“ was used. 1 and 4 mL L” respectively. For the T6 treatment, the concentration of 3 
mL L was adopted. 

To prepare the coagulant solutions, the following procedures were performed: for 
the coagulant extraction solution based on Moringa Oleifera seed, 10 g of the seed was 
crushed and diluted with 1 M NaClin 1 L of distilled water. The solution was processed in 
a blender, and then the mixture was strained through a cloth sieve to remove the coarser 
organic matter, leaving only the coagulant liquid. Aluminum Sulfate, Ferric Chloride, and 
Bufloc 5122 coagulants were prepared by dissolving 1 g of the coagulantin 1 L of distilled 
water. Tannin and Bufloc 5158 were prepared by diluting 1 mL of each coagulant in 1 L of 
distilled water, and for Bufloc 5158 it was necessary to heat the water so that its complete 
dissolution was possible. 

In the Jar-test equipment, the experiment was executed in duplicate and after adding 
the coagulants in their respective concentrations, the equipment was activated and the 
coagulation process was carried out for 3 minutes at 150 rpm, then reducing the rotation to 
15 rpm. for 10 minutes in order to start the flocculation process, then turn off the device to 
characterize the sedimentation process (adapted from THEODORO (2012)). 

For the construction of the filter, pre-tests were carried out where 3 filters were tested 
aiming to choose the best one for the development of this work and so the chosen filter was 
produced with sand with a granulometry from O to 0.425 mm (G1) in which it was filled 15 
centimeters of the filter composing the top layer and a base layer of nonwoven geotextile. 

The analyzes during the sedimentation process were performed at an interval of 10 
in 10 minutes, with the first collection being performed with 3 minutes of sedimentation and 
the last with 33 minutes. After collections, the pre-treated water from the jars was directed 
to the filters and after the passage of the water, the samples were collected for the last 
analysis. 

The parameters analyzed were pH, electrical conductivity, and turbidity, all of which 
were analyzed by the methodology of APHA (2012). 
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31 RESULTS AND DISCUSSION 


Initially, the parameters observed for water in its raw state were pH 8.6; electrical 
conductivity 89.55 yS em” and turbidity 107.5 uT, and Total Dissolved Solids 43 ppm. 
Figure 1 shows the behavior of pH during the test. 
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Note: T1: tannin; T2: aluminum sulfate; T3: Moringa oleifera; T4: ferric chloride; T5: Bufloc 5122; T6: 
Bufloc 5122 and Bufloc 5158. 


Figure 1 — Graph of pH behavior during the test. 
Source: Own Authorship (2021). 


When analyzing Figure 1, it is observed that the pH values fell as time progressed, 
remaining after filtration with values lower than the initial gross value of 8.6, a behavior 
corresponding to that presented by Borges and Pereira (2020) in their studies. Notably, of 
all treatments, treatment T3 had the lowest pH values, resulting in a pH value of 6.2 at the 
end of the test. It is worth mentioning that at the end of the experiment, during the passage 
through the filter medium, the T6 treatment showed the highest pH compared to the other 
treatments, showing values of 6.5. Despite this, the results showed that all treatments did 
not show great pH variation between them at the end of the test. 

Figure 2 shows the values obtained from the variation of electrical conductivity during 
the experiment. Observing the results, it can be seen that the T3 treatment presented, 
taking into account the gross value of 89.55 yuS cm”, high levels of electrical conductivity 
throughout the experiment, which was also observed by Borges and Pereira (2020). 

Itis noted that at 23 minutes the treatment reached its maximum peak when reaching 
values of 269.8 yS cm and this uneven change to the other treatments is due to the fact 
that the coagulant based on Moringa oleifera seed extraction has the presence of NaCl in 
its composition, and once the salt releases Na* ions into the solution, there is a considerable 
increase in electrical conductivity, as also pointed out by the study by Vizibelli et al. (2019). 
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Note: T1: tannin; T2: aluminum sulfate; T3: Moringa oleifera; T4: ferric chloride; T5: Bufloc 5122; T6: 
Bufloc 5122 and Bufloc 5158. 


Figure 2 — Graph of the electrical conductivity of values obtained during the test. 
Source: Own Authorship (2021). 


The electrical conductivity of the other treatments shown in Fig. 2 showed no 
significant differences in the gross water value, however, it is noticeable that in the T1 
treatment there was an increase in the value obtained in the first analysis at 3 minutes of 
70.3 yS cm”, with a value of 87 yS cm” at the end of filtration. Current Brazilian legislation 
does not present limit values for electrical conductivity. Figure 3 represents the turbidity 
removal efficiency results during the experiment. 
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Note: T1: tannin; T2: aluminum sulfate; T3: Moringa oleifera; T4: ferric chloride; T5: Bufloc 5122; T6: 
Bufloc 5122 and Bufloc 5158. 


Figure 3 — Percentage of turbidity removal efficiency during the test. 
Source: Own Authorship (2021). 
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As can be seen in Fig. 3, the T3 treatment showed low turbidity removal efficiency 
during the first and second collections and this was due to the large amount of organic 
matter ffom Moringa oleifera. The T3 treatment only obtained considerable result with time, 
at 43 minutes, and the same treatment showed an advance in the turbidity removal during 
the course of the sedimentation test, presenting values greater than 68%. 

Itis also observed that at the end of the passage of water through the filter medium, 
the turbidity removal efficiency of this same treatment presented values above 99%, thus 
emphasizing the filter's ability to retain the high amount of organic matter present in this 
organic coagulant. 

Arantes et al. (2015) verified in their studies, using Moringa oleifera associated with 
slow filtration, removal of up to 95% for turbidity, corroborating these results. 

The T4 treatment had a rapid effect, from the first 3 minutes, where it is possible to 
identify in Fig. 3 that the turbidity removal efficiency value reached the mark of 89%, and this 
value increased during the other collections. At 43 minutes, the T4 treatment reached 96% 
of turbidity removal, surpassing the other treatments. Higashi et al. (2016) point out in their 
study turbidity removal of up to 96%, using ferric chloride in water treatment. 

Note that at the end of the test, all coagulants showed values greater than 99%, 
caused by the filter's ability to retain colloidal particles. Prado et al. (2007) observed in their 
studies using aluminum sulfate in water treatment turbidity removal efficiency higher than 
80%. Demonstrating the possibility of using organic coagulants in water treatment, since 
they presented values equal to or more expressive than inorganic coagulants. 

When the ANOVA analysis of variance, Table 1, is applied, it is observed that there is 
a significant variation in turbidity removal, both between treatments and over time, as both 
sources of variation have a P-value of less than 5% of significance, and F distribution value 
above the values delimiting the acceptance and rejection areas given by F critical. 


Source of so 


Variation gl MQ F P value Fcrit 
Time 5460,899 4 1365,225 17,87069  2,19822F-06  2,866081 
Treatment 3876,208 5 775,2417 10,14786 5,93117E-05 2,71089 
Error 1527,893 20 76,39463 
Total 10865 29 


Table 1. ANOVA analysis of variance to double factor without repetion. 
Source: Own Authorship (2021). 


Thus, the hypothesis that at least one of the sources of variation is statistically 
different from the others is validated. However, in order to determine the level of influence 


of the sources of variation, a greater volume of sample data would be necessary, a 
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recommendation for future work. 


41 CONCLUSIONS 


After the filtration process, it was possible to observe that the organic coagulant 
extracted from the Moringa oleifera seed stood out together with the ferric chloride since 
both presented results above 99% of turbidity removal. It was also observed that associating 
the coagulant made from the extraction of the Moringa oleifera seed to the water treatment 
with filtration was essential for an increase to 99% of turbidity removal efficiency. It is also 
noted that the T5 treatment showed better values than T6, evidencing the no need for the 
Bufloc 5158 auxiliary in the treatment of the water in question. 

The results for the pH parameter did not show great inequalities to the gross 
value, showing that there is no need for correction on the water treatment. The electrical 
conductivity values did not show a great difference from the gross value, except for the 
use of Moringa oleifera seed, which is due to its saline solution in its composition, showed 
values much higher than the others. 

Therefore, it demonstrated the possibility of an organic coagulant presenting 
equivalent or associative efficiency to an inorganic coagulant in the treatment of water in a 
system with filtration, being able, in this way, to bring benefits to the environment for being 
a less polluting component and being able to present greater sustainability to the process. 
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RESUMEN: El pasado 5 de marzo del 2021 
y durante ocho semanas dio início este 
experimento en San Martin Caltenco, ubicada en 
el municipio de Tochtepec en el estado de Puebla 
que consistió en separar 5 cerdos destinados 
para engorda, cada uno en un corral diferente, e 
identificados, como C1, C2, C3, C4 Y Cs, estos 
cerdos fueron alimentados con la línea Premium 
de Desarrollo y Poder Growpig!, con la finalidad 
de demostrar la hipótesis correlacional que 
dice: a mayor suministro del alimento comercial 
Poder Growpig! Y Desarrollo Growpig! en 
cerdos de engorda; mayor incremento de 
masa muscular durante su crecimiento. El 
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procedimiento experimental de validación fue 
aplicar la regresión y correlación lineal simple y 
la ración de comida suministrada fue diferente 
para cada cerdo, puesto que solo varió en cada 
uno 100 gramos y al cabo de una semana a 
todos se les aumentaba 200 gramos, así hasta 
concluir el experimento y determinar en base a 
los parámetros del coeficiente de correlación que 
fue 0.99 que en sus valores de interpretación (r) 
0.8<r< 1, indica que la correlación es muy alta 
y que efectivamente existe una correlación entre 
la cantidad de alimento suministrada de la marca 
comercial mencionada arriba y la ganancia de 
peso. La hipótesis planteada es aceptada por la 
correlación obtenida al final del proyecto. 


ANTECEDENTES 


La etapa de vida o de producción de 
los cerdos, se puede definir como un período 
de vida del animal donde necesita una 
determinada cantidad de nutrimentos para 
cumplir con sus funciones de mantenimiento y 
máxima producción. Además, tiene la capacidad 
según su desarrollo digestivo de utilizar los 
alimentos con diferentes grados de eficiencia. 
No solo tiene el efecto positivo de maximizar la 
utilización eficiente de nutrimentos y alimentos, 
sino también un efecto económico pues se evita 
un faltante o desperdício de nutrimentos que 
afecta los rendimientos de los cerdos y como 
consecuencia la rentabilidad económica. 

La alimentación eficiente de los cerdos 
es una de las prácticas más importantes, ya 
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que de ella dependen los rendimientos productivos de los cerdos, así como también la 
rentabilidad de una granja. La alimentación representa entre un 80 a un 85% de los costos 
totales de producción. Por esta razón es importante conocer ciertos conceptos importantes 
relacionados con la alimentación eficiente de los cerdos (Campabadal-Carlos; 2009). 
El cerdo posee un gran poder digestivo y de asimilación, por lo que, de acuerdo con el 
alimento suministrado, será su capacidad de aumento de peso y conversión alimenticia 
(Grepe, 2001; Roldán y Durán, 2006). 

El crecimiento incluye uno o más de los tres procesos siguientes: multiplicación 
celular, ampliación celular e incorporación de material extraído del medio ambiente (Brody, 
1945; citado por Gu et al., 1992). Hammond (1966) sistematizó los conocimientos sobre 
crecimiento y desarrollo, y definió al crecimiento como: “El aumento de peso vivo (PV) 
experimentado por un individuo desde su concepción hasta su estabilización en la edad 
adulta”; e indicó que durante el crecimiento ocurren dos fenómenos: 1) el crecimiento en sí, 
donde el animal va aumentando de peso, hasta alcanzar el tamafio adulto o madurez; y 2) 
el desarrollo, donde el animal modifica su conformación corporal, alcanzando la plenitud de 
sus funciones y facultades. La medida de crecimiento más usual en las distintas especies 
es la variación del peso corporal en el tiempo. Sin embargo, en el caso particular del cerdo, 
el proceso de crecimiento se juzga en términos de masa corporal de tejido magro y no a 
través de la masa total, debido a que los niveles de tejido adiposo en el cerdo maduro son 
variables y sumamente dependientes de la nutrición y del genotipo (Whittemore, 1996). 

Diversos factores afectan los requerimientos de los cerdos para un nutriente 
específico. Los requerimientos están influenciados por una combinación del potencial de 
crecimiento y el consumo voluntario, los cuales requerirán cambiar la concentración del 
nutriente en la dieta y conocer los requerimientos de los cerdos en base a una cantidad por 
día. Algunos de estos factores son (Dritz et al., 1997): 


* | Medio ambiente (temperatura, humedad, instalaciones y competencia por el 
alimento). 


* | Raza, sexo y origen genético de los cerdos. 

*— Estatus de salud de la piara. 

* | Presencia de hongos, toxinas o inhibidores en la dieta. 

* Disponibilidad y absorción de los nutrientes en la dieta. 

*— NVariabilidad del contenido de nutrientes y disponibilidad en el alimento. 
* Nivel de aditivos alimenticios o promotores de crecimiento. 

* | Concentración energética de la dieta. 


. Nivel de alimentación como son la alimentación limitada contra la alimentación 
a libre acceso. 
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La producción de cerdos no sólo ha avanzado en la obtención de líneas genéticas 
más precoces, con mejores índices de conversión de alimento, sino también y hacia la 
obtención de cerdos con canales mucho más magras. Este avance ha sido motivado 
primeramente a la necesidad de incrementar los rendimientos obtenidos en el desposte de 
las canales, de manera de proporcionar una mayor cantidad de carne con la consecuente 
mejora en la rentabilidad, y a una mayor demanda de este tipo de carnes en base a 
exigencias relacionadas con la salud de los consumidores. Adicionalmente a la mejora 
genética, se han desarrollado aditivos no nutricionales, los cuales son capaces de inducir 
adicionales respuestas en la calidad de las canales, como por ejemplo el compuesto 
ractopamina (RAC), un agonista B-adrenérgico de la familia de las fenoletanolaminas, 
que actúa sobre los receptores B-adrenérgicos de las células adiposas y del músculo 
esquelético, promoviendo la lipólisis, con el consecuente incremento del magro en la canal 
(Garcia-Torres et al, 2001; García-Rojas et al., 2007). 


VENTAJAS DEL USO DE LA RACTOPAMINA 


Mejora en el rendimiento en canal 


Aumenta el peso del canal: 2 — 3 puntos porcentuales 
Disminuye la grasa dorsal: 0.7 a 1 mm 
Aumenta la profundidad del lomo: 0.8-1.2mm 


GANANCIAS ADICIONALES CON EL USO DE LA RACTOPAMINA 


Mayor volumen de carne producida. (Aumenta los Kg carne/cerdo:2.5 -3.5 Kg (1,2). 

Menor uso de instalaciones. 

Mejora la rentabilidad. 

Chae et al. (1997), destacan el hecho de que el alimento en pellet es usado 
extensamente en producciones comerciales de cerdos. También se reporta en diversos 
estudios que las dietas en forma de pellet mejoran los promedios de ganancia diaria de 
peso, en comparación a dietas harinosas (Wondra et al., 1992, Aravena-Utreras; 2013) . 

GROWPIG! PODER. Es un alimento para crecimiento de cerdos en presentación 
de pellet, que ayuda a consumirlo adecuadamente y aprovechar todos los nutrimentos que 
su formulación provee. Indicado parara cerdos durante las últimas 4 semanas de engorda 
(aproximadamente en un rango de peso entre 75y 100 kg. 

Esta formulado y elaborado con proteínas de primera calidad, granos de cereales 
procesados, minerales, vitaminas y aditivos que proporcionan los nutrimentos esenciales; 
con alta biodisponibilidad para una rápida y eficiente finalización de los cerdos Dieta 
balanceada finalizadora con agradable sabor a frutas que ayuda a garantizar la óptima 
salida al mercado de cerdos en engorda. 
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Contiene Ractopamina, que ayuda a mejorar el desempefio productivo de los cerdos 
y el rendimiento y calidad del canal. 


JUSTIFICACIÓN 


La educación virtual provocada por la pandemia en el afo 2020 también fue una 
ventana de oportunidad para lograr el trabajo práctico de cada uno de los estudiantes en 
sus lugares de origen. Mientras el facilitador aportaba con el apoyo de una pantalla los 
conocimientos de herramientas estadísticas que permitieran evaluar los diferentes tipos 
de hipótesis; los educandos en su comunidad ejidal tenían la oportunidad de hacerlo de 
manera práctica; éste es uno de los casos. Algunas veces inapropiadamente se utiliza el 
termino correlacionar; sin tener la herramienta estadística básica que se encarga de manera 
objetiva de lograr ese cometido. El trabajo de nutrición de cerdos que se desarrolló durante 
8 semanas permitió aplicar y evaluar un alimento comercial Premium de Growpig!, con éste 
se varió la cantidad de alimento (variable independiente) y con esto ver sus efectos en la 
ganancia de peso de los cerdos (variable dependiente). El solo hecho de Ilevar a la práctica 
los conocimientos teóricos estadísticos que se han compartido de manera virtual genera 
una justificación hacia el presente trabajo 


OBJETIVO 


1.-Determinar la ganancia de peso obtenida a partir del suministro de alimento en 
diferentes cantidades 


2.-Evaluar la correlación que existe entre el consumo de alimento de la línea 
Premium de Growpig! para cerdos de engorda y peso obtenido al finalizar el proyecto 
aplicando la regresión lineal simple y el coeficiente de correlación lineal. 


HIPÓTESIS 


A mayor suministro del alimento comercial Poder Growpig! Y Desarrollo Growpig! 
en cerdos de engorda, mayor incremento de masa muscular durante su crecimiento. 


METODOLOGÍAS 


El experimento se realizó en la comunidad de San Martin Caltenco, ubicada en el 
municipio de Tochtepec en el estado de Puebla La localidad se encuentra a una altura de 
1980 metros sobre el nivel del mar. Con una duración de ocho semanas que va del día 05 
de marzo y concluyo el dia 30 de abril de 2021; se evaluó la ganancia de peso en cinco 
cerdos de engorda, previamente separados y con un lugar acondicionado para cada uno, 
variando la cantidad de alimento de engorda PODER GROWPIG suministrado al azar para 
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cada uno de los cinco cerdos, la cantidad suministrada a cada uno de los cerdos tenía una 
diferencia de 100 gramos, y cada semana a todos se le aumentaba a su ración 2009. 

Durante las siguientes 7 semanas el suministro de alimento aumentaba para todos 
200g por semana. En la semana 4 se les iba cambiando de manera gradual mezclando 
ambos alimentos esto con la intención de evitar alguna diarrea causada por el cambio 
brusco de alimento. El alimento siguiente fue PODER GROWPIG! de la misma línea. La 
cantidad total de alimento suministrada siguió siendo la misma 

Para este experimento se ocuparon cinco cerdos, los cuales fueron separados cada 
uno en diferente corral debidamente acondicionado para su estancia, fueron alimentados 
con la línea DESARROLLO y PODER GROWPIG! de la marca MALTA CLEY TON, teniendo 
un libre acceso al agua. 

Se les hacía medición de pesos los días viernes de cada semana, para esto se 
utilizó una báscula ganadera de la marca Braunker; los pesos tomados durante cada 
semana, el registro de porciones de alimento los capture en una bitácora de datos en el 
programa Excel y por ultimo para el análisis de datos se aplique la regresión lineal simple 
y el coeficiente de correlación. 

Procedimiento del análisis de datos mediante regresión lineal según Rocha y Avila 
(2017): 

Fórmula 1: Regresión lineal simple. 


y=a+b(x) 
Fórmula 2: Pendiente de regresión lineal 
p=? Es de Dé 
NE =(50)" 


Fórmula 3: Constante de regresión lineal. 


ad E = bl) 


Fórmula 4: Valor promedio de Y. 


=2Y 


F 
Pa 


Fórmula 5: Numerador de la fórmula del coeficiente r 


Sce-r=50-vr?-5o-ro? 


Fórmula 6. Fórmula del coeficiente de correlación r 
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Lc — Yp)? 
DO —Yp)? 


Donde: 

Y= Variable dependiente. 

X= Variable independiente. 

b= Pendiente = Coeficiente de correlación lineal. 
a= Constante de regresión lineal. 

n= Número de datos en la tabulación. 
>=Símbolo de sumatoria. 

Yc = Valor calculado de Y 

Yo = Valor de Promedio de Y. 


RESULTADOS 


Como se puede observar en tabla 1y tomando como referencia el cerdo 1; su peso 
inicial fue de 51.9 kilos y al final de la última semana de engorda su peso logrado alcanzó 
los 99.6 kilogramos; es decir un incremento por encima del 52%. 


Cerdo | Semana 1 Semana 2 Semana 3 Semana 4 Semana 5 Semana 6 Semana 7 Semana 8 
Con/ | P. Con/ | P Con/ | P Con/ | P Con/ | P Con/ | P Con/ | P Con/ | P 
día (Kg) dia (Kg) día (Kg) día (Kg) día (Kg) dia (Kg) día (Kg) día (Kg) 
(Kg) (Kg) (Kg) (Kg) (Kg) (Kg) (Kg) (Kg) 
C1 2 51.9 2.2 58.4 2.4 65 2.6 71.9 2.8 78.7 3 85.6 3.2 92.5 3.4 99.6 
c2 1.9 52.3 21 58.7 2.3 65.3 2.5 724 2.7 78.9 2.9 85.7 3.1 92.5 3.3 99.6 
c3 1.8 52.8 2 59.1 2.2 65.6 24 72.3 2.6 79 2.8 85.7 3 93.4 3.2 100 
c4 2.2 52.9 24 59.6 2.6 66.4 2.8 734 3 80.4 3.2 87.5 3.4 94.6 3.6 102 
C5 21 53 2.3 59.6 2.5 66.3 2.7 73.3 2.9 80.1 3.1 87.1 3.3 94.1 3.5 101.2 


Tabla 1: Pesos obtenidos en Kg/semana y su Consumo/dia en Kg. 
En la tabla 2 se puede encontrar el consumo promedio de alimento en kilogramos 


de los 5 cerdos por un lapso de 8 semanas; así como su peso promedio ganado en el 
tiempo de referencia con las mismas unidades de medida. 
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No do Samanas Consumo EA ad semana do a 

1 2 52.6 

2 2.2 59.12 
3 24 65.78 
4 2.6 72.62 
5 2.8 79.44 
6 3 86.34 
7 3.2 93.43 
8 3.4 100.48 


En la gráfica 1 se puede observar el comportamiento lineal positivo entre la variable 


X y la variable Y. Además se puede apreciar la ecuación de regresión generada con su 


instrumento de validación r con un valor de 0.99; lo que nos indica una correlación muy 


buena entre el alimento suministrado y el peso ganado por los cerdos. 


El resultado obtenido del coeficiente (r) 


entra dentro del rango 0.8< r <1, esto 


quiere decir que existe una correlación MUY ALTA entre el consumo de diario de alimento 


investigado y la ganancia de peso 


PESO PROMEDIO 
Kg/SEMANA 
8 


REGRESIÓN LINEAL Y CORRELACIÓN r 


y=34.23 
r= 


0.5 1 1.5 2 Z 


RA 


5 


9x-16.218 q 
0.99 e 
e 
“o 
3 3.5 4 


CONSUMO PROMEDIO Kg/SEMANA 


Gráfica 1: Ecuación de regresión lineal y coeficiente de correlación obtenidos 


Otro resultado que compara el éxito del experimento es lo que se puede apreciar 


en la tabla 3; tomando como referencia los pesos de la semana 1 y comparando con la 


semana número ocho; se ha logrado un incremento de peso por encima del 52%. Esto en 


términos de productividad es bastante bueno. 
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No de cerdo Peso en Kg semana 1 Peso en Kg semana 8 % de incremento de 

peso 

Ci 51.92 99.6 52 

c2 52.85 99.64 58 

c3 52.91 100 53 

c4 52.98 102 52 

c5 52.60 101.2 52 
PROMEDIO = 52.5% 


DISCUSIÓN 


Estadísticamente la investigación fue un éxito; pues la correlación logrado es casi 
perfecta es decir; la linealidad lograda es perfectamente coincidente con los pesos iniciales 
que promediaron 52.6 kg y los finales con una media de 100.48 kg. Un peso ganado en 
promedio para los 5 cerdos de 47.88 Kg en los 2 meses del experimento es decir; en cada 
uno de los días los cerdos ganaron en promedio 800 gms. 

Relacionado con la venta; estos se acomodaron en el mercado local de Puebla a 
$43.00/ Kg logrando un ingreso de $ 21 603.00 pesos; muy por encima de los $3000 pesos 
de Inversión que se erogaron en la compra de los lechones; que fue a razón de $600.00 
pesos por lechón; mismos que se adquirieron al inicio del semestre Enero —Junio de 2021 
pensando en el trabajo práctico que se evaluaría en la parte media del semestre. 


CONCLUSION 


El sistema de alimentación de los lechones y de los cerdos en desarrollo y engorde 
es un factor muy importante sobre los rendimientos productivos de los cerdos (Church y 
Pond; 1996). La línea de alimentos evaluada que fueron Desarrollo y Poder GROWPIG! 
De la marca comercial Malta Cleyton nos arrojaron en promedio una ganancia de peso de 
primera semana a octava semana por encima del 52%, en términos de kilogramos logrados 
son en promedio arriba de 49 kilogramos es decir; una ganancia de más de 6 kilos por 
semana. Por lo que la hipótesis planteada se acepta, ya que existe una correlación MUY 
ALTA entre el consumo de diario de alimento en (kg) y la ganancia de peso en (kg). 
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RESUMEN: Enesteestudiosecomparóelviruloma 
de una cepa atenuada de Babesia bigemina, la 
cual se ha mantenido bajo condiciones de cultivo 
in vitro durante varios aos en el laboratorio, con 
reducida virulencia para el hospedero bovino, e 
incapacidad de ser transmitida por la garrapata 
vector, con el viruloma de una cepa virulenta de 
B. bigemina transmitida en México por garrapatas 
Rhipicephalus microplus a bovinos. Se realizó un 
mapeo de genes de virulencia contra un genoma 
de referencia: Cepa Bond de Australia. El análisis 
de los resultados preliminares obtenidos con la 
herramienta de búsqueda de alineación local 
básica (BLAST) demostró que de 27 genes 
de virulencia descritos para Babesia spp. e 
identificados en la cepa mexicana virulenta de B. 
bigemina transmitida por la garrapata, solo cinco 
fueron completamente identificados en la cepa 
atenuada de laboratorio. En todos los casos, 
los porcentajes de cobertura e identidad de los 
genes identificados en la cepa virulenta de B. 
bigemina fueron más altos que los identificados 
en la cepa atenuada de laboratorio. Este hallazgo 
está supuestamente asociado con la pérdida 
parcial continua de genes de virulencia en la 
cepa de laboratorio después de múltiples pasajes 
de la población de parásitos bajo condiciones 
óptimas de cultivo in vitro. La pérdida de algunos 
de los factores de virulencia podría reflejarse en 
la ausencia de síntomas de la enfermedad en 
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ganado inoculado con la cepa atenuada, a pesar de la presencia de infección en los eritrocitos 
delhospedero bovino. Esto, tendrá que ser corroborado por el análisis transcripcional en futuros 
estudios de investigación. El análisis del repertorio de genes expresados diferencialmente en 
la cepa atenuada de B. bigemina y estudios del transcriptoma están en curso para determinar 
si se mantiene la expresión génica diferencial o si varía con respecto a la cepa virulenta. 

PALABRAS CLAVE: Babesia bigemina; cepa virulenta; cepa atenuada; genes de virulencia. 


COMPARATIVE GENOMIC STUDY OF ATTENUATED AND VIRULENT STRAINS 
OF Babesia bigemina 


ABSTRACT: In this study, the viruloma of an attenuated strain of Babesia bigemina, which 
has been maintained under in vitro culture conditions for several years in the laboratory, 
with reduced virulence for the bovine host, and inability to be transmitted by the tick vector, 
was compared with the viruloma of a virulent strain of B. bigemina transmitted in Mexico by 
Rhipicephalus microplus ticks to cattle. Virulence gene's mapping was performed against a 
reference genome: B. bigemina Bond strain from Australia. Analysis of preliminary results 
obtained with the Basic Local Alignment Search Tool (BLAST) showed that out of 27 virulence 
genes described for Babesia spp. and identified in the virulent Mexican tick-borne strain of 
B. bigemina, only five were fully identified in the laboratory attenuated strain. In all cases, 
the percentages of coverage and identity of the genes identified in the virulent strain of B. 
bigemina were higher than those identified in the laboratory attenuated strain. This finding 
is putatively associated with continued partial loss of virulence genes in the laboratory strain 
after multiple passages of the parasite population under optimal in vitro culture conditions. 
The loss of some of the virulence factors could be reflected in the absence of symptoms of 
the disease in cattle inoculated with the attenuated strain, despite the presence of infection 
in the erythrocytes of the bovine host. This will have to be corroborated by transcriptional 
analysis in future research studies. Analysis of the differentially expressed gene repertoire 
in the attenuated strain of B. bigemina and transcriptome studies are ongoing to determine 
whether differential gene expression is maintained or varies from the virulent strain. 

KEYWORDS: Babesia bigemina; virulent strain; attenuated strain; virulence genes. 


11 INTRODUCCIÓN 


La babesiosis del ganado es una enfermedad transmitida por garrapatas de elevada 
importancia socioeconómica, la cual es causada por protozoarios Apicomplexa del género 
Babesia considerados como parásitos intraeritrocíticos obligados (Bock et al., 2004). 
En México, la babesiosis bovina es causada por Babesia bigemina y Babesia bovis. La 
enfermedad se distribuye en zonas tropicales y es transmitida por garrapatas Rhipicephalus 
microplus y R. annulatus (Rojas et al., 2018). La presentación de la babesiosis bovina 
causada por Babesia bigemina (B. bigemina) es 12-16 días después de que las garrapatas 
se alimentan del bovino, coincidente con la presencia del parásito intraeritrocítico en sangre 
periférica. La temperatura se eleva paralelamente al incremento en parasitemia hasta 41- 
42ºC en2 03 días. Los animales infectados están tristes, anoréxicos y con pelo hirsuto. Se 
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presenta hemoglobinemia y hemoglobinuria seguida por palidez de las mucosas junto con 
otros signos que pueden incluir constipación, deshidratación, temblor muscular, debilidad, 
postración y, si no se instaura el tratamiento adecuado, la muerte (Dalgliesh et al., 1981). 
Se han generado diversos estudios sobre el bovino hospedero, el agente y sus 
relaciones tratando de reconocer sus características biológicas, así como su distribución, 
con el objeto de establecer programas de prevención y/o control de la enfermedad. Estas 
actividades han requerido la replicación del ciclo biológico logrado mediante el crecimiento 
artificial de estos parásitos en condiciones in vitro (Levy y Ristic, 1980; Vega, et al., 1985). 
Se han reportado diversos métodos para disminuir la virulencia de las especies patógenas 
en bovinos: Pases rápidos por inoculación parenteral en becerros esplenectomizados 
(Callow et al., 1979), pases lentos por inoculación parenteral en becerros intactos (Dalgliesh 
et al., 1981), y la metodologia de cultivo in vitro para las especies de Babesia que afectan 
a los bovinos (Levy y Ristic, 1980; Figueroa et al., 1984; Vega et al., 1985; Fish et al., 
2008). La atenuación de patógenos intracelulares por cultivo in vitro es un fenómeno bien 
conocido y patógenos como B. bovis y B. bigemina se han tornado avirulentos o menos 
patógenos mediante prolongado crecimiento en cultivo celular (Yunker et al., 1987; Kutiler 
et al., 1988; Schuster, 2002). Trabajos de nuestro grupo de investigación han permitido 
demostrar que una cepa de B. bigemina mantenida por pases continuos en cultivo in 
vitro, se ha comportado como población atenuada al no afectar en forma importante los 
valores hematológicos en animales inoculados (Hernández et al., 1990; Figueroa et al, 
1998). La cepa atenuada de B. bigemina ha inducido protección al desafio heterólogo con 
sangre infectada o con garrapatas infectadas bajo condiciones controladas y de campo 
(Hernández et al. 1990; Figueroa et al, 1998; Álvarez et al, 2004). Las cepas vacunales 
han sido desarrolladas para conferir una adecuada protección sin ser transmitidas por 
garrapatas, un rasgo considerado ecológicamente deseable ya que previene la aparición 
de casos clínicos debido a la transmisión de la cepa vacunal por garrapatas (Mangold et 
al., 1993). Esta capacidad ha sido demostrada después de un pase en ganado susceptible 
(Dalgliesh et al., 1983; O'Sullivan y Callow, 1966; Hernández et al., 1990; Timms et al, 
1990; Mangold et al., 1993). Además, este tipo de vacunas vivas atenuadas tiene la gran 
ventaja conferida por su producción in vitro: un bajo riesgo de contaminación con otros 
agentes infecciosos (Pipano, 1995; Rojas et al, 2018a; 2018b); una suspensión de glóbulos 
rojos altamente parasitada, y la posibilidad de producción de la vacuna en relativamente 
gran escala (Timms y Stewart, 1989). Se demostró que la cepa mexicana de Babesia 
bigemina atenuada no fue transmitida por garrapatas después de varios pases en bovinos 
susceptibles. Se demostró ausencia de reversión a la virulencia de esta cepa después de 
tres pases sucesivos en bovinos susceptibles y se confirmó que los parásitos mantenidos 
en cultivo in vitro perdieron la capacidad para multiplicarse en la garrapata vector durante 
un segundo pase por sub-inoculación en bovinos susceptibles (Rojas et al, 2011). En este 
estudio se confirmó que los parásitos mantenidos en cultivo in vitro perdieron la capacidad 
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para multiplicarse en la garrapata vector durante un segundo pase por sub-inoculación en 
bovinos susceptibles. Esto, evidenciado por la ausencia de quinetos en la hemolinfa de 
garrapatas alimentadas en bovinos con parasitemia patente. Además, se demuestra que 
aun cuando es sujeta a tres pases sucesivos en bovinos altamente susceptibles, la población 
de parásitos mantiene su característica biológica de reducida virulencia manifestada por 
los valores clínicos determinados en los bovinos sub-inoculados, los cuales no mostraron 
cambios relevantes. 

Por otro lado, se desconocen hasta el momento los componentes del parasito 
responsables de la inducción de la hemolisis y otros efectos fisiopatológicos en los bovinos 
infectados. Resulta evidente que se requiere de estudios adicionales e investigación 
de frontera con metodologías de nueva generación para determinar los mecanismos 
involucrados en la determinación de virulencia para el desarrollo del parásito dentro del 
eritrocito infectado (y provocar la enfermedad en un bovino inoculado), así como para 
el desarrollo de estadios de Babesia subsecuentes y necesarios para la infección del 
intestino de Rh. microplus (formas sexuales) que permiten al parasito adquirir nuevamente 
la capacidad infectiva hacia la garrapata vector y, consecuentemente, ser transmitidos de 
formatrans-ovárica. Solo a través de la utilización de innovadores métodos de secuenciación 
masiva y el análisis del genoma y del transcriptoma, se podrá probar la hipótesis que 
algunos cambios genéticos (rearreglos cromosomales, deleciones o inserciones génicas) 
han ocurrido en el genoma de la cepa atenuada sujeta de estudio y mantenida en cultivo in 
vitro, los cuales se ven reflejadas en un fenotipo atenuado y carente de la capacidad de ser 
transmitido por la garrapata vector R. microplus. El análisis comparativo del genoma de la 
cepa virulenta mexicana con el genoma de la cepa virulenta australiana denominada Bond 
y de otras cepas (Jackson et al, 2014), permitirá demostrar que los cambios observados 
con las cepas mexicanas no se derivan de polimorfismos en genoma asociados al origen 
geográfico de las cepas. La genómica comparativa permite identificar genes únicos de 
ciertas cepas patógenas los cuales pueden estar ausentes en las cepas apatógenas. 

La patogenicidad de Babesia para el ganado bovino está determinada en parte 
por la interacción con el sistema inmunológico del huésped y la presencia de genes de 
virulencia del parásito. En este estudio, se comparó el viruloma de una cepa atenuada de B. 
bigemina, la cual se ha mantenido bajo condiciones de cultivo in vitro durante varios afios 
en el laboratorio, con reducida virulencia para el hospedero bovino, e incapacidad de ser 
transmitida por la garrapata vector, con el viruloma de una cepa virulenta de B. bigemina 
transmitida por garrapatas Rhipicephalus microplus a bovinos. 
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21 MATERIALES Y MÉTODOS 


2.1 Parásitos 


La cepa atenuada de B. bigemina es una población de parásitos originalmente 
derivada de un aislado virulento de México recolectado de un caso clínico de babesiosis. La 
cepa se adaptó a un cultivo in vitro usando un sistema estacionario microaerófilo (Vega et 
al, 1985). Una vez establecida en cultivo in vitro se cultivó continuamente durante casi cinco 
afios consecutivos en laboratorio con pases cada 72 o 96 hrs (al menos 500 subcultivos en 
cultivo). Posteriormente, el cultivo in vitro de la B. bigemina se realizó de forma discontinua 
y con un número indefinido de pasadas en cultivo; se cultivó por períodos cortos de tiempo 
(uno o dos meses), justo hasta obtener material suficiente para los experimentos que se 
realizaron durante el desarrollo de la vacuna viva atenuada (períodos 1991-1998 y 2001- 
2006) y se congeló inmediatamente en líquido nitrógeno hasta que se realizaron otros 
experimentos en el período 2011-2017 en los que se cultivó para los últimos experimentos 
(transmisibilidad de garrapatas) y validación de inmunogenicidad (Figueroa et al, 1998; 
Álvarez et al, 2004; Rojas-Martínez et al, 2018a:2018b; Rojas-Ramírez et al, 2011) hasta 
ahora, cuando se realizó la secuenciación. Más importante para tener en cuenta es que 
el material biológico atenuado no ha pasado a través de un hospedero bovino desde 
su aislamiento original a diferencia del aislado virulento (del que se seleccionó la cepa 
atenuada) y que se replica en animales siempre que se requiere. Morfológicamente, la cepa 
atenuada no ha cambiado; sigue infectando eritrocitos tanto in vivo como in vitro, pero no 
provoca enfermedad cuando se inocula en el ganado y se ha mantenido en forma alterna 
en cultivo continuo y crioconservación desde entonces (Figueroa et al, 1998; Álvarez et al, 
2004; Rojas-Martínez et al, 2018a;2018b; Rojas-Ramírez et al, 2011). La cepa virulenta, 
originalmente aislado de un caso clínico de campo, se ha mantenido a través de pasajes 
de garrapatas en susceptibles animales y criopreservación en nitrógeno líquido (Vega et al, 
1985; Figueroa et al, 1998; Rojas-Ramírez et al, 2011). 


2.2 Extracción de ADN genómico 


Se utilizaron eritrocitos bovinos infectados con las diferentes poblaciones de B. 
bigemina para llevar a cabo la extracción de ADN genómico por métodos convencionales 
(Strauss 2001; Genis et al, 2008). Se utilizaron entre 20 y 30 ug de ADN genómico para 
preparar las dos bibliotecas clonales requeridas para cada cepa. 


2.3 Secuenciación del genoma de B. bigemina, cepas atenuada y virulenta 


La secuenciación del genoma completo se realizó con la plataforma Illumina de 
secuenciación de próxima generación Mi-Seqg (Ilumina, San Diego, CA, EE. UU.) en el 
Instituto de Biotecnologia de la Universidad Nacional Autónoma de México ubicado en 
Cuernavaca, Morelos. Las lecturas de secuencias obtenidas con el sistema Illumina se 
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ensamblaron de novo utilizando el Programa SPAdes (versión 3.1.1, Laboratorio de Biologia 
Algorítmica, Universidad Académica San Petersburgo, San Petersburgo, Rusia) (Bankevich 
et al, 2012). Los contigs generados se ordenaron con la herramienta MUmmer (Kurtz et 
al, 2004). El esfuerzo de secuenciación del genoma de B. bigemina permitió obtener una 
cobertura de 264 X y 89 X para la cepa atenuada y virulenta, respectivamente. La lista de 
genes seleccionados genes se basó inicialmente en la consideración de artículos anteriores 
sobre inmunología clásica y estudios del transcriptoma donde se menciona que son genes 
atribuidos a la virulencia o posibles factores de virulencia (Bastos et al, 2013; Eichenberger 
et al, 2018). El mapeo de secuencias de genes virulentos contra un genoma de referencia 
(genoma australiano de B. bigemina) se realizó con BLAST (Altschul, 1990). Esto permitió 
un primer análisis de la composición del genoma de cada cepa e identificó posibles 
cambios (deleciones o inserciones) o reordenamientos cromosómicos que ocurrieron en 
los genomas de las cepas virulentas y atenuadas de B. bigemina. 


31 RESULTADOS Y DISCUSIÓN 


La secuenciación del genoma y el análisis del ensamble mostraron que la cepa 
atenuada de B. bigemina contenía 9,180,241 pb, mientras que la cepa virulenta 11,852,459 
pb, una diferencia de 22.54% (Cuadro 1). El esfuerzo de secuenciación obtuvo una cobertura 
de 264X y 89X para las cepas atenuada y virulenta, respectivamente. Las secuencias de 
genes de virulencia se obtuvieron a partir de una búsqueda BLAST y las diferencias entre 
las cepas fueron notables; mientras que la cepa virulenta de B. bigemina transmitida por la 
garrapata presentó 27 genes con un promedio de cobertura del 98.06%, la cepa atenuada 
sólo presentó cinco de esos genes, con 81.2% de cobertura (Cuadro 2). En todos los casos, 
los porcentajes de cobertura e identidad de los genes identificados en la cepa virulenta 
de B. bigemina fueron más altos que los identificados en la cepa atenuada de laboratorio 
(Cuadro 3). Este hallazgo virulencia en la cepa de laboratorio después de múltiples pasajes 
de la población de parásitos bajo condiciones óptimas de cultivo in vitro. 

Las secuencias de los genes de virulencia se obtuvieron de una búsqueda de 
alineación local básica con la herramienta BLAST, y las diferencias entre las cepas fueron 
notables; mientras que la cepa virulenta presentó 27 genes con una cobertura promedio 
de 98.06%, la cepa atenuada solo presentó cinco de esos genes con 81.2% de cobertura 
(Cuadro 2). Eltipo más común de genes identificados en ambas cepas fueron genes parciales 
putativos; 18 para el tipo virulento y tres para la cepa atenuada. El único gen hipotético 
identificado con la secuencia más corta observada y sin diferencias en las bases pertenece 
a la cepa virulenta (107 pb; Cuadro 2). El gen más grande (6117 pb) presente en la cepa 
virulenta mostró el mayor número de diferencias en bases (216 pb; Cuadro 2). Se hicieron 
intentos para predecir si esos genes se traducirían en proteínas. Para las predicciones de 
los genes eucariotas, generalmente se utiliza el programa Augustus, que predice el marco 
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de lectura abierta de las proteínas en las seis posibles formas. Desafortunadamente para 
nosotros, el ensamblaje del genoma se realizó con la secuencia de lecturas aún conteniendo 
demasiados contigs, lo que dificultó la realización de un adecuado análisis bioinformático. 
Esto se resolverá mediante la realización de genoma de alto rendimiento secuenciación 
en una plataforma de secuenciación de lecturas mas largas como la denominada “PacBio 
Single Molecule Real Time (English et al, 2012) o la tecnología de Oxford Nanopore, con el 
objeto de cerrar huecos o extender los extremos de los contigs obtenidos en el ensamblaje 


de Illumina. 


Datos de Secuenciación | Cepa B. bigemina atenuada B. bigemina virulenta 


1,537 


*% Contigs 


4,017 
Longitud total genoma 9,272,266 11,852,459 
Contenido GC (%) 52.3 50.74 
Cobertura (X) 264 89 


GenBank No. de acceso PRJNA685856 PRJNA685857 


Cuadro 1. Cepas de Babesia bigemina analizadas en el estudio. 


Yo 


identidad Valor E 


No. De Acceso 


Identificación de genes 


No. de Contig 


Cepa virulenta de Babesia bigemina 


calcium-dependent protein kinase 4 XM 012911530.1 contig00169 99.08 0 
calmodulin-domain protein kinase 2 XM 012910710.1 | contig00134 98.75 0 
CAMP-dependent protein kinase XM 012914446.1 contig00026 98.88 0 
casein kinase | XM 012913338.1 contig00174 98.62 0 
cGMP dependent protein kinase XM 012914849.1 | contig00056 99.21 0 
eyclin 4 XM 012914402.1  contig00018 98.62 (o) 
diphosphate kinase family, putative XM 0129137591 contig00012 99.34 1E-166 
dnaJ C terminal region domian XM 012910446.1 contig00220 98.74 0 
glycerol kinase XM 012914443.1 contig00026 98.99 4E-160 
glycogen synthase kinase-3 alpha XM 012910925.1 ' contig00076 98.86 4E-141 
Hypotetical gene XM 012911455.1 contig00674 100 2E-54 
merozoite surface glycoprotein AF298630.1 contig00098 92.31 2E-106 
mitogen-activated protein kinase XM 0129124991 contig00013 98.23 0 
phosphatidylinositol 3-and 4-kinase XM 012912148.1 | contig00156 97.66 0 
phosphatidylinositol 4 kinase XM 012911128.1 contig00059 97.19 0 
pa ei XM 0129144971 contigooo6s 9910 0 
ee XM 012910512.1 contigoo138 96.47 ) 
Probable fructokinase XM 012910446.1 contig00220 98.74 0 
protein kinase domain XM 012914171.1 contig00006 98.32 0 
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putative rhoptry protein NC 027216.1 contig00167 96.66 0 
RAP-1 related antigen NC 027216.1 contig00003 97.46 0 
related serine/threonine protein kinase XM 012911485.1 | contig00028 95.70 0 
Ser/Thr protein kinase XM 012910695.1 ' contig00064 96.96 0 
serine/threonine kinase XM 012914264.1 contig00171 98.78 0 
serine/threonine kinase 1 XM 012914521.1 contig00041 98.23 0 
transcription factor TFIIB XM 012912762.1 | contig00016 98.16 0 
transcription initiation factor TFIIB XM 012910684.1 contig00064 98.67 0 
98.06* 


*Cobertura promedio 


Cuadro 2. Genes de virulencia identificados en el genoma de la cepa virulenta de B. bigemina 


Identificación de Genes 


cGMP dependent protein kinase XM 012914849.1 contig00032 84.55 1E-15 

calcium-dependent protein kinase 4 XM 012911530.1 contig00066 78.68 6E-32 

phosphatidylinositol 4 kinase XM 012911128.1 contig00009 80.52 2E-46 

serine/threonine kinase XM 012914264.1  contig00035 83.54 4E-24 

transcription factor TFIIB XM 012912762.1 contig00023 78.73 2E-41 
81.20* 


*Cobertura promedio 


Cuadro 3. Genes de virulencia identificados en el genoma de la cepa atenuada de B. bigemina 


Además, considerando que las diferencias fenotípicas entre la cepa atenuada de 
cultivo y la cepa virulenta incluye la falta de transmisión por garrapatas de la cepa cultivada, 
se realizó un análisis comparativo de la codificación de genes específicos de quinetos para 
proteínas involucradas en el desarrollo de la etapa sexual, u otros genes requeridos para 
la transmisión de Babesia por garrapatas. Varios de los genes que están involucrados en 
el desarrollo del parásito en la garrapata vector están bien descritos en la literatura, tales 
como los genes de la familia hap2 y CCp (Cuadro 3) (Bastos et al, 2013; Alzan et al, 2016; 
Camacho-Nuez et al, 2017; Bohaliga et al, 2018; 2019). Se encontró que estos genes 
estaban bien conservados en ambos genomas analizados con una secuencia ligeramente 
superior identidad en la cepa atenuada. Curiosamente, los genes CCp2 y BBBOND fueron 
identificado en dos contigs diferentes y aparentemente faltaban algunos nucleótidos 
(406 pb y 61 pb, respectivamente). El gen CCp1 descrito en B. bovis (Alzan et al, 2016) 
aparentemente falta tanto en el genoma de B. bigemina de la cepa virulenta como de la 
atenuada (Cuadro 3). En cuanto al desarrollo genes descritos en el Cuadro 3 que parecían 
tener una homologia similar sin importar si pertenecían al aislado avirulento o al virulento, 
estudios previos han demostrado que los miembros de la familia del gen CCp de Babesia 
se expresan diferencialmente a lo largo del ciclo de vida y los datos presentados en esos 
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estudios demostraron que los genes CCcp de Babesia se expresan predominantemente 
durante la replicación del parásito en la garrapata vector (Bastos et al, 2013; Alzan et al, 
2016; Camacho-Nuez et al, 2017; Bohaliga et al, 2018; 2019). Si este es el caso de los genes 
CCp identificados en este estudio de genómica comparativa necesita más investigación y 
actualmente estamos preparando estudios apropiados in vivo y experimentos in vitro para 
probar eso por RT-PCR. 


41 CONCLUSIONES 


La pérdida de algunos de los factores de virulencia podría reflejarse en la ausencia 
de síntomas de la enfermedad en ganado inoculado con la cepa atenuada, a pesar de 
la presencia de infección en los eritrocitos de un hospedero bovino. Esto, tendrá que 
ser corroborado por el análisis transcripcional en futuros estudios de investigación. Será 
interesante analizar el repertorio de genes expresados diferencialmente en la cepa atenuada 
de B. bigemina y estudios del transcriptoma están en curso para determinar si se mantiene 


la expresión génica diferencial o si varía con respecto a la cepa virulenta de B. bigemina. 
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RESUMEN:  Doscientos cincuenta yaques 
(Leiarius marmoratus) de 14 meses de edadfueron 
distribuidos en cinco estaciones piscícolas del 
departamento del Meta (Colombia), para realizar 
seguimiento durante 6 meses de su crecimiento 
y condición corporal, midiendo mensualmente la 
longitud total (cm), longitud anal (cm), ancho de 
la cabeza (cm) contorno del cuerpo (cm) y peso 
corporal (Kg). La observación continuó durante 6 
meses adicionales para establecer el momento 
de su primera madurez sexual y determinar 
así la edad a la pubertad. Cuando presentaron 
por primera vez signos de madurez sexual, la 
maduración final de las gónadas, la ovulación y 
la espermiación fueron inducidas hormonalmente 
(Extracto de Hipófisis de Carpa — EHC y Ovaprim 
— OVAP), evaluando en total cinco protocolos 
de administración de las hormonas. La mayor 
ganancia de peso y de talla se presentaron en 
las fincas 2 (0,41K9g) y 5 (6,1cm), observândose 
que las hembras presentaron mayor desarrollo 
corporal que los machos. Teniendo en cuenta 
el diâmetro ovocitario post-inducción (868,0 + 
98,6 um) y las fecundidades absoluta (638.248 
+ 129.381 ovocitos/nembra) y relativa (207.479 + 
27.550 ovocitos.Kg!), la mejor respuesta ovárica 
al tratamiento hormonal se observó con los 
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protocolos basados en EHC. Igualmente, el protocolo más efectivo para aumentar el volumen 
seminal (3,3 + 0,8mL), la movilidad espermática (58,0 + 13,9 %), la duración de la motilidad 
(62,2 + 7,4 seg) y la viabilidad espermática (95,6 + 1,8 %) fue también aquel basado en EHC; 
sin embargo, la concentración espermática más alta se observó cuando se administró OVAP 
(4.430 + 915,9 x 10º espermatozoides/yL). En conclusión, bajo condiciones de cautiverio, 
los primeros signos de madurez sexual de la especie se pueden presentar entre los 22 y 24 
meses de edad, lo cual puede considerarse como la edad a la pubertad de la especie. 
PALABRAS CLAVE: Gónadas; inducción hormonal; Leiarius marmoratus; pubertad; silúridos 
nativos. 


EVALUATION OF THE REPRODUCTIVE BEHAVIOR AND DETERMINATION OF 
YAQUE'S (Leiarius marmoratus) MALES AND FEMALES PUBERTY AGE UNDER 
CAPTIVE CONDITIONS 
ABSTRACT: 250 fourteen-month-old yaques (L. marmoratus) were split in five fishing farms 
of the department of Meta (Colombia). This distribution was performed in order to monitor 
monthly their growth and body condition by measuring their total length (cm), anal length 
(cm), head width (cm), body girth (cm) and body weight (Kg) during a period of six months. 
The observation was performed for six more months so as to establish the time for their first 
sexual maturity and determine the age of puberty. Once the first signs of sexual maturity were 
observed, the gonads final maturation, ovulation and spermiation were hormonal induced 
(Carp Pituitary Extract - CPE and Ovaprim - OVAP), assessing five hormone administration 
protocols in total. The greatest gain in weight and length were observed in the fishing farms 2 
(0,41Kg) and 5 (6,1cm) where females showed greater body development than males. Based 
on the after-induction ovarian diameter (868,0 + 98,6 um) and the absolute fecundity (638.248 
+ 129.381 oocytes/female) and relative fecundity (207.479 + 27.550 oocytes.Kg'!), the best 
ovarian response to the hormonal treatment was observed in protocols under CPE. Similarly, 
the protocol under CPE was also the most effective in increasing the seminal volume (3,3 + 
0,8mL), the sperm motility (58,0 + 13,9 %), the duration of motility (62,2 + 7,4 sec.) and the 
sperm viability (95,6 + 1,8 %); nonetheless, the highest sperm concentration was seen when 
OVAP (4.430 + 915,9 x 10º spermatozoa/uL) was administered. In conclusion, the first signs of 
sexual maturity in L. marmoratus may appear around the age of 22 and 24 months in captivity 

conditions, which can be considered as the species puberty age. 
KEYWORDS: Gonads; hormonal induction; Leiarius marmoratus; puberty; reproduction; 
silurid fish. 


INTRODUCCIÓN 


El yaque (Leiarius marmoratus), es un pez reofílico suramericano que pertenece a 
la familia siluridae y habita generalmente ambientes acuáticos demersales de las cuencas 
de los ríos Amazonas, Essequibo y Orinoco (Castillo, 2001; Ramírez Gil & Martínez Ajiaco, 
1997). Su reproducción es estacional, ocurriendo ésta entre los meses de abril y junio 
(época de Iluvia), alcanzando los animales adultos una longitud total promedio de 100 cm y 
un peso corporal de 12 kg, aproximadamente. Ostenta una carne blanca muy apetecida en 
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los mercados locales debido a sus características organolépticas (textura, sabor y pocas 
espinas), convirtiêndolo en un pez de cuero de interés para la acuicultura nacional (Ramírez 
Y Ajiaco, 1997). 

Algunos trabajos preliminares realizados en el Instituto de Acuicultura y Pesca de 
la Universidad de los Llanos — IALL (Colombia), han mostrado una buena adaptación de 
la especie al confinamiento logrando adaptar paulatinamente sus hábitos alimenticios 
piscívoros (propio de la especie en ambientes naturales) (Layman Craig A., Winemiller Kirk 
O., Arrington D. Albrey, & Jepsen David B., 2005) u omnívoros (Novoa R, 2002; Ramírez 
Gil & Martínez Ajiaco, 1997) al consumo específico de dietas secas comerciales (Cruz 
Casallas et al., 2008). De igual forma han realizado investigaciones encaminadas hacia la 
estandarización de protocolos para la inducción hormonal de machos de L. marmoratus bajo 
condiciones específicas de cautiverio, obteniendo resultados favorables para la piscicultura 
(Tatiana Mira, Medina, & Cruz, 2010). 

No obstante, resulta poco probable fomentar su cultivo a escala comercial si no 
se conoce con amplitud el comportamiento reproductivo de las hembras bajo las mismas 
condiciones de cautividad, teniendo en cuenta que la reproducción de la especie puede verse 
afectada por diferentes factores tales como las condiciones climáticas y/o medioambientales 
diferentes a las que encuentran en su entorno natural, a la nutrición y al estrés producido 
por la captura y/o confinamiento en que se encuentren los animales, afectando de manera 
directa su reproducción (Sumpter et al., 1994; Pankhurst y Van der Kraak, 1997; Zohar y 
Mylonas, 2001). Por lo anterior, un grupo de investigaciones se han desarrollado con el fin 
de establecer la edad en que los peces criados en cautiverio inician su madurez sexual, 
dentro de los cuales se encuentra el realizado por (Davis, 2009) donde determinaron si la 
maduración sexual del Ictalurus punctatus es dependiente de la talla y edad del animal, o del 
número de ciclos de temperaturas anuales al que son sometidos. De igual forma, (Cubillos 
S., Canales A., Bucarey S., Rojas S., & Alarcón M., 1999) efectuaron en Chile un estudio 
evaluando la estacionalidad del período reproductivo de la sardina común (Strangomera 
bentincki) y anchoveta (Engraulis ringens) y sus cambios interanuales, como también la 
determinación de la talla de la primera madurez sexual de ambas especies. Es por ello, 
que este trabajo tiene como objetivo principal conocer el comportamiento reproductivo 
del Yaque (L. marmoratus) cuando es producido en cautiverio, con el fin de optimizar su 
periodo reproductivo. 


MATERIALES Y MÉTODOS 


Área de estudio: La investigación se Ilevó a cabo entre los meses de Agosto del 
2010 a Julio del 2011 en cinco estaciones piscícolas del Departamento del Meta (Colombia), 
así: Villa Lucy (Finca 1), ubicada en el municipio de Puerto López (178 m.s.n.m; temperatura 
media de 26,5ºC); Munuanú (Finca 2) y Aguallano (Finca 4), pertenecientes al municipio 
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de Restrepo (570 m.s.n.m; temperatura media de 25,8ºC); Alevinos del Meta (Finca 3), del 
municipio de Guamal (525 m.s.n.m; temperatura media de 26ºC); y El Paraíso (Finca 5) 
perteneciente al municipio de Cumaral (452 m.s.n.m; temperatura promedio de 21ºC). 

La evaluación ovocitaria y seminal de los animales seleccionados para reproducción 
artificial se realizó en el Laboratorio de Reproducción y Crioconservación del Instituto 
de Acuicultura y Pesca de la Universidad de los Llanos (IALL), ubicado a 12 kilómetros 
del municipio de Villavicencio (Meta), a una altura sobre el nivel del mar de 418 metros, 
pluviosidad de 4050mms anuales, humedad relativa del 75% y temperatura ambiental 
media de 27ºC. 

Material biológico y manejo de los animales: Los animales (machos y hembras) 
se obtuvieron de tres parejas reproductoras silvestres desovadas en la estación piscícola 
del IALL entre los meses de Mayo y Junio del afo 2009. Un stock de 250 animales fue 
distribuido entre las cinco estaciones piscícolas vinculadas a la investigación (50 peces/ 
finca), a los cuales se les realizó marcaje de manera individual a través de microchips 
implantados en la región dorsal de los mismos. La densidad de siembra utilizada fue de 200 
gr.m?! y su alimentación correspondió a dieta comercial con 30% de proteína bruta (PB) a 
razón del 2% de la biomasa total de cada estanque, variando la cantidad de alimento según 
el peso promedio registrado durante los muestreos mensuales. 

Variación del crecimiento y condición corporal de machos y hembras desde los 
12 hasta los 18 meses de edad: Para la evaluación del crecimiento y condición corporal 
de los animales, se realizó un muestreo mensual al 50% de los individuos (25 ejemplares 
seleccionados al azar) en cada estación piscícola, desde el mes de Septiembre/2010 a 
Marzo/2011, determinando el peso corporal (Kg) y medidas morfométricas como longitud 
total (cm), longitud anal (cm), ancho de la cabeza (cm) y contorno del cuerpo (cm) de los 
animales muestreados. Para establecer el peso individual, se utilizó una báscula electrónica 
(Jaz International Trading LTD, Colombia) con 0,001 de precisión. Entre tanto, las medidas 
morfométricas de los animales se tomaron mediante la utilización de un ictiómetro. 

Determinación de la edad a la pubertady evaluación del desarrollo y maduración 
gonadal de hembras y machos entre los 12 y 24 meses de edad: Para la determinación 
de la edad a la pubertad y evaluación del desarrollo y maduración del ovario mediante 
características sexuales secundarias, se realizaron biopsias ováricas (Garcia, 1989) a 
hembras seleccionadas aleatoriamente durante los muestreos mensuales realizados. La 
biopsia consistió en la introducción de un catéter plástico en la papila urogenital, obteniendo 
una muestra representativa de ovocitos (+70 ovocitos/hembra) y almacenándolos en un 
tubo eppendorf con cloruro de sodio al 0,9% para mantener las muestras estables hasta 
la evaluación de los mismos. Posteriormente, se dispusieron en cajas de petri para su 
observación bajo un estereoscopio binocular (Nikon SMZ 800, Japón) y determinar así 
el diâmetro ovocitario (um). Luego, los ovocitos fueron impregnados con solución Serra, 
compuesta de 85% de Alcohol, 10% de Formol y 5% Ácido Acético Glacial, la cual tiene 
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como efecto aclarar el vitelo del ovocito para poder determinar la posición de la vesícula 
germinal (central, migrando o periférica). Con base en el diámetro ovocitario y posición de 
la vesícula germinal se estimó el estadio de maduración ovárico. 

Para evaluar el estado de desarrollo y maduración testicular se tuvo en cuenta 
características sexuales secundarias como lo es la presencia o no de semen en la papila 
urogenital después de realizar un leve masaje en sentido cráneo-caudal de la región 
celómica del animal. A los animales que presentaron emisión de semen se les fue colectado 
para realizar el posterior análisis del mismo, determinando volumen seminal (ml), movilidad 
masal (%), duración de la motilidad espermática (seg), viabilidad espermática (%), 
osmolaridad (mOsm.Kg!) y concentración espermática (sptz x 10%yL) a cada individuo. 
Para la evaluación del volumen seminal se utilizaron tubos eppendorf aforados de 1,5 ml; 
la movilidad masal y la duración de la motilidad se determinaron mediante la colocación de 
20yL de semen sobre una lámina excavada (Micro-Slides Premiere, China) con 180yL de 
agua destilada (1:20) para la activación del mismo y observado a través de un microscopio 
óptico (Axiostar-Carl Zeiss, USA) con magnificación de 10X; la duración de la motilidad 
se estimó desde la mezcla del agua destilada con el semen hasta el momento en que el 
90% de las células cesaran su movimiento. La viabilidad espermática fue determinada 
mediante la realización de un extendido (mezcla de 5 yL de semen + una gota (20yL) de 
eosina-nigrosina (Sigma Chemical Co., USA)> sobre una lâmina portaobjetos (Waldemar 
Knittel, Germany) y observada a través de un microscopio óptico (Axiostar-Carl Zeiss, 
USA), magnificación 10X. La osmolaridad espermática se evaluó mediante un osmómetro 
(Osmette Ill — Precision System INC, USA), utilizando 10yL de semen. Por último, la 
concentración espermática fue determinada a través de una câmara de Neubauer de 
0,100mm de profundidad (Optik Labor, Germany). 

Evaluación a la inducción hormonal: Se realizó inducción hormonal a treinta 
y cuatro (34) ejemplares, entre hembras (2,79 + 0,1K9g) y machos (1,32 + 0,03Kg) que 
mostraron signos secundarios de madurez gonadal (enrojecimiento y protrusión de la 
papila urogenital en hembras y emisión de semen en machos luego de realizar masaje 
abdominal). Estos ejemplares fueron trasladados 24 horas antes de dar inicio a la 
inducción hormonal a piletas circulares de cemento con 2,0mº de agua, disponiendo de 
aireación constante. Para la inducción a la maduración gonadal se utilizaron dos inductores 
hormonales, Extracto de Hipófisis de Carpa (EHC) y una mezcla comercial de análogo 
superactivo de Hormona Liberadora de Gonadotropina + Bloqueador de los receptores de 
D2 de Dopamina - OVAPRIMº (OVAP), mediante aplicación intramuscular (IM) en la región 
dorsal de los animales. 

Los protocolos utilizados para las hembras fueron: EHC = 2,75mg.Kg'' (0,25mg. 
Kg” (0h), 0,5mg.Kg" (9h), 2mg.Kg! (18h)); OVAP = 0,5ml.Kg! (Única inyección). Entre 
tanto, los machos seleccionados se indujeron mediante aplicación de tres protocolos bajo 
dos inductores hormonales: EHC' = 2,7mg.Kg! (0,2mg.Kg"! (Oh), 0,5mg.Kg"" (9h), 2mg. 
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Kg” (18h) EHC? = 7,5mg.Kg” (0,5mg.Kg” (Oh), 2mg.Kg"! (9h), 5mg.Kg” (18h)) y OVAP = 
0,25ml.Kg” (Única aplicación); Grupo control (machos) = 0,5mlI de solución salina fisiológica 
(Única inyección aplicada con la última dosis de los animales tratados). 

Para la obtención de los gametos (ovocitos y espermatozoides), los animales 
previamente fueron tranquilizados mediante inmersión en una solución con 2-Fenoxietanol 
(J.T. Baker, USA) a una dosis entre 300-400ppm, hasta observar la pérdida del eje de nado. 
Para el caso de los machos, las variables a evaluar fueron volumen seminal (cm), movilidad 
masal (%), tiempo de activación (seg) y viabilidad espermática (%), osmolaridad seminal 
(mOsm/Kg) y concentración espermática (sptzx10%yL). Entre tanto, a las hembras se les 
evaluaron variables tales como diâmetro ovocitario pre y post inducción (um), fecundidad 
absoluta (fovocitos/hembra) y relativa (fovocitos/Kg). 

Monitoreo de la calidad del agua y de condiciones ambientales: Parámeiros 
físico-químicos del agua como temperatura (ºC), oxígeno disuelto (mg/L), saturación de 
oxígeno (%), conductividad eléctrica (us/cm), pH y sólidos disueltos (mg/L) fueron medidos 
mensualmente durante los muestreos realizados. Para la medición de la temperatura, 
oxígeno disuelto, saturación de oxígeno, conductividad eléctrica y sólidos disueltos se 
utilizó una sonda multi-paramétrica YSI (YSI Professional Plus Ohaus, USA). El pH fue 
medido mediante la utilización de un pH-metro (Waterproof Hanna, Mauritius). 

Así mismo, variables hidrometeorológicas como temperatura ambiental (ºC), brillo 
solar horaria (horas), nubosidad horaria (octas) y precipitación (mms) fueron consultadas 
mensualmente (Fuente IDEAM) con el fin de evaluar la posible influencia de las mismas 
sobre el comportamiento reproductivo de los peces (hembras y machos). 

Análisis estadístico: Los datos fueron expresados como media + error estândar 
de la media (SEM). La información obtenida fue procesada mediante análisis de varianza 
(ANOVA), verificando con anterioridad los supuestos de normalidad (prueba de Kolmogorov- 
Simirnov). Para comparar las medias entre muestreos y estaciones piscícolas, se utilizó 
prueba de Tukey (una vía) y t-student para los casos donde solo evaluábamos dos variables. 
En todos los casos p<0,05 fue utilizado como criterio estadístico para revelar diferencias 
significantes. La totalidad de los procedimientos estadísticos fueron realizados mediante el 
Software GraphPad Prism 5,01 para Windows (GraphPad Software, USA). 


RESULTADOS 


Variación del crecimiento y condición corporal de machos y hembras desde 
los 12 hasta los 18 meses de edad: Dentro de la variable peso corporal (Kg) (Fig. 1), en 
las fincas 1 y 3 no se observaron diferencias significativas entre muestreos al interior de 
cada finca (p>0,05), presentândose un peso inicial (Pi) de 1,62+0,1Kg y peso final (Pf) de 
1,75+0,1Kg para la finca 1 y de 1,60+0,1Kg de Piy 1,90+0,1Kg de Pfpara la finca 3. Entre 
tanto, en las fincas 2, 4 y 5 los valores difirieron estadísticamente entre los muestreos 
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al interior de cada finca (p<0,05). Para la finca 2, el Pi fue de 1,42+0,1Kg y el Pf de 
1,83+0,1Kg; en la finca 4, el Pifue de 1,43+0,07Kg y el Pfde 1,75+0,1Kg; mientras que los 
datos obtenidos en la finca 5 mostraron un Pi de 1,70+0,1Kg y Pfde 2,05+0,1Kg. Como se 
observa en la figura 1, en todas las fincas se presentó ganancia de peso corporal (GP = 
peso final — peso inicial. 
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Fig 1. Peso corporal (Kg) de individuos de Yaque (L. marmoratus) mantenidos en estanques de tierra 
durante los muestreos mensuales realizados. Los valores son mostrados como media + SEM. 


Fig 1. Body weight (Kg) of Yaque individuals (L. marmoratus) kept in earthen ponds during the monthly 
samplings carried out. The values are shown as mean + 